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_wertproblem bestimmt: Es sei w die wohldefinierte reelle und stetige Funktion, ~ 


Uber cin simultanes Diffetenzenverfahren zur Abschawang | 


der Torsionssteifigkeit und det Kapazitat nach beiden Seiten 


Von JOsEPH HERscH, ALBERT PFLUGER und ANDREAS SCHOPF, Ziirich?) 


Die Ausgangspunkte dieser Arbeit sind einerseits eine von Prof. G. POLYA 


im Wintersemester 1953/54 an der ETH. gehaltene Vorlesung iiber Variations- - 


methoden bei physikalisch-technischen Aufgaben, andrerseits ein Artikel [5] a) 


von Prof. A. PFLUGER tiber das Typenproblem einer Riemannschen Flache. 


§1. Charakterisierung der Torsionssteifigkeit als Maximum 
und als Minimum 


1. 1. Wir betrachten das Problem der Torsion eines prismatischen Stabes, 
dessen Querschnitt wir zunachst als einfachzusammenhiangend voraussetzen. 


Zwischen der spezifischen Torsionssteifigkeit P, dem Schubmodul G, dem_ 


Drehmoment M und dem spezifischen Verdrehwinkel # besteht die Beziehung 


G%P=M. Die Torsionssteifigkeit P hangt nur von der geometrischen 


Gestalt des Querschnittes Q ab und ist bekanntlich [12] durch folgendes Rand- 


welche in Q der Gleichung 


Au = Une + Uyy = —2 (tarda t ies 
geniigt und auf dem Rande J" verschwindet. Dann ist 


Pie ai dx dy mba udxdy (Dimension cm‘). (1% 2) 


Gleichzeitig liefert uv die Spannungsfunktion ® = G #u, das heisst t= ice D,) a8 


ist der Schubspannungsvektor. 
Beispiel: Q ist ein Kreis vom Radius 7. Dann ist 


w= > (9%), p25, T=G#B(—y, x). 


ay Mathematisches Seminar der ETH, 
th Die Ziffern in eckigen Klammern verweisen auf das Literaturverzeichnis, S. 112. 
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Wir setzen zur Abkiirzung : [| 1. y) dx dy = f, [vr + f?) dx dy = Dif), 

‘0 Q 
|[e@ dx dy = (v,#), (v, 0) = N(®), wo f eine stetige bzw. stiickweise stetig” 
Q 


differenzierbare Funktion, v und w stiickweise stetige Vektorfelder sind. 
Die Torsionssteifigkeit gentigt zwei bemerkenswerten Variationsprinzipien [8]. 


von CASTIGLIANO. Es beruht auf der Variation der Lésungsfunktion u (bzw. der 
Spannungsfunktion gy) unter Festhalten der Randbedingung. Es sei also p 
eine beliebige stetige und stiickweise stetig differenzierbare Funktion auf dem 
Querschnitt Q, die am Rande I" verschwindet. Wir bilden die von p abhangige 


Grosse te 
Py (fp) = 4p — Dip) 


1. 2. Das Dirichletsche Prinzip oder, in mechanischer Deutung, das Prinzip 
: 


und finden mit Hilfe der Greenschen Formel, indem wir f = u + A setzen und 
(1.1. 1) und (1. 1. 2) beachten, dass P; (p) = P — D(h) und daher 


Pr (fp) SP (1.2. 


: 
E 
| 
ist. Um also eine untere Schranke von P zu erhalten, braucht man nur mit 
einer Funktion # zu «probieren», das heisst die zugehérige Grosse Py (pf) zu — 
berechnen. Deshalb nennen wir # eine «Probierfunktion», Um eine méglichst 
gute Abschatzung zu bekommen, wird man natiirlich eine solche Probierfunk- 
tion wahlen, deren P; (p) méglichst gross wird. Gleichheit tritt in (1. 2. 1) dann 
und nur dann ein, wenn # schon die Lésungsfunktion w ist. ; 
Fiir jede reelle Zahl A (+ 0) ist mit # auch A p eine Probierfunktion. Hat 
man also schon ein f gewahlt, so wird man iiber den Parameter / noch so ver- — 
fiigen, dass P7 (Ap) = 412 6 — A? D(p) méglichst gross wird. Dies ist der Fall, — 
wenn A = 2 $/D(f) ist, und das Maximum betragt 4 62/D(p). Damit haben wir 
folgendes Resultat (P6rya und SZEG@): 
Fiir alle auf Q stetigen und stiickweise stetig differenzierbaren Funktionen p, — 
die am Rande I’ verschwinden (Probierfunktionen), ist | 
1 


P-(p) = 4 $9/D(p) SP, oder Max P-(p) = P. (1.2; 2) a 
P : 
Das Gleichheitszeichen gilt hier bzw. das Maximum wird erreicht fiir alle Viel- | 


fachen du (A + 0) der Lésungsfunktion u. Denn es ist P-(p) homogen vom Grade — 
Null in p, das heisst P-(A p) = P-(p). 


1. 3. Das Thomsonsche Prinzip oder, in mechanischer Deutung, das Prinzip 
der minimalen potentiellen Energie. Es beruht auf der Variation des Vektor- _ 
feldes grad w bzw. der Dehnung unter Festhalten der Bedingung (1. 1.1) in 
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der Form div v = —2.Der Querschnitt Q sei durch endlich viele glatte Bogen I’, 

in Teilgebiete G; unterteilt, und das Vektorfeld v= (vz, Uy) geniige folgenden 

Bedingungen: 

1. v ist in jedem Teilgebiet G; stetig und stetig differenzierbar und erfiillt dort 
die Gleichung div v = —2. 

2. Langs der Bégen J’; ist die Normalkomponente von 7 stetig, jedoch nicht 
notwendig die Tangentialkomponente. 

‘Ein solches v nennen wir Versuchsfeld. 
Fiir die von v abhangige Grésse 


P+(i) = [| [ (v2 + v2) dx dy = NG) 
6 


finden wir mit h = 7 — grad w unter Beachtung von (1. 1. 2): 


=> =e 


P+(v) = P+ N(h) + 2(h, grad u) . 


Wegen div h=0 ist h grad u = div (u h) und daher das letzte Glied rechts 
gleich 


[[ div wit ax ay = X | [ div (ui) dvdy = XY p uhaas, 


Q y ] RdG; 


wo in bezug auf das jeweilige Gebiet G; im positiven Sinne herumintegriert 
‘wird und n jeweils die aussere Normale beziiglich G; bedeutet. Wegen der 
‘Stetigkeit der Normalkomponenten von h langs der J’; und der Randbedingung 
‘u = 0 auf J” verschwindet die letzte Summe. Daher ist P+(v)=P+N (h), und 
wir haben das Resultat: 

 Fiir alle Versuchsfelder v ist 


P+(v) =P oder Min P+(0) =P. (ies. 0 
(v) 


Das Gleichhertszerchen gilt bzw. das Minimum wird erreicht nur fiir v = grad u. 
Jedes Versuchsfeld v liefert in der Grésse P+(v) eine obere Schranke fiir P. 
Um eine gute Abschatzung zu erhalten, wird man Pt(v) méglichst klein zu 


machen suchen. g 
Hat man eine Probierfunktion # und ein Versuchsfeld v gew&hlt und die 


zugehorigen Gréssen P-(p) und P+(v) berechnet, so gilt 
| P-(6) < P< P*(). 


Ist die Differenz P+(v) — P-(f) klein, so hat man eine gute Approximation fiir 
P. Danach beurteilt man die Giite der Wahl von und v. 


a 
- 
7 


Re Cre ree De 
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9 

1.4. Es sei noch kurz auf den Fall eines mehrfach zusammenhangenden 
Querschnittes Q hingewiesen [9]. Man erhalt Q, indem man in einem einfach- 
zusammenhingenden Gebiet G (mit dem Rand J’) einige «Locher» G, mit den’ 
Randern y; (i = 1,2,...,2) anbringt, es ist Q = Gs a G,. F, sei der Flachen- 


inhalt des «Loches» G,. Die Lésungsfunktion u ist in G stetig, genligt in Q 
wieder der Poissonschen Differentialgleichung Au = —2, verschwindet auf ry 
ist auf den Lochern G, und deren Randern y; je konstant und erfiillt die Bedin- 
gungen f du/dn ds = 2F, (i=1,2,...,), wobei im negativen Sinne um G; 
herum integriert wird und 7 die dussere Normale beziiglich Q bedeutet. Dann ist’ 


P=2[fudxdy — ff grad? u dv dy. ’ 
¢ Q | 


P geniigt wiederum den Variationsprinzipien 1. 2 und 1. 3, sofern einerseits 
die in ganz G definierte Probierfunktion # auf den G; je konstant ist und das 


Integral | if p dx dy = p iiber ganz G erstreckt wird; anderseits das in Q defi- 


nierte Versuchsfeld v die zusadtzlichen Bedingungen f vnds=2F,;,1+=1, 2, 
Se inert lt Yi 


. 
§2. Das Differenzenverfahren fiir die Abschatzung der Torsions- : 


steifigkeit nach unten 


Zur Konstruktion geeigneter Probierfunktionen # fiir das Variationspro- 
blem 1. 2 verwendet P6Lya [7] eine stiickweise lineare bzw. bilineare Funktion. 
Hier wird nur die letztere, etwas schneller konvergierende Methode angegeben. 
Wir setzen voraus, dass der Querschnitt Q von achsenparallelen Linien 
begrenzt und in lauter kongruente Quadrate Q; von der Seitenlange h (Ma- 
schenweite) eingeteilt sei. Hat der Querschnitt nicht diese Form, so wird er 
durch eine solche von innen approximiert. Wir schreiben jeder Innenecke des 
Quadratnetzes einen beliebigen Wert #; und jeder Randecke den Wert 0 zu. 
In jedem Quadrat Q, und dessen Rand ist p diejenige bilineare Funktion 
a+bx+cy+dxy, die in den vier Ecken die vorgeschriebenen Werte p; 
annimmt. Haben fiir ein Quadrat Q) die Ecken ¢,, é9, és, ¢, zum Beispiel die 
Koordinaten (0, 0), (h, 0), (4, 2) und (0, A), so ist in Q, 


p=ali-F)(l-F)+ag(i-p+az- 244 (1-2) 2 
Wegen der Linearitat auf allen Kanten stimmen die Bilinearfunktionen zweier 
benachbarter Quadrate langs der gemeinsamen Seite iiberein; p ist somit stetig 
und stiickweise stetig differenzierbar in Q und Null auf I’, ist also eine Probier- 
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funktion im Sinne von 1. 2. Eine einfache Rechnune ergibt 


p=h? S* d; 


DO) = 3 [Xbb + Y 1), 


D 


und 


wo auf der rechten Seite der letzten Gleichung die erste Summe iiber alle 


Quadratseiten und die zweite tiber alle Diagonalen zu erstrecken ist. Gemiass 
(1. 2.1) ist dann 


4m Sa 5([X-pt+ Lb d] <P, ey 


Ss 


_ und es entsteht die algebraische Aufgabe, die Werte fy, py, .-. , Py, So zu wahlen, 
_ dass der linksstehende Ausdruck méglichst gross wird. 


Diese Aufgabe ist ein algebraisches Analogon des Variationsproblems (1. 2). 
Die Lésungswerte w,, #3, ..., %, werden durch folgendes (zum Beispiel durch 
Differentiation sofort erhaltenes) Gleichungssystem bestimmt: Seien E,, eine 
beltebige innere Ecke der Unterteilung, E,, E,,, ..., E;, die acht «benachbarten» 


Ecken, dann soll sein 


= 2 
8 Uz, “a U; — U;, pa U;, baal Uj, Te U;. = Uy, od Uy, Sac U,, ad 6 h . (2. 2) 


Bei x Innenecken haben wir im allgemeinen ein System von » Gleichungen 
in x Unbekannten, die sich aber bei Symmetrie stark reduzieren kénnen. Das 
System hat stets eine und nur eine Lésung: denn das entsprechende homogene 
System besitzt (Mittelwerteigenschaft, daher Maximumprinzip) nur die triviale 
Lésung wu, = ut, = --- = up, = 0. Die erhaltenen Werte p; = u, werden in (2.1) 


: eingesetzt und liefern die fiir diese Maschenweite bestmégliche untere Abschat- 


zung. Fiir diese Werte lasst sich aber das linke Glied von (2. 1) unter Bentitzung 


_ des Gleichungssystems (2. 2) umformen: 


n 8 n 
> u- uy)? + Dd (us — 1)? = Dy Mi, (s m,— 3%] = 6h? d” u,; 
D e=1 a 1 


S =a 
wir erhalten also schliesslich 


an Su, < P. (2. 3) 
! . 


§ 3. Das Differenzenverfahren zur Abschiatzung der Torsionssteifigkeit 
nach oben 


3. 1. Der Konstruktion eines geeigneten Versuchsfeldes v fiir das Variations- 


| problem 1. 3 legen wir dasselbe Quadratnetz wie in § 2 mit der Maschenweite h 


Pee ee Sev 


ay 


zugrunde. Hat der Querschnitt nicht diese Form, so wird er von aussen durch 
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eine solche approximiert. Die horizontalen Seiten orientieren wir nach rechts, 
die vertikalen nach oben. 


In je einem beliebigen Quadrat 0, (i =1,2,..., N) betrachten wir nun} 
wirbelfreie Vektorfelder v;, die langs jeder Seite je eine konstante Normalkom- 
ponente besitzen und im Innern die Gleichung div v; = —2 erfiillen. Den Fluss” 


der Felder v; durch die Seite k;, und zwar von der Nachbarzelle, zu der kh; 
negativ orientiert ist, zur zweiten Nachbarzelle, bezeichnen wir mit h S;. Sind 
die Zahlen S,; (Normalkomponenten) fiir die Kanten h;,,, R;,, h;,, Ri, des Qua 
drates Q; gegeben (der Koordinatenanfang liege in der Ecke links unten), so. 
hat v; = (vz, v,) die Gestalt 
, 


%9=-S,(1-F)-S, 5 w= S,(1-F) + Suz | 


h Pe | 
Wegen der Gleichung div v; = — 2 in Q; kénnen die S;, nicht beliebig sein, son- 
dern miissen die Saacecee | 
eke! oaks + ya 5 Si,=—2h (3.: 16 1) )j 

erfiillen. 


Legen wir nun allgemein den Seiten k; solche Werte S; bei, dass die Bele 
gung (3. 1. 1) fiir jedes Quadrat Q, erfiillt ist, so existiert in jedem Q; ein Vektor-_ 
feld der angegebenen Art. Sind Q; und Q, die Nachbarzellen der Kante k;, so. 
haben die Felder v,; und v, langs k; dieselbe Normalkomponente. Die Felder 9; 
definieren also zusammen in Q ein Vektorfeld v, das den Bedingungen in 1. 3 
gentigt, also ein Versuchsfeld ist. 


3. 2. Eine kleine Rechnung ergibt : 
ao 2 2 2 2 
I | BF do dy ahs (Shop SPREISR AGP aE Gules inane 
Qj 


und wegen : 
S; Si = (Si + SF — (ee 5.12) 


und eines entsprechenden Ausdrucks fiir S, S, 


ff decay = (558 4, SPE Ge ie eee Ae 
Q; : 


Beachtet man, dass jede Innenkante an zwei und jede Randkante an ein 
Quadrat grenzt, und setzt man e,= 1 bzw. 1/2, je nachdem k; Innen- oder 
Randkante ist, so liefert die Summe dieser Betrage der einzelaen Quadrate 


P+(0) = N@) = 9 ey eS ES S/)a) (3. 2.1) 
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Dabei wird die erste Summe tiber alle Kanten wnd die zweite Summe iiber alle 
Gegenseitenpaare der Quadrate erstreckt. — Natiirlich braucht man aber bei 
krummlinigen J’ die Integrale nur innerhalb G zu erstrecken. 

Die Wirbelfreiheit von v in jedem™Quadrat hat tibrigens zur Folge, dass 


_ unter allen Versuchsfeldern mit den Normalkomponenten S; lings der k; das 


eben konstruierte Feld die kleinste Norm N hat. Um fiir das gegebene Quadrat- 


‘netz mit der Methode des stiickweise linearen Vektorfeldes ein giinstiges Ver- 


suchsfeld zu bekommen, bleibt also noch die Aufgabe, den Ausdruck (B.2e 2} 
unter den Nebenbedingungen (3.1.1) zu minimalisieren. Eine rechnerische 
Durchfiihrung dieser Aufgabe ist bei kleiner Kantenzahl nicht allzu mithsam: 
sie fiihrt aber bald zu einem recht komplizierten Gleichungssystem, da die 
Unbekanntenzahl (im allgemeinen) gleich der Kantenzahl ist. (Es sei jedoch 
bemerkt, dass diese Aufgabe mit einer ahnlichen Methode wie der hier zu 
beschreibenden exakt behandelt werden kann.) 

Es ist also bei feiner Rasterung zweckmassiger, eine etwas weniger scharfe 
Extremalbedingung zu stellen, dafiir aber ein einfacheres Gleichungssystem zu 
erhalten. Dazu fihrt die folgende Bemerkung: Die Bedingungen (3. 1. 1) 
fiihren zu 

(SS. )24Se> 5,7 2a 
und [vgl. (3. 2. 1)] 


ss Mit SR 
em (5 - ak. a 3 


wenn N die Zahl der Quadrate, F = N h? den Flacheninhalt des Querschnittes 
Q bezeichnen. Also ist 


P*(s) < Pi @ eal (Dion z): (3, 2:2) 


Bei grober Rasterung kann N(v) wesentlich kleiner sein als der rechtsstehende 


_ Ausdruck, wahrend bei feiner Rasterung dieser Unterschied, wenigstens bei 


den durchgerechneten Beispielen, geringfiigig wird. Es ist daher zweckmassig 
(vor allem auch wegen der Simultanmethode, § 4), den Ausdruck 2 e; S? unter 
den Nebenbedingungen (3. 1. 1) zu minimalisieren. 


3. 3. Diese Minimalaufgabe fiihrt zu einem algebraischen Analogon des 


_ Variationsproblems 1. 3. Bei dem von achsenparallelen Linien begrenzten und 


in lauter kongruente Quadrate Q; (i = 1, 2, ..., N) eingeteilten Querschnitt Q, 
der unserer Untersuchung zugrunde liegt, Dervae ted wir jetzt neben den Ecken 


und den Kanten k; (das sind die Quadratseiten) vor allem auch die Innen- 


? 


flachen, das heisst die Quadrate Q, (Zellen). Die ganze Figur ist ein zweidimen- 
sionaler Zellenkomplex. Die Kanten sind wie in 3.1 orientiert. Eine innere 
Kante k; hat zwei Nachbarzellen. Sie heisst zu einer Nachbarzelle positiv 
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(negativ) orientiert, wenn sie auf der Zelle einen positiven (negativen) Umlauf- 
sinn induziert. Diese zwei Nachbarzellen bezeichnen wir mit Q;+ und Q;-. Eine 
Randkante k; hat dagegen nur eine Nachbarzelle Q;+ oder Q,- in Q. . 

Kehren wir zur Minimalaufgabe am Schluss von 3. 2 zuriick. Wir fassen die . 
Werte S; = G(k,) als Werte einer Kantenfunktion G auf. Als «Norm» von G — 
(beziiglich der Gewichte e;) definieren wir 


NAS) = Df & SF 3.3: 1 
1 


_— 


diese soll unter den Nebenbedingungen (3. 1. 1) minimalisiert werden. 

Zur Lésung fiihren wir Lagrangesche Multiplikatoren 4; = 2 U;/h ein, mit — 
welchen die Gleichungen (3. 1. 1) multipliziert werden sollen. Da jede dieser / 
Gleichungen einem Quadrat Q; entspricht, so liegt es nahe, auch die U; = U(Q;) 
als Werte einer Zellenfunktion UW aufzufassen. Gemadss der Lagrangeschen Regel © 
haben wir nun folgenden Ausdruck zu minimalisieren: 


D6 G+ DA (HS, — S, + Si, + S;, +24) 
iF 0; 


= [SP+ 5 5; -U. + G.)] +O (5 SPE 5; Ur), 


; 1 
innere kj Rand-k; 


| 
. 
: 
a 
wobei wir nunmehr iiber die S; frei verfiigen kénnen. Differentiation nach S; 
ergibt unmittelbar fiir den gesuchten Wert B der Kantenfunktion G: 
. 


fiir innere Kanten kj: V; = —_—+ ; 
fiir Randkanten k;: Vee 44 Uj , : 


das heisst + oder —, je nachdem k; zum inneren Nachbarquadrat Q,+ bzw. Q;- 
positiv oder negativ orientiert ist. | 
Es gilt also fiir alle Kanten k; | 


U., — U 


irae (3. 3. 2) 


wenn wir abmachen, jeder Randkante k; mit innerem Nachbarquadrat Q, ein 
fiktives dusseres Nachbarquadrat Q’ (ebenfalls von der Seitenlange /) anzufiigen, | 
dem wir den Wert U/ = —U; zuordnen. Dies bedeutet, dass U iiber die Rand- : 
kante hinaus durch Spiegelung fortgesetzt wird. Dasselbe gilt, wenn Q; mit mehr 
als einer Kante an die Aussenflache grenzt. Der Fall, wo ein spied Quadrat > 
mit mehreren inneren Q; Seiten gemeinsam hat (einspringende Ecke) macht — 
bei dieser Vorschrift keine Schwierigkeiten: denn der Wert —U, bezieht sich auf 
dieses nur als N achbarn von Q; und wird nur fiir die Randlcants k; verwendet. 


Vol. VII, 1956  Simultanes Differenzenverfahren zur Abs chatzung der Torsionssteifigkeit 97 


3. 4. Eine lésende Kantenfunktion % fiir das Variationsproblem 3. 3 muss 
also auf die erwahnte, zur Gradientbildung analoge Weise von einer Zellen- 
funktion UW ableitbar sein. Damit aber B8 den Bedingungen (3. 1. 1) geniigt, muss 
YU die eindeutig bestimmte Lésung des folgenden Gleichungssystems sein: 

A eee Uy Ce 2h? (3. 4. 1) 


Aus diesem System von N Gleichungen mit den Unbekannten U,, U,, ... , Uy 
wird dre Zellenfunktion WU berechnet, und aus ihr leitet man gemdiss (3. 3. 2) die 
gesuchte Kantenfunktion & ab. 

Bilden nun die U; die Lésung des Gleichungssystems (3. 4. 1), so lasst sich 
gemass (3. 3.2) der Ausdruck (3.3. 1) folgendermassen umformen: 


h? N(®) = 3 ¢ (Ys — U,-)? 
Rj 


= >) U,[4U,-U,-U,—-U,-UJ=2" YU. 
Qj, 2; C0 


Gemiass (3. 2. 2) haben wir also 


Bele), U,- =) =P. (3. 4. 2) 
Gee 
Die Bestimmung der U; bildet die eigentliche Rechenarbeit. Ist diese ge- 
leistet, so erhalt man die rechte Seite von (3. 4. 2) denkbar einfach. — An und 
fiir sich ist (3. 4. 1) eine sehr iibliche Differenzengleichung (vgl. zum Beispiel [1], 
S. 358) ; sie ist hier in natiirlicher Weise gekommen und liefert nicht nur einen 
Naherungswert fiir P, sondern die Abschatzung P,’. 
Einen etwas besseren Wert liefert der Ausdruck (3. 2.1), wenn man dort 
fiir die S; die Werte V; einsetzt. Fiir das Quadrat Q;, wird dann 


2 (S,— S))*=2U,—U,—U,)* dew. (2U,-U,- 0), 


auch am Rande, wenn man durch Spiegelung von Q; in der Aussenflache 
fiktive Quadrate Q/ anbringt und dort den Wert —U; vorschreibt. Wir haben 
dann an Stelle von (3. 4. 2) die etwas bessere Ungleichung 


6. 


= ; ty =1 


N N 1 N 
PRE ROS DLA OLS UE)" engl aed 
1 


wenn wir beachten, dass zufolge der Differenzengleichung (3. 4. 1) 


QEGAU Ui Sh A BO Upa ity de)? 


ist. 


ZAMP VII/7 
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Bemerkungen 


1. Die Vorschrift der «Spiegelung der Werte U; am Rande» kann heuristisch 
so erklirt werden: Man denke sich die Lésungsfunktion z = u(x, y) in 1. 1 als 
Fliche im dreidimensionalen Raum dargestellt; die U; h? sind Naherungswerte 


fiir das Volumen [[ dx dy iiber dem Quadrat Q,, die V; sind Naherungen 
Q; 


fiir die Steilheit (Wert der Normalableitung) langs &;. Im Grenzfall h > 0 ist 
am Rande V = -- 2 U/h, was durch die Fiktion eines aussern Quadrates mit 
dem Wert —U erwirkt wird. 

2. Die durch die Minimalaufgabe 3. 3 bedingte Einfiihrung der Zellenfunk- 
tion U{ sowie die Vorschrift der « Spiegelung am Rande» fiir dieselbe liessen sich 
ebensogut mit Hilfe einer algebraischen Betrachtung des konstruierten Zellen- 
komplexes, dann aber ohne Heranziehung Lagrangescher Multiplikatoren be- 
griinden. . 

3. Schliesslich sei erwahnt, dass die Abschatzung (3.4.2), welche aus 
(3. 4. 3) durch Abschwachung hervorgeht, auch direkt mit Hilfe eines in Drei- 
ecken (statt in Quadraten, vgl. 3.1) linearen Versuchsfeldes fir das Varia- 
tionsproblem 1. 3 erhalten werden kann, und zwar dann als genaue Lésung der 
entsprechenden diskreten Minimalaufgabe. 


§ 4. Simultane Methode 


4.1. Die Methode von Porya (§ 2) liefert fir P untere Schranken P-; die 
in §3 dargestellte Methode obere Schranken P+. Jedoch ist der Umstand, 
dass man fiir jede Maschenweite h zwei verschtedene lineare Gleichungssysteme 
(zur Berechnung von P; bzw. P;) zu lésen hat, fir die Anwendungen unan- 
genehm. Andrerseits stellt sich naturgemass die Frage nach der Konvergenz 
von P; und P; fiir h >0 gegen P; die Beantwortung dieser Frage wird 
durch die Verschiedenheit der Werte u; und U; wesentlich erschwert. 

Nun soll ein simultanes Verfahren dargestellt werden, das erlaubt, aus 
einem eimzigen Gleichungssystem Abschatzungen fiir P nach beiden Rich- 
tungen zu bekommen. Dadurch werden die beiden oben genannten Nachteile 
beseitigt ; die Konvergenzfrage wird dann durch die Gleichheit der Werte p; und 


U; prinzipiell erleichtert, indem sie zu einer rein algebraischen (aber keineswegs 
leichten!) Frage wird. 


4.2. P,; wird wieder genau so berechnet, wie das in § 3 erlautert wurde. 
Das Gleichungssystem (3.4.1) bestimmt die U, eindeutig. Dagegen andern 
wir ein bisschen fiir Pj; die Methode von Pérya. Es wird eine andere Eintei- 
lung von G definiert (Figur 1): Jedes Zentrum E ; eines friiheren Quadrates Q; 


age ae 


ta et i tail 


ee 


wird nun zum Eckpunkt von Quadraten der neuen Einteilung; die neuen Qua- | 
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drate haben dieselbe Seitenlange 4 (Maschenweite) wie die alten; am Rande 
bleiben aber jetzt kleinere, naher zu untersuchende Zellen lbrig. 

Unter Verzicht auf die bestméglichen Werte f; {das heisst diejenigen, die 
P~(p) maximal machen] wihlen wir einfach (E;) = U(Q, :), das heisst , = U;. 


Dieser Verzicht erweist sich an Beispielen als numerisch unwesentlich, also 
gerechtfertigt, wenn man die Vereinfachung der Rechnung beachtet. 


Figur 1 


Es treten in der neuen Einteilung folgende vier Arten von Zellen auf: 
I: Im Innern wiederum Quadrate der Kantenlange h. 

Il: Am Rande Rechtecke mit den Seiten / und h/2. 

III: An den ausspringenden Ecken kleine Quadrate der Kantenlange h/2. 

IV: An den einspringenden Ecken winkelférmige Zellen, welche nochmals in 
drei kleine Quadrate der Kantenlange //2 unterteilt seien (in Figur 1 
gestrichelt gezeichnet). 

Unsere Probierfunktion # sei dann folgendermassen stetig und aces 
bilinear definiert: 

1. Am Rande sei f = 0. 

2. In den Punkten E; sei p = P; (p; vorlaufig beliebig). 

3. Innerhalb der einzelnen Zellen I, II, III, und im Fall IV innerhalb der 

einzelnen Teilquadrate sei # bilinear (mit stetigem Anschluss an die vorge- 
schriebenen Eckenwerte). 


4. 3. Die Ausrechnung ergibt dann fiir die Integrale 
=|[? dx dy und D, = |} (62+ p) dx dy 
Z Z 
(2, He=?: aad 
Z Z 


fiir die einzelnen Zellen Z folgende Werte (mit den speziellen Indizes von 
Figur 1): 
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Fall I (fiir das Quadrat mit den Ecken E, OT OPED oy i 


H, = = (p+ bat Pot Ba)» 


ZAMP 


Dy = (P+ P+ PR pi Pi bs bs Ps) — + (Pi Pat Pods + ba Pat Bahr) - 


Fall II (fiir das Rechteck mit den Eckpunkten Fj, E,): 
h? il 4 
yy ae (by + Pe) oi Pu = - (5 pi + 5 p3 + 2 py pe) - 


Fall III (fiir das kleine Quadrat mit dem Eckpunkt F)): 


h? 
Ay = GE Pi» 


Fall IV (fiir die winkelférmige Zelle mit den Eckpunkten E,, E,, E,): 


2 42 
Dy = 3 pi: 


h2 
Eye Se (25+ 3 p6+ 2 dy), 


Dyy == (SPR + 5 PR + 5 £7 + bs Pe + Ba Pr — Ps br) - 


Bevor wir aufsummieren, wollen wir noch im Hinblick auf spater notwen- — 
dige Umformungen die Ausdriicke fiir Dz in eine Differenz 


DiI 


zerlegen, und zwar: 


D; = pi + po + £3 + pf — bi be — fobs — Psha— Pili» 


/ 3 ee! ane 
Dy 7 PB Pr Pa» 


Diy = 2h, 


he 


/ 3 2 2 2) 2 
Dy = 7 Pat Pi+ = bi — bs be — Pete - 


1 
Ry = 3 (b1— ba t+ bs — pa)?®, 
2 
Ry == (b= p2)*, 
4 2 
Ry mene Pi» 
1 
Ry as Gi 


(4 25 + £6 +4 £2 — 7 bs by — 7 fe by + Bs Pr) - 


—————— >, — ——_—___ 


(4. 3. 1) 


(4. 3. 2) 
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Ubnigens wiirde man im Falle I bei Verwendung einer auf zwei Teildreiecken 
linearen Probierfunktion an Stelle unserer auf dem ganzen Quadrate bilinearen 
Probierfunktion gerade den Wert Dj an Stelle von D, erhalten. 


Fir den ganzen Querschnitt ergibt sich nun durch Summation der ein- 
zelnen Zellenwerte: 


mart : 2 a 
p= Way Pit 4 Xi Pi +46 Dbt ed by 


oder also 
SL A BP PEED ye ED i) ir ati a wir mt eS 
bea BaD, be — adi tat wD, bs z (bs + & bi) (4. 3. 3) 
und 

Ip) = Do — Rs 
wobel 


Diag, DABS DS OPEC tie eS Or, 
SE R A K 
R =) Ry DP Ry tO Rint 2) Rv: 
T Tl Ii IV 
Dabei sollen die Symbole unter den Summenzeichen diejenigen Klassen von 


Eckpunkten, Kanten oder Zellen unserer Zelleneinteilung angeben, tiber 
welche jeweils summiert werden soll. Die Klassifizierung der Zellen ist dabei 


(4. 3. 4) 


Figur 2 


bereits bekannt; diejenige der Eckpunkte ist in Figur 2 dargestellt. Und zwar 
umfasse die Klasse 
J: Alle ganz im Innern des Querschnittes liegenden Eckpunkte mit vier Nach- 
barn, ausgenommen die Eckpunkte der Klasse E. 
R: Alle beim Rande des Querschnittes liegenden Eckpunkte mit drei Nachbarn. 
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A: Alle bei den ausspringenden Querschnittsecken liegenden Eckpunkte mit 
zwei Nachbarn. 

E: Alle bei den einspringenden Querschnittsecken liegenden Eckpunkte mit 
vier Nachbarn. 

Q=J+R+A+E: Alle im Querschnitt liegenden Eckpunkte (die auf dem 
Rande selbst liegenden Eckpunkte unserer Zelleneinteilung zahlen dabei 
nicht mit; auf ihnen ist ja ohnehin p = 0). 

Und schliesslich umfasse 

K: Alle Kanten (in Figur 2 ebenfalls dargestellt; die zum Rande fihrenden 
Kanten der Lange h/2 zahlen dabei nicht mit). 


4.4, Die Maximalisierung von P-(p) = 4 — D(f) fiihrt auf ein ahnliches 
lineares Gieichungssystem fiir die p; wie (2. 2). Die den Eckpunkten EF, der 
Randklassen R, E, A entsprechenden Gleichungen sind dabei naturgemdass zu 
modifizieren. Auf eine explizite Angabe dieses Systems sei indessen verzichtet, 
denn wir wollen es, wie schon gesagt, gar nicht auflésen, sondern uns mit aoa 
Werten #,; = U; als Naherung begniigen. 

Setzt man nun diese speziellen Werte in (4. 3. 3) und (4. 3. 4) ein, so ge- 
statten die Gleichungen (3.4.1), welche die U; definitionsgemiass erfiillen, 
analog wie in 3.4 vereinfachende Umformungen fiir # und D(f). Zunachst 
liefert Summation aller Gleichungen von (3. 4. 1): 


N 2a 2 Fa Sb Um Uo’, ~ Up Uh 2 U, de 
Q;, C2 7 R A 


oder also 


F= SU +25U, (4.4.1) 


was bei der numerischen Auflésung von (3. 4.1) als willkommene Kontrolle 
dienen kann. Ebenso ergibt die Multiplikation der i-ten Gleichung von (3. 4. 1) 
mit U; und anschliessende Summation die bereits in 3. 4 beniitzte Gleichung: 


2h* D0 UD) Uy U;,, — U,, — U, — U,); 
Q;c0 Q;C0 ; 


dies ergibt 
2h? S'U;,=6 SUP +5 J’ UR +4 Y’UF—-2 UU. (4.4.2) 
of Q A R E,J 8 
Anwendung von (4. 4. 1) auf (4. 3. 3) und von (4. 4. 2) auf (4. 3. 4) liefert: 


2p-"2SU- = SU +s LU aia (4.4. 3) 
Q E A 
und 
D(p) = 21? S'U, R= B, (4. 4. 4) 
5 | 
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und man erhalt damit gemiiss (1. 2. 1) und (1. 2. 2) die unteren Schranken 
P>P-(t)= St sae (4. 4. 5) 
und 
PO > BEG) 2 A Ba: Po (4. 4. 6) 
oder, explizite: 


P, = h8)2 yu, 


Q 


meas 


He Yule ry Fhe Ry ea 
+ 7 = J 


wobei also R die durch (4. 3. 2) und (4. 3. 4) definierte quadratische Form in 
den U; ist. Aus (4. 3. 2) ersieht man, dass R; und Ry positiv semidefinit und 
Ry positiv definit ist. Ryy ist aber indefinit, wie folgende Zerlegung in Qua- 
drate zeigt: 


4 


Bases = (Ut U,- = Us) ge (Ue th) we 


BF 


2 
U2. 


Damit erhalt man durch Unterdriickung positiv semidefiniter Bestandteile von 
R folgende untere Schranke fiir R: 


FOB se 7D U2 — ey UP 
E 
und damit eine noch etwas einfachere, aber auch schlechtere untere Schranke 
for P* 
Pe 2m {2 7 U3 US pore FI 
Q 
4 ‘gre ROO Shas 
hu ae tak le 4 | 
in welcher quadratische Glieder nur noch iiber die Querschnittsecken summiert 


werden. Dabei werden offenbar die iiber A summierten Glieder gegentiber den 
iibrigen Gliedern mit abnehmendem / verschwindend klein werden. 


(4. 4. 8) 


4. 5. Betrachtungen zur Konvergenz. Fiir die Differenz P — P, erhalt man 
gemass (3. 4. 2) und (4. 4. 8) den Wert 


2 1 1 20 4 ; 
BaP lta ua yt aan a oe 


und davon die obere Schranke (der Einfluss der iiber A summierten Glieder 


ist ohnehin unbedeutend): 


er ay 


+5 SU+5 ee (4. 5.1) 


Py — Py <3 F 


a 


104 J. Herscu, A. Prrucer und A. Scuorr ZAMP 


Diese Formel gibt einen gewissen Aufschluss tiber das Konvergenzverhalten 
unserer verschiedenen Schranken 


p< Poe Pe ae er 


bei progressiver Verkleinerung der Maschenweite h. Zu ahnlichen Konvergenz- 
aussagen gelangt G.E. ForsitHE bei einem Differenzenverfahren fir den 
. Grundton einer Membran [2]. Insbesondere weist der Term 


> UP? (4. 5. 2) 
E 


auf den kritischen Einfluss der einspringenden Querschnittsecken hin. Denn : 
bei diesen wird bekanntlich die Normalableitung der exakten Spannungsfunk- 
tion u unendlich, somit ist u(Z,;) auf den Eckpunkten E£; der Klasse E von 
tieferer Ordnung als O(h) [und zwar, wie eine genauere Betrachtung zeigt, in 
unserem Falle rechtwinklig einspringender Ecken von der Ordnung O(h?/*)]. 
Fiir die U; auf denselben Punkten F; ist ein ahnliches Verhalten zu erwarten 
und wird bei numerischen Beispielen auch beobachtet. Vom kritischen Term — 
(4. 5. 2) muss man daher eine tiefere Ordnung als O(h?) annehmen; allerdings 
ist es durchaus méglich, dass er erst bei sehr kleinem / gegentiber dem ersten, 
mit /? konvergierenden Anteil von 6 in (4. 5. 1) gross wird. 

Unser kritischer Term (4. 5. 2) riihrt itbrigens nicht lediglich von der even- 
tuell ungentigend genauen Abschatzung der Differenz P+ — P- durch (4. 5. 1) _ 
her; denn ein ahnliches, konvergenzschwachendes Glied bleibt auch fiir P — P-_ 
allein bestehen, selbst dann, wenn die #; nicht = U,, sondern direkt als Extre-— 
male zu maximalem P-(p) gewahlt werden [10]. 

Die Differenz P~ — P, dagegen konvergiert mindestens so gut wie 62, das 
heisst es gilt : 


|P" =P, | =0(6%). (4.5.3) 


Denn einerseits folgt aus (4. 4. 2) bis (4. 4. 6): 


ide A—B=R+0(h?); B=O(I); ote 


ip uate: 2 


0< Pi — Py =O(h?) — R<6 
somit 
|R| <6 + O0(h2) = 0(6) 
und damit 


E 

1 

. 

Andrerseits folgt aus (3. 4. 2) und (4. 4. ihe | | 
| 

. 

|A — B|=0(8),. also |P- — P| = 063). : 
. 
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§ 5. Beispiel 


Als Querschnitt wurde ein aus fiinf Quadraten der Kantenlange 1 zentral- 
symmetrisch zusammengesetztes Kreuz gewahlt. 

Die fiir h=1, 1/2, 1/3, 1/4 berechneten Resultate und die wichtigsten 
Zwischenresultate sind in der folgenden Tabelle wiedergegeben. 


0 (extra- 
poliert) 


SE, Se ok aac See 78 2,548 2,244 2,114 1,945 
Ma ea ei ee + nal S907 2,194 2,040 1,981 1,905 
| SAS AR EES SP ae Ar re 1,574 1,676 1,734 1,808 
ol a eae 1,229 1,549 1,675 1,837 
ME ed ey. cul ck 1,024 1,414 1,575 1,783 
eee | “5167 2,964 2,430 2,218 
Q 
isd (2U, —U;,—U, —h*)2 | 0,593 0,354 0,205 0,133 
ee ered eh 0,554 0,415 0,341 
oS aga 0,161 0,085 | 0,053 
pay 2 ite aes ee 0,908 0,511 0,345 
3/8 = SU oi, — 14/9 0,2. 0,097 0,051 0,036 
pea (Hs; 0,)3 Siceihihs deie 0,109 0,067 0,042 
me et UE Us) 0 0,018 | 0,022 
eee OP ee hah SE ED 0,276 0,077 0,030 
p | 


Im Diagramm in Figur 3 sind die Werte von P,, Ps. Lae ewes 
Funktion von h? ausgetragen und geradlinig verbunden. Die Strecken zwischen 


h = 1/3 und h = 1/4 sind geradlinig bis h = 0 weitergefiihrt (extrapoliert). Wie 


nach den Ausfiihrungen von 4.5 zu erwarten ist, werden die Verbindungs- 
strecken mit abnehmendem h allmahlich steiler. Fiir die vorliegenden Werte 


- geschieht dies sogar auf monotone (konvexe) Weise. 


een se 


Von G. Lein sind Schranken fiir eine ganze Klasse von kreuzformigen 
Querschnitten auf andere Weise berechnet worden [4], insbesondere fiir unseren 


ae 1,9017 < P< 1,9042. 


Merkwiirdig ist dabei, dass obige Schranken fiir alle kreuzformigen Quer- 
schnitte, welche aus dem vorliegenden durch Translationen der beiden Kreuz- 
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balken hervorgehen, ebenso gelten sollen. Im Diagramm ist der Wert 1,90 
durch eine horizontale strichpunktierte Linie dargestellt. 

Aus dem Diagramm ersieht man auch, dass der extrapolierte Wert von P- 
schlechter ist als derjenige von P,, was durchaus im Einklang mit (4. 5. 3) 


steht. 


Okeke 7 


Figur 3 


Schhiesslich sei noch auf den besonderen numerischen Vorteil hingewiesen, 
welchen die Extrapolation erfahrungsgemass bei solchen Differenzenverfahren 
bietet, wenn man bedenkt, dass die Anzahl der Unbekannten mit A-? und die 
Rechenarbeit mit der dritten Potenz der Anzahl Unbekannten wachst. Aller- 
dings besteht die Hauptarbeit nur noch in der Auflésung eines Systems fiir die 
U,; der Randklassen R, E, A allein, da sich die tibrigen U; der Klasse J in ein- 
facher, durch die Zelleneinteilung selbst nahegelegter Weise eliminieren lassen. 
Aber auch so wachst die Rechenarbeit immer noch mit A-; bei Annahme einer 
Konvergenz mit h® verlangt damit eine Vergrésserung der Genauigkeit um eine 
Dezimalstelle immer noch einen rund 30mal grésseren Rechenaufwand. 
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Fir 4 = 1/4 war in unserem Beispiel ein System mit drei Unbekannten auf- 
zuldsen, wobei nattirlich auch die Symmetrie des Querschnittes zur Erniedri- 
gung der Anzahl der Unbekannten beitrug. Der aus Pyyjg und Py), in h? 
linear extrapolierte Wert lieferte bereits: 


16P* .,—9P* 
Pho = tatlt #718 19047. 


§ 6. Elektrostatische Kapazitat 


6. 1. Im dreidimensionalen Raum beschaftigen uns zwei Arten von Zellen- 
funktionen, namlich drei- und zweidimensionale; diese wollen wir «Wiirfel- 
funktionen» bzw. « Quadratfunktionen» nennen. 

Wir betrachten einen Kondensator, der etwa durch den Zwischenraum G 
zweier geschlossener Flachen (Elektroden) Ey und E, gebildet wird, wobei E, 
die andere Elektrode E, vom Unendlichen trennen mége. 


Bemerkungen 


1. Im Falle GS oo ist unsere Methode wegen der Anwesenheit unendlich 
vieler Unbekannten ungiinstig. 

2. Eventuelle adiabatische Teile des Randes von G liessen sich dagegen 
ohne Schwierigkeit beriicksichtigen. 

Es sei u die in G harmonische Funktion mit den Randwerten 0 auf Ey und 
1 auf E,. Dann ist die elektrostatische Kapazitat C des Kondensators gegeben 
durch 


42C =D(u) = [| [grag u dt = tp ou OF (eat 
uC 


Ey 


(dx = Volumenelement, dF = Flachenelement, n = dussere Normale). 


6.2. Ist p eine beliebige stetige und stiickweise stetig differenzierbare 
Funktion («Probierfunktion») in G mit den Randwerten 0 auf E, und 1 auf Fy, 
so gilt nach dem Dirichletschen Prinzip C S C+(p), wo 


4 C+(p) = D(p) (6. 2. 1) 


ist. Setzen wir ndmlich = u+h, so wird 4aCt(p) =4aC0+ D(h); denn 


If gradu gradh dr = (hh 5% dF — |[[h du dr—0. 
G G 


E, +E, 


Nur fiir p= u ist Ct+(p) = C. 
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6.3. G sei durch endlich viele glatte Flachenstiicke J’; in Teilgebiete G; 
unterteilt und v = (v,, v,, ¥,) ein in allen G, definiertes Vektorfeld, welches den 
folgenden Bedingungen gentgt: ; 
1. v ist in jedem Teilgebiet G; stetig, stetig differenzierbar und quellenfrei 

(div v = 0). . ae 
2. Langs der Flachenstiicke J’; ist die Normalkomponente von » stetig, jedoch 

nicht notwendig die Tangentialkomponente. (Mit anderen Worten, es ver- 

schwindet dort die Flachendivergenz: Div v = 0.) 


3, That durch E, den Fluss {p rhea 
EY - 
Ein solches v nennen wir jetzt Versuchsfeld. Es gilt C 2 C~(v), wo 


[4.00 C-(0) = [{f* dv= NG) (6:3: 0a 


ist. (Etwas allgemeinere Fassung des Thomsonschen Prinztps.) 
Setzen wir namlich v = [4 C]-! grad u + h, so ist 


[40 C-(v) 4 = [42 C}-4+ N(h); 


Hel h gradu dx = X /f h gradu dr = SY fp wh nar 
(ex J G; J 1 


denn 


wegen divi = 0; nun verschwindet aber die letzte Summe iiberall: iiber den J’, 


wegen der Stetigkeit der Normalkomponente von v (also auch von h); iiber — 


E, und E,, da u dort konstant ist und der Fluss von h verschwindet. 


6.4. Unsere Methode zur unteren Abschatzung der Kapazitaét ist ganz 
analog der von 3.4 fiir die obere Abschatzung der Torsionssteifigkeit. Das 


Raumgebiet G wird in Wiirfel der Kantenlange # zerlegt (bzw. von aussen — 


approximiert). — Der Kiirze halber beschranken wir uns hier auf Quadrat- — 


funktionen, welche von Wiirfelfunktionen ableitbar sind. Diese Einschrankung 


sowie die Vorschrift der «Spiegelung am Rande>» liessen sich analog wie in — 


3. 2 und 3. 3 rechtfertigen. 


Jedem Wiirfel W; wird eine reelle Zahl U; zugeordnet; diese Zahlen sollen — 


folgendes Gleichungssystem befriedigen: Fiir jeden Wiirfel W,, und die 6 be- 


nachbarten W,,, W;,, W;,, W;,, W;,. Wj, soll die Gleichung gelten: 


2? 


OO Ur Ue U; U;, U;,=0. (6. 47 


Wir halten die Vorschrift der «Spiegelung am Rande» (vgl. 3. 3) in folgender 


Form aufrecht: Gehért ein Seitenquadrat des Wiirfels W, der Randflache E, 


an, so soll dem dusseren Wiirfel, der mit W, dieses Quadrat gemeinsam hat, © 


} 


vo re 


aA ey ey 
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der Wert —U; zugeordnet werden. — Dabei ist wieder gemeint, dass der Wert 
—U,; dem ausseren Wiirfel als Nachbarwiirfel von W,; zugehort; ist etwa auch 
W, mit ihm benachbart, so bekommt er auch, Wistselich W,, den Wert —U,. - 
Gehort aber ein Seitenquadrat von W, der-Randjfldche E, an, so erhalte der benach- 
barte dussere « Bildwiirfel» den Wert 2A—U,, wobei A eine spiter zu be- 
stimmende Konstante ist. - Wir ordnen noch dem gemeinsamen Seitenquadrat 
Qi(x = a) zweier Wiirfel W; (x <a) und W, (x> a) den Wert V; = (U, — U;) |h 
zu, ebenso fiir y bzw. z statt x. 

Wir konstruieren nun in jedem Wiirfel W,, (0S 4, y,2 Sh; die Seiten- 
quadrate Q;: x = 0 bzw. x=h, y=0, y=h, z= 0, zg =h bilden den Durch- 
schnitt von W;, mit W;, bzw. mit W;,, W;,, W;,, W;,, W;,) ein quellenfreies Vektor- 
feld v = (v,, vy, "”,) mit de Rieceiconponeate v, = V; auf dem Quadrat Q;: 


- (ee 
= AS. NaS. OR bg 
ra bs: 


ae 


h 
U. 
v, a Jo 7 Js (1 7 + M4 : Jo +; 


Man verifiziert sofort: div v= —(6 U;,, — U;, — U;,—--- — U,,)/h® = 0 wegen des 
Gleichungssystems (6.4.1), dem die U; geniigen miissen. — Auf allen Seiten- 
quadraten der Wiirfel hat v eine stetige Normalkomponente; damit v ein 
Versuchsfeld ist, soll noch der Fluss durch £, nach aussen = 1 sein; dieser ist 
aber wegen der Quellenfreiheit von v in G gleich dem Fluss durch Ey nach 
innen; das heisst wir bekommen die zusatzliche Bedingung: 


h >” (2A4—2U)=1 oder h 2U,=1. 


0; CE NW; 0; CE NW; 


(6. 4. 2) 


Sei N die Anzahl der Wiirfel in G; die N Gleichungen (6. 4. 1) und eine der 
Gleichungen (6.4.2) bilden ein System zur Bestimmung der N +1 Grdssen 
U,,..., Uy und A. Mit den so bestimmten Grdssen erhalt man also ein zu- 
cae Versuchsfeld fiir das Thomsonsche Prinzip (6. 3. 1). 

Die eindeutige Lésbarkeit des Gleichungssystems wird dadurch sicherge- 
stellt, dass das entsprechende homogene Ss [das heisst mit = 0 statt = 1 
in (6. 4. 2)] nur die triviale Lésung U, = --- = Uy = A =O besitzt: denn ein 
maximaler Wert + 0 kann keiner dieser N + 1 Gréssen zukommen. 


6. 5. Fiir den oben betrachteten Wiirfel Wj, berechnet man leicht 


if Wee h 
f[fwanmz dG, Ci ae 
W; ge 
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wobel 
A, = (2U, —U,— U,)*+ 2Y, — U, — O84 (20 Us U;)?. (6.5.1) 
Summieren wir iiber alle Wiirfel W; C G und schreiben wir tiberall 


(U, — U,)? = U; (U;, — Uy) + U, (U, — Y) 
so erhalten wir dank (6. 4. 1) 
aN > " h3 . he 
[4.2 C-(v)]2 = N(v) = W8 DS Wp > > Wie i prc 
oi 


innere Q; Q;C EVE: 
h 
0; CALNW; W; 
das heisst bei Beriicksichtigung von (6. 4. 2) 
[42 C-@)} 1 = A - He DD Aas . (6.5.2) 
; ; 
oder etwas gréber 
C2G= aaa (6. 5. 3) 


Fiir letztere Abschatzung geniigt es also, aus dem Gleichungssystem nur den 
Wert von A zu berechnen. 


6.6. Simultane Methode. Wir verfahren ganz analog wie im § 4. Statt der 
bis jetzt gebrauchten Wiirfeleinteilung fiihren wir fiir die untere Abschatzung 
der Kapazitat eine andere ein, indem jedes Zentrum Z; eines alten Wiirfels W; 
zum Eckpunkt von Wiirfeln der neuen Unterteilung wird. Es bleiben noch am 
Rande Zellen iibrig, die aus Wiirfeln der Kantenlange 4/2 zusammengesetzt 
sind. — Nun ordnen wir jedem Punkt Z; den Wert 4; = U; zu, der in 6.4 dem 
Wirfel W; zukam; auf der Randflache Ey setzen wir tiberall = 0; auf der 
Randfliche E,: = A. In jedem neuen Wiirfel konstruieren wir sodann die 
trilineare, das heisst in jeder der drei Variablen x, y und z lineare Funktion 4, 
die an den 8 Ecken Z; die vorgeschriebenen Werte p; annimmt. Da # in allen 
Seitenquadraten bilinear ist, ist der stetige Anschluss zwischen zwei benach- 
barten Wiirfeln gesichert; so dass p/A als stetige Probierfunktion fiir das 
Dirichletsche Prinzip zugelassen ist. — So bekommt man aus denselben Werten 
p; = U; auch eine obere Abschatzung C+(p/A) fiir die Kapazitat. 


Bemerkung 


j; Diese Verwendung einer stiickweise trilinearen Probierfunktion ist wieder 
eine etwas abgeschwachte Fassung eines (zum § 2 analogen) Verfahrens von 
Porya, das eine Bestimmung der Eckenwerte p; durch die Extremalbedingung 
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C*(p) = min vorsah. Erst unser Verzicht auf diese Bedingung erméglicht in 


relativ einfacher Weise die Gewinnung der Schranken C- und C+ aus einem 
einzigen Gleichungssystem. 


6. 7. Dasselbe Problem der Kapazitat lisst sich auch fiir 2 statt 3 Dimen- 


sionen ebenso lésen (Kapazitait pro Langeneinheit). Das zweidimensionale 


Gebiet G (mit den Randkurven Fy und E,) wird (wie beim Problem der Tor- 
sion, §§ 3 und 4) in Quadrate eingeteilt, denen Werte U; zugeordnet sind. Die 
Gleichungen (6.4.1) werden durch 4 U;, — U;, — U;,— U; — Uj, = 0 ersetzt, 
die Vorschriften von 6.4 iiber die «Spiegelung am Rande» iibertragen sich 
wortlich (« Quadrat» statt «Wiirfel», «Kante» statt «Quadrat» lesen); so wird 
wieder eine Konstante A eingefiihrt. Das Feld v wird wie in 3. 1 konstruiert; 
es ist jetzt quellenfrei. Die beiden aquivalenten Gleichungen (6. 4. 2) lauten 
jetzt 

(2A —2U)=1 bew. as Z2U.=1. (6. 7. 1) 

ky CEQ; ky CEONQ; 


Ist N die Anzahl der Quadrate Q; im Gebiet G, so haben wir wieder ein 
System von N + 1 linearen Gleichungen in den N + 1 Unbekannten U;,, A. 


Man erhalt [42 C-(v)|-1=N(v) =A YY (2A—2U)—1/6 YA, wo nun 
Re CEWNQ; Q; 


Jo 


BA, = (2U, —U, —U,)?+ (2U, —U, —U,)*=2(2U, —U, —U,)3; (6.7.2) 


| wegen (6. 7. 1) haben wir also 


> i 
[4 IU C-(v) au = A = iat A; y (6. ‘Te 3) 
oder etwas groéber wieder 
es 2 1 


Fiir die Simultanmethode wird die stiickweise trilineare Funktion von 6. 6 
durch eine bilineare ersetzt. 


§ 7. Schlussbemerkungen 


Es wurde verschiedentlich auf die Dualitét der Prinzipien von DIRICHLET 


und von THomson hingewiesen, die Abschadtzungen in verschiedenen Rich- 
tungen liefern. Die hier dargestellten Methoden zur Gewinnung oberer und 


unterer Schranken fiir Torsionssteifigkeit und Kapazitat bestehen in einer 


" Algebraisierung des Dirichletschen Prinzips durch Zuhilfenahme einer diskreten 
-Eckenfunktion #; (0-dimensionale Funktion [G. Porya]) und einer Algebraisie- 
rung des Thomsonschen Prinzips durch Verwendung von Zellenfunktionen U; 


Ss sa) Bo 


_ (2- baw. 3dimensionale Funktionen im Falle der Ebene bzw. des Raumes). 
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Zwischen diesen Prinzipien und ihren Algebraisierungen besteht eine sehr 
weitgehende formale Analogie, auf welche hier nur vereinzelt hingewiesen — 
wurde, die aber den mathematischen Hintergrund der angegebenen Methoden 
darstellt. Unter diesem Gesichtspunkt tritt hinter der erwahnten Dualitat der 
beiden Prinzipien jene der algebraischen Topologie eines n-dimensionalen Zel-— 
lenkomplexes hinsichtlich der komplementaren Dimensionen & und  — k in 
Erscheinung. | 
Die Grundgedanken der vorliegenden gemeinsamen Arbeit wurden 1954 an- 
lasslich der Jahresversammlung der Schweiz. Mathematischen Gesellschaft — 
in Altdorf bekanntgegeben [6]. | 
Obere Schranken fiir die Torsionssteifigkeit sowie untere Schranken fiir die 
Kapazitat kann man auch mit Hilfe von Zellenfunktionen ohne Benutzung des — 
Thomsonschen Prinzips gewinnen [3]; die Dualitat der komplementaren Di- | 
mensionen & und n — k tritt auch in Abwesenheit eines Thomsonschen Prinzips 
in Erscheinung. 
Eine nahere Untersuchung tiber Konvergenzfragen, die zum Teil in 4.5 
auftauchten, wird in einer Arbeit von A. ScHopF [10] erscheinen. | 
Fiir viele Untersuchungen iiber Differenzenverfahren wende man sich an 
das Standardbuch von L. CoLLatz [1] und die dort angegebene Literatur. 
Schliesslich sei auf die mit unserem Gedankenkreis verwandte schéne 
Methode von J. L. SyNGE [11] hingewiesen, die auf anschaulichen Betrachtun- 
gen im Funktionenraum beruht. 
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Summary 


We consider the torsion of a prism. Let its cross-section Q be covered by a 
square network (lattice) of mesh size h. PoLya defined a function u, on all lattice 
points by means of a partial difference equation, and then a function p on Q by 
_ bilinear continuation of uw; in every square. To he applied DrricHLEt’s principle 
and obtained thus lower bounds for the torsional rigidity P of Q. — Here we define 
(§ 3) a cell function U;, attributing a real value to each square; using U;, we 
construct a vector field v, to which we can apply THomson’s principle and get 
upper bounds for P.— A simultaneous method for upper and lower bounds is 
then indicated and discussed (§ 4) and an example is given (§ 5). — Evaluations for 
capacity may be calculated in a very similar way (§ 6). 


(Eingegangen: 11. Januar 1955.) 


Numerische Behandlung der Diffusionsgleichung 
mit Konvektionsterm’) 


Von EpuarpD BATSCHELET?) und FRANZ GRUN*), Basel 


1. Einleitung 


In einer kiirzlich erschienenen Arbeit hat M. ABRAMOWITZ [1]*) ein speziel- 
les Problem der Warmekonvektion durch Separation der Variablen und durch 
Reihenentwicklung gelést und gezeigt, wie die Lésung fiir die numerische Aus- 
wertung verwendet werden kann. Wir beschreiben im folgenden bei einem ver- 
wandten Diffusionsproblem ein andersartiges Verfahren, das ebenfalls zu einer 
brauchbaren numerischen Lésung fiihrt. Nur diirfte der Rechenaufwand we- 
sentlich geringer sein. Unser Vorgehen ist itberdies nicht an die Voraussetzung 
einer parabolischen Geschwindigkeitsverteilung gebunden. 

Wir ersetzen zuerst die Differentialgleichung und die Randbedingungen des 
Problems durch ein System von linearen Differenzengleichungen (Abschnitt 3). 


1) Die vorliegende Mitteilung handelt von einem Teil eines Forschungsprojektes «Probleme 
des Gasaustausches in der Lunge» der Arbeitsgemeinschaft K. Bucuer, Pharmakologische An- 
stalt, und F. Grtn, Physikalisch-Chemische Anstalt. Die Untersuchungen werden in verdankens- 
werter Weise subventioniert von der F.-Hoffmann-La-Roche-Stiftung zur Forderung wissen- 
schaftlicher Arbeitsgemeinschaften in der Schweiz. 

2) Mathematische Anstalt der Universitat Basel. 

3) Physikalisch-Chemische Anstalt der Universitat Basel. i 

4) Die Ziffern in eckigen Klammern verweisen auf das Literaturverzeichnis auf Seite 120. 
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Im Gegensatz zu den meisten Anwendungen des Differenzenverfahrens aut 
partielle Differentialgleichungen vom elliptischen Typus ergeben sich hier bei 
der Auflésung zwei besondere Schwierigkeiten: Erstens sind unendlich viele 
Gitterpunkte vorhanden. Zweitens konvergiert die iterative Auflésung der 
Gleichungen nicht. Das sonst so bequeme Relaxationsverfahren ist daher nicht 
anwendbar. 

Umgekehrt wird unsere Aufgabe durch die speziellen Voraussetzungen er- 
leichtert, die wir in Abschnitt 2 iiber das Integrationsgebiet treffen: Es soll die 
Form eines Kreiszylinders besitzen. Damit wird es méglich, einen Ansatz mit 
Exponentialfunktionen zu machen und die Differenzengleichungen ohne Itera- 
tion exakt aufzulésen. Fiir den Fehler, der beim Ubergang von der Differential- 
zu den Differenzengleichungen entsteht, konnen Abschatzungen ermittelt wer- 
den (Abschnitt 6). 


2. Formulierung des Diffusionsproblems — 


Eine Fliissigkeit, die eine Substanz A in der Konzentration u* geldst ent- 
halt, strémt, aus einem Reservoir kommend, durch ein unendlich langes, kreis- 
zylindrisches Rohr vom Radius R. Das Geschwindigkeitsfeld sei stationar und 
zur Rohrachse symmetrisch. Der Geschwindigkeitsvektor v(7) sei nur Funktion 
des Abstandes von der Rohrachse. Wir lassen die Rohrachse mit der x-Achse 
eines Koordinatensystems zusammenfallen (Figur 1). 


Figur 1 


Das Geschwindigkeitsfeld der str6menden Lésung. 


Wir nehmen an, die Rohrwand sei fiir negative x-Werte fiir die geléste 
Substanz undurchlassig. Sie sei aber fiir positive x-Werte aus einem fir die 
Substanz A durchlassigen Material angefertigt. Ausserhalb des Rohres sei die 
Konzentration der Substanz A durchweg gleich Null. Unter diesen Annahmen 
findet ein standiger Transport von Substanz A aus dem Rohr nach aussen 
statt. Die Substanz diffundiert zundchst in der Fliissigkeit mit der Diffusions- 
konstanten D von der Rohrachse gegen die Wand und von da durch die Wand 
mit der Diffusionskonstanten Dy nach aussen. Durch die Strémung wird be- 
standig Substanz nachgeliefert. Gesucht ist die Konzentration u der Substanz A 
in allen Punkten des Rohres, nachdem sich der stationdre Zustand eingestellt hat. 
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In Zylinderkoordinaten lautet die Differentialgleichung fiir die unbekannie 
Funktion u = u(x, 7) 


O70 02u s v(r) Ou 
oxt ort te ae De Oe 9: (1) 


i An der Stelle x = 0 gelte die Randbedingung®) 


“(O, 7) = a* (2) 
) und an der Rohrwand 
Ou 
(Seleep = RM (3) 


_ wo k eine von der Rohrwand und von D abhangige Konstante ist®). Schliesslich 
_ muss gelten 3 
lim u(x, 7) =0. (4) 
#—>0O 
Aus Symmetriegriinden kénnen wir uns in Figur 1 auf den Streifen «ober- 
_halb» der x-Achse beschranken. Lings der x-Achse entsteht so ein neuer 
«Rand», fiir den die Bedingung 


(or ),-0 =° (5) 
gilt. 


3. Ubergang zu den Differenzengleichungen 


Zur Vereinfachung der Darstellung wahlen wir ein ziemlich grobes quadra- 
tisches Gitter (Figur 2). Bei Verwendung eines feineren Gitters kann man in 
gleicher Weise vorgehen. 


Figur 2 
Approximation durch ein Gittergebiet. 


5) Wir sehen vom Randfehler ab, der durch Zuriickdiffundieren der gelésten Substanz in das 


Gebiet negativer x-Werte entsteht. 


8) An der Grenze zwischen Fliissigkeit und Rohr ergibt die Erfiillung der Kontinuitatsbedin- 
gung fiir die Diffusionsstr6me die Bedingung 
u 
(D grad ") Normalkomp. Dy 2 


. 


wo 6 die Dicke der Wand bedeutet. Daraus wird k = Dy/D0. 


; 
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Fiir die Gitterseite wahlen wir den Wert 4 = R/3. Die Funktionswerte - 
u(x, ”) approximieren wir durch die Gitterwerte U, am Rande 7= R, durch” 
u, langs der Parallelen ry = 2 R/3, durch w, langs der Parallelen 7 = R/3 und 
durch W,, lings der x-Achse (m = 1, 2, 3,...). Bei Verwendung von zentralen 
Differenzen lauten die Differenzenapproximationen langs der Geraden r= 2 k/3 
und fir += nh: | 


Ou Un ex Un 1 R Ou U n <> Wy RK 
0% 2 a eee et i + ett ah “A bicep aa ae Bee + R 
“Ox? 07 h 


Darin bezeichnen §,, 85, Rz die spater im Abschnitt 6 angegebenen Restglieder. 
Entsprechende Ausdriicke haben wir langs der Geraden 7 = R/3 anzusetzen.— 
Setzen wir (6) in (1) ein und vernachlassigen wir die Restglieder, so erhalten 
wir die Differenzengleichungen: 


Fung = 2h: 

h Voh h Von ss { 
(1 Es st) n+1 sta 4 U,, st (1 aig 3 | Un—1 =r gid oh “as Un = 0 , (7) 
Pir y= | 

fee Avy, 3 h Up, i 
( oth) nat 7 tn + (1+ 53) Wart > Wa—4,=0. 
Durch (2) sind die Gitterwerte fiir x = 0 vorgegeben: 
UY =% ===; (8) 


Zur Befriedigung der Randbedingungen fiir y= 0 und y = R verwenden wir 
die folgende Dreipunktformel der numerischen Differentiation’) : 


( af 3:7.) it feet hh) = 2h 
ae us ate i ) = fH% + ) + Ry. (9) 


Das Restglied 8, wird in Abschnitt 6 ausgeschrieben werden. Mit Hilfe von 


(9) lauten die den Randbedingungen (3) und (5) entsprechenden Differenzen- 
gleichungen 


—3+2h8) Ut 4, —w,=0, —3W, 44,4. 20, ie 


Zur Auflésung der simultanen Gleichungen (7) und (10) eliminieren wir zunichst 
U,, und W, und erhalten so 


Una t Punt yun_1+6w,=0, EWnry + CWn+ HYWy_1t+0u,=0, (11) 


) Me t} oden der Randappr O% ition verwe I 1 Coruna TZ ho S: 428 le T Ox 4 
Ander e€ x1ma 
_ nde . > ? [ J 
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worin wir zur Abkiirzung 


Z h voy 
l= Atahe > a= 1 <5 =v pe ow) = 4 Y es i! “+ oe > 
| —5l hi 10 i ee 
(ae re) : es Tyee h Oe, hv, pea tt 
4 Be Bete pee, Me phat mag 


gesetzt haben. 


4. Die exakte Auflésung der Differenzengleichungen 


g Die Gleichungen (11) stellen fiir 7 = 1, 2, ... ein System linearer homogener 
Differenzengleichungen mit konstanten Koeffizienten dar. Zur exakten Auf- 
_ lé6sung macht man den Ansatz 


t= AO”, vw, = BO; 
Dies fiihrt zu den Gleichungen 
A(a6?+BO+y)+Bd0=0, AdO+ B(cO?+6O047)=0, 

die nur befriedigt werden kénnen, wenn 
elo ty 60 | 

. vO e084 CO4+n 
=a604+ «06+ 62) 084+ aqnt+bl+yve—8 0) O 

+ (Bntyl)O+yn=90 


(13) 


ist. Die Formel (13) stellt die charakteristische Gleichung des Problems dar. 
Ihre Wurzeln seien mit 0,, O,, Os, O, bezeichnet’). Dann lautet bekanntlich 
die allgemeine Lésung der Gleichungen (11) 


4 ca 
ge ef ns Dis Ba Se p= Oo, 2.2. Je (14) 


Die Koeffizienten A,, B, lassen sich aus den Randbedingungen bestimmen. 
Die Bedingung (4) kann nur erfiillt werden, wenn alle charakteristischen Wur- 
_zeln vom absoluten Betrag => 1 unterdriickt werden, wenn also die zugehérigen 
- Koeffizienten in (14) gleich Null gesetzt werden. Die restlichen Koeffizienten 
ergeben sich bei der Erfiillung von (8). 


8) Wir betrachten hier den Fall mehrfacher Wurzeln nicht, da er im vorliegenden Problem 


ohne Bedeutung ist. 
} 
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5. Zahlenbeispiel 
Wir legen die folgenden Daten zugrunde’®) : 
R=6-10-*cm: D=10-5cm?s-!; k=0,5- 10"coi- #3 


Fiir die Stromung der Flissigkeit im Rohr nehmen wir eine parabolische Ge- 
schwindigkeitsverteilung an: 


v(7) = (1- ee) | mit v = 0,54 cm s-. (15) 
Die Differenzengleichungen (11) lauten dann 
s Zo 
— 2 thn gy — en + 4 thang Pee Wi = Dis 


19 10 29 4 
area Cs 3 Wn 5 Wn-1+ a Un= 0. 


Daraus gewinnt man die charakteristische Gleichung 
22,8 04 + 47,7875 0? — 57,4625 O? — 82,4125 O + 69,6 = 0 7 


mit den Wurzeln @, = 0,9911378..., 0, = 0,8213722..., 0,= —1,69.:.3 
@,= —2,21.... Da die Lésung fiir 7 >co, wie oben erwahnt, nach 0 konver- 
gieren muss, kann sie nur die beiden Wurzeln 9,, @, vom absoluten Betrag < 1 


ie 


W 82 Wb 896 


G64 878 835 828 88 80 b2 794 
100 8 998 G5 954 MY G5 25 Ib 37 899 890 $82 874 BB 


100 339 $87 $§3 $89 $83 978 $it $5 358 951 Hd G8 G8 921 


100 100 0 00 100 996 92 


Bl 35 IE Ibt 954 IE 939 


Figur 3 


Die exakte Losung der Differenzengleichungen. Die Werte stellen die Konzentration des diffun- 
dierenden Stoffes dar. 
®) Das in dieser Mitteilung behandelte Diffusionsproblem stellt sich bei der Betrachtung eines 
Modells, das wir uns machen, um den Gasaustausch aus einer einzelnen Lungenkapillare in den 
umgebenden Alveolarraum theoretisch zu behandeln. Die fiir die Berechnungen im Text verwen- 
deten Zahlenwerte sind die fiir dieses Modell anzunehmenden Daten. Fiir Naheres siehe Publika- 
tionen, die in den Helvetica Physiologica et Pharmacologica Acta erscheinen werden. 
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enthalten. Die Lésung lautet: 
U, = A, O{ + AO), w, = BL Ot + BeO% (n=0,1,2,...). 


Setzt man die Anfangskonzentration-w, = Wy = 100, so gewinnt man fiir die 
Koeffizienten die Gleichungen A, + A, = 100; B, + B, = 100 und aus (16) fiir 
n= 1 zwei weitere lineare Gleichungen, aus denen die Koeffizienten zu 
A, = 97,90101 ... , A, = 2,09898....) By= 104,61210..., B, = —4,61210... 
ermittelt werden kénnen. Die Lésung der Differenzengleichungen ist in Figur 3 
fiir n= 0,1,..., 14 dargestellt. 


6. Die Methode der Fehlerabschatzung 


Durch den Ubergang von der Differentialgleichung und ihren Randbedin- 
gungen zu einem System von Differenzengleichungen entsteht der sogenannte 
Abbrechfehler (truncation error). Wir skizzieren hier das Prinzip der Abschat- 
zung und beniitzen dazu das oben verwendete grobe Gitter. Es ist zu betonen, 
dass fiir die tatsachliche Durchfiihrung einer Fehlerberechnung feinere Gitter 
beniitzt werden miissen. 

Wir bezeichnen den Abbrechfehler mit A(x, 7) und verstehen darunter die 
_ Differenz aus der exakten Lésung u(x, 7) der Gleichungen (1) bis (5) und der 
Lésung der Differenzengleichungen (11). Die in (6) mit ®,, R,, Rs, bezeichne- 
ten Restglieder erhalt man aus der Taylor-Entwicklung: 

2 2 2 
ee ee ee 
wo ®,, ®, Mittelwerte der 3. Ableitungen von uw nach x und 7, Y, und ¥, ent- 
sprechende Mittelwerte der Ableitungen 4. Ordnung sind. Ferner ist in (9) 
h2 


Rye i 2k. 


5 (18) 


Hierin sind //’ und f,’ gewisse Mittelwerte der 3. Ableitung von /(x). 

Die genannten Mittelwerte der héheren Ableitungen von u(x, 7) sind zwar 

unbekannt. Sie kénnen jedoch aus der Lésung des Gleichungssystems (11) ge- 

schatzt werden. Man braucht dazu nur die héheren Differenzen zu bilden und 

geeignete Formeln der numerischen Differentiation anzuwenden. Die Schatzung 
wird um so zuverlassiger, je kleiner die Gitterseite / ist?®). 

! Leitet man die Formeln (11) wie in Abschnitt 3, jedoch ohne Vernachlassi- 


gung der Restglieder her, so erhalt man die Fehlergleichungen 
Dut, + But tyut + duk=S, cwt,,+ouktqwt,+Out=T, (9) 


10) Dass der Abbrechfehler gegen Null strebt, kann auf die gleiche Weise wie in [2] gezeigt 
werden. 
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in denen u* und w; die Werte der exakten Losung von (1) und S, T lineare 


Ausdriicke in den Mittelwerten der Formeln (17) und (18) bezeichnen. Subtra- 


hieren wir die Gleichungen (11) von den entsprechenden Gleichungen (19), so 


ergibt sich ein inhomogenes System linearer Gleichungen fiir die Werte A(x, 7) 
in den Gitterpunkten. Die Fehler A(x, 7) geniigen also demselben Gleichungs- 


be he f ; AEE, : 
system wie die Naherungswerte #,, w,, nur mit anderen rechten Seiten!!). Zur 


Auflésung dieser Fehlergleichungen kénnen wir wiederum Gebrauch von der 
Theorie der linearen Differenzengleichungen machen: Wir addieren zu einer 
willkiirlichen partikularen Lésung die allgemeine Lésung des homogenen 


Systems, die wir schon in Abschnitt 4 ermittelt haben, und passen die Koeffi- 


zienten den folgenden Randbedingungen an: 


A(O,.7) = 0, kim A(x, 7) = 0. (20) 


x%—>0O 


Herrn PETER TSCHICHOLD sprechen wir fiir die Durchfiihrung der numeri- 


schen Rechnungen unseren besten Dank aus. 
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Summary 


The paper describes a method for numerically solving an elliptic type diffe- 
rential equation connected with a diffusion problem with superimposed convec- 
tion. The method consists in replacing the differential equation by a set of 
simultaneous difference equations which are solved exactly, thus avoiding 


difficulties inherent to the relaxation method. An example is treated in detail and 
the truncation error is discussed. 


(Eingegangen: 17. Marz 1955.) 
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— 


Vol. VII, 1956 121 


Ein Grenzproblem 
einer technisch wichtigen nichtlinearen Differentialgleichung 


Von ADOLF GiGER, Ziirich!) 


Einleitung 


Die nichtlineare Differentialgleichung 


q@2 
@ 40 


3 +sinO =p (1) 


dt 
spielt bei verschiedenen Problemen der Technik eine grosse Rolle. So gehorcht 
die Auslenkung (@) eines Pendels, das durch fliissige Reibung («) gedampft ist 
und auf das ein konstantes Drehmoment () wirkt, dieser Gleichung (t ist eine 
dimensionslose Zeitgrésse) [1]?). Normalerweise fiihrt das Pendel gedimpfte 
Schwingungen um seine Gleichgewichtslage @, = arcsinf aus. Ist aber das 
Drehmoment zu gross, namlich / B | > 1, so findet das Pendel keine Ruhelage 
mehr. Es dreht sich dann dauernd um seinen Aufhangepunkt. Es kann aber 
auch vorkommen, dass das Pendel bei gewissen Anfangsbedingungen trotz 
| 8| < 1 nicht in seine Gleichgewichtslage zuriickkehrt, sondern sich fortwah- 
rend dreht. Je grdésser die Dampfung «, desto weniger wahrscheinlich wird 
dieses Verfehlen der Ruhelage sein. Es gibt nun fiir jedes | | <1 einen 
Wert «,, so dass fiir « > «) das Pendel unabhangig von den Anfangsbedingun- 
gen stets in die stabile Gleichgewichtslage zurtickkehrt. Diese wird aber unter 
Umstanden erst nach etlichen Pendelumlaufen, also erst nach langerer Zeit, er- 
reicht. Es ist der Zweck dieser Arbeit, die Werte « (8) analytisch zu berechnen. 

Aus der Theorie der elektrischen Maschinen ist des weiteren bekannt, dass 
die Phasenverschiebung zwischen dem Rotor eines Synchronmotors und seinem 
Drehfeld obiger Gleichung gehorcht [2], [3]. Dann rihrt « von der Dampfer- 
wicklung des Rotors her, und £ stellt das angehangte Lastmoment dar. Das 
Riickstellmoment (synchronisierendes Moment), das durch das Drehfeld auf 
den Rotor ausgeiibt wird, ist auch hier wieder sinusférmig. Wenn die Belastung 
~ des Motors zu gross wird, kann er ausser Tritt fallen, und O(z) wird periodisch. 
Auch auf dem Gebiete der Hochfrequenztechnik ist die angegebene Diffe- 
-rentialgleichung von Bedeutung. Sie tritt beim Problem der Synchronisation 
eines Hochfrequenzgenerators durch eine fremde Frequenz auf. Dabei wird die 


1) Institut fiir Hochfrequenztechnik der ETH. Sic’ 
2) Die Ziffern in eckigen Klammern verweisen auf das Literaturverzeichnis auf Seite 129. 
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Synchronisation aber nicht durch direkte Einwirkung des Fremdsignals auf den — 
Oszillator bewirkt, sondern mit Hilfe einer Gegenkopplung (Figur 1). 

In einem Phasendetektor wird die gesteuerte Frequenz /, mit der Steuerfre- 
quenz f, verglichen. Der Detektor liefert eine Spannung proportional dem Sinus 
der Phasendifferenz zwischen den beiden Signalen. Die Spannung wird in einem 
einstufigen Verstirker, der nur ein einziges RC-Tiefpass-Glied enthalt, ver- 
stirkt und steuert darauf die Frequenz des Oszillators auf lineare Art und 
Weise. Ist die Frequenzdifferenz zwischen den beiden Signalen nicht zu gross 
(| 6 | < 1), so kann der Oszillator durch /, synchronisiert werden. Ob er wirklich 
synchronisiert wird, hangt von den Anfangsbedingungen ab. Befindet sich das 
System zum Beispiel anfanglich schon in Synchronismus und wird f, Jlangsam 


Fhasen berstarker Steverbarer 
Detek/or (7 xR) Oszilator 


Synchronisation eines Oszillators durch Phasensteuerung. 


Figur 1 


verschoben, so kann der Oszillator bis auf | |= 1 in Tritt bleiben. Dieses © 
Gebiet | 8| <1 wird mit «Halte-Bereich» (Hold-in range) bezeichnet. Bei 
gentigend grosser Bandbreite (prop. «?) des Tiefpasses, fiir « > %), wird nun 
aber unabhangig von den Anfangsbedingungen stets nach einiger Zeit der syn- 
chronisierte Zustand erreicht. Dieses Gebiet nennt man «Fang-Bereich» (Pull- 
in oder Lock-in range). 

Das Regelsystem nach Figur 1 wird in neuerer Zeit in Empfangern fir 
Farbfernsehsendungen nach dem Punktsprungverfahren (dot interlaced colour 
television) verwendet. Der kontrollierte Oszillator liefert dabei die Farbenum- 
schaltfrequenz (colour subcarrier) von 3,89 MHz, deren Phase méglichst wenig 
von derjenigen eines vom Sender gelieferten Synchronisiersignales abweichen 
soll [4]. 

Auch in guten Schwarz-Weiss-Fernsehempfangern wird zur Synchronisa- 
tion des Zeilenablenkgenerators (15625 Hz) die Anordnung mit automatischer 
Phasenregelung angewendet. Nur haben wir es hier nicht mehr mit sinusfér- _ 
migen Schwingungen zu tun, sondern mit Rechteckimpulsen und sagezahn- 
formigen Spannungen [5]. 

Bei dem mit Impulsen synchronisierten Oszillator (Impulse governed oscil- 
lator, IGO) [6] wird die sinusférmige Oszillatorschwingung im Phasendetektor 
mit einer frequenzstabilen Folge von schmalen Rechteckimpulsen verglichen. 
Der Vorteil dieses Systems besteht darin, dass der Oszillator dabei auch auf 
Oberwellen der Impulsfolgefrequenz synchronisiert werden kann. 
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Zur Mittelfrequenzstabilisierung von frequenzmodulierten Sendern dient 
eine Methode, bei der vorerst die Sendefrequenz stark untersetzt wird (um. die 
durch die Modulation verursachten Phaseninderungen gentigend klein zu 
machen). Alsdann wird diese geteilte Frequenz in einem Phasendetektor mit 
einer quarzstabilisierten Referenzfrequenz verglichen. Der Phasendetektor gibt 


dann die (nichtlineare) Steuergrésse zur Regelung des Sendeoszillators ab [he 


Da die Amplitudenmodulation bei Magnetron-Oszillatoren stets von einer 
unerwiinschten Frequenzmodulation begleitet ist, wurde ein in der Frequenz 
steuerbares Magnetron entwickelt, dessen Mittelfrequenz durch eine Regelung 
nach Figur 1 konstant gehalten werden kann [8]. 

Am Institut fiir Hochfrequenztechnik an der ETH. ist die beschriebene 
Methode verwendet worden, um ein Reflexklystron durch ein frequenzmodu- 
liertes Signal zu synchronisieren. Der Wert £ ist in diesem Fali noch zeitlich 
veranderlich. 

Die Aufgabe, die Werte a) in Funktion von # zu finden, ist bis jetzt noch 
nicht mit der gewiinschten Vollstandigkeit und Genauigkeit gelést worden. Es 
wird nun im folgenden eine analytische Methode angegeben, mit welcher «, 
beliebig genau berechnet werden kann. 


Die mathematische Lésung des Problems 


Wir fiihren vorerst an unserer Differentialgleichung (1) die Transformation 
@ = $+ 2 —arcsin# durch und erhalten die neue Gleichung 


bh + o sa — sin(# — arcsinf) = p. (2) 


Mit d9/dt = y wird alsdann (2) von erster Ordnung in y, namlich: 
y se +ay— sin(#— arcsinf) =f. (3) 


Mit Hilfe des Isoklinenverfahrens (dy/d = konstant aufzeichnen!) kénnen wir 
Lésungen der Gleichung in der sogenannten Phasenebene (y gegen #) auf gra- 
phischem Weg erhalten. Figur 2 zeigt einige Lésungen bei « < a ftir verschie- 
dene Anfangsbedingungen. £ ist dabei, wie auch im folgenden, eine positive Zahl. 

Auf der 9-Achse befinden sich in periodischen Abstanden von 2 Sattel- 
punkte S(—22, 0,2 usw.) und Zentrumspunkte Z(a + 2 arcsin #). Von den 


; Sattelpunkten strahlen nach 4 Richtungen besondere Lésungskurven aus, so- 


genannte Separatrizes. Beginnt die Lésungskurve beispielsweise im Punkt A,, 


~ der sich in dem durch zwei Separatrizes eingerahmten und schraffierten Gebiet 


befindet, so endet sie nach einiger Zeit in der stabilen Gleichgewichtslage eA, 
Sind die Anfangsbedingungen jedoch anders (Punkt A,), so geht die Lésungs- 
kurve nach und nach in die sogenannte periodische Losung zweiter Art y, tiber. 


124 ADOLF GIGER ZAMP 


1 CZ 
t—Farameter ‘ 


Figur 2 


Lésungskurven in der Phasenebene. 


Diese Kurve y, mit der Periode 2 kann durch eine Fourier-Reihe dargestellt 
werden. Setzen wir diese in (3) ein und integrieren iiber die Periodenlange, so 
liefert nur das konstante Glied der Reihe einen Anteil an das Integral (nebst A). 

Als Resultat folgt fiir den Mittelwert der periodischen Losung: : 


za] 0) dg =f. : 


4 


Wenn wir nun den Wert von « grésser machen als in Figur 2, so riickt die 
Lésung y, mehr und mehr gegen die Abszisse hinunter, sie schmiegt sich aber 
auch mehr und mehr an die aus den Sattelpunkten nach oben laufenden Sepa- 
ratrizes an. Im Grenzfall, fiir den Wert « = a, verlauft y, von einem Sattel- 
punkt in den andern. Bei « > a ist schliesslich keine periodische Lésung mehr 
moglich, und alle Lésungskurven bewegen sich unabhangig von den Anfangs- 
bedingungen in die Zentrumspunkte Z. 

Fir die Bestimmung von %) verwenden wir nun die Tatsache, dass die nach 
oben gehenden Separatrizes zusammenfallen und von einem Sattelpunkt in den 
andern laufen miissen. Diese Lésungskurve y,,, die zum Beispiel von # = 0 
nach # = 2 lauft, hat einen durch (4) gegebenen Mittelwert. Zudem lassen 
sich aus (3) die Ableitungen der Separatrizes in den Sattelpunkten berechnen. 

Wir kénnen dann y,, durch eine Potenzreihe y,, approximieren, wobei die 
Naherung stets durch die Punkte = 0 und #= 2a gehen und in diesen 


Punkten in 7 Ableitungen mit y,, iibereinstimmen soll. Schliesslich gelte ent- 
sprechend (4) das Integral 


aaa d= J, (5) 


das uns die Bestimmungsgleichung fiir «,(/) liefert. 
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Je mehr Ableitungen in der Rechnung mitgenommen werden, desto genauer 
wird der Wert von a. Eine gute Konvergenz von a» fiir alle 6 < 1 wird aber 
nur bei Beachtung gewisser Vorsichtsmassnahmen erreicht. So zeigt es sich, 
dass die héheren Ableitungen der nach rechts oben laufenden Separatrix fiir 
grossere 6-Werte divergenten Charakter haben. Das ist bei der nach links oben 
gehenden Separatrix nicht der Fall, obschon auch hier die Konvergenz nicht 
allzu gut ist. Aus diesem Grunde wird in y,, nur die erste Ableitung bei # > 0 
der Differentialgleichung angepasst. Hingegen werden bei % < 2a beliebig 
viele Ableitungen aus der Differentialgleichung entnommen und in die Naherung 
eingesetzt. 

Vorerst wollen wir nun die Ableitungen der Separatrizes in den Sattelpunk- 
ten berechnen. Die Separatrix in dem sich bei } = 0 befindenden Sattelpunkt 


wird dazu durch eine Potenzreihe 7 = Dd’ k, 0” dargestellt und in (3) einge- 
n=1 


setzt. Nachdem der Sinus ebenfalls in eine Reihe entwickelt worden ist, wird 
der Koeffizientenvergleich durchgefiihrt. Er liefert vorerst : 


ao a5 + VG) +Vi-P 
= -%-|/(2) + Vip, 


wobei die Indizes + bzw. — andeuten sollen, dass k, zu der nach rechts oben 
_ bzw. links oben laufenden Separatrix gehdrt. Fiir die weitern Koeffizienten k,, 
erhalten wir die folgenden i a eae 


a B/n! + (—1)"? (mn + 1) [kg hk, 4 + Be Pt a ae Ryo 41 Rye] 
hy = (-1)"2*2 | Peay (6a) 
(m = gerade) , 


und 


| ky, = (—1)""0° 


eo ; 2 
V(1 — B?)/n! + peur (mn + 1) [Ry kya + Fg Pc Nal 1/2 kin say (6b) 
for AP Ry +e 


(m = ungerade) . 


: k,, erhalt entsprechend der betrachteten Separatrix den Index + oder —. Und 
schliesslich berechnen sich die Ableitungen im Sattelpunkt aus den obigen 
- Koeffizienten gemiss d"7/d9”" = n! ky. 


n+2 


Fir y,, setzen wir die Reihe x a, 0’ an. An der Stelle # = 0 ist y,, von 
"selbst gleich Null. Eine erste Gleichdfiy macht den Koeffizienten a, gleich der 


rH 
: 
4 
z 
mi 
4 
FS 
‘ 
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ersten Ableitung im Nullpunkt, also a, = k;*. Eine zweite Gleichung setat Ve. 
bei 0 = 2 gleich Null. Dann folgen weitere lineare Gleichungen, die die n 
ersten Ableitungen (= 2! k,) im Punkt } < 2a festlegen. Somit haben wir 
n + 2 lineare Gleichungen fiir die m + 2 Koeffizienten a,. Fiir die n +1 Unbe- 
kannten a, bis a, 2 erhalten wir nach einigen Umformungen das lineare Glei- 


chungssystem: 
— (22) kf =a,(22)*?+a4,(2m)3+ a(2m)*+ a5 (27) 
+ (2a)k>4+ (20) kp ee az (2)? + 2 a, (22)4+ 3 4, (2m)? 


— (2m) kit —2(20) ky + (2m)? ky = ay (2)* +34; (22)? : 


Die m-te der insgesamt 1 + 1 Zeilen ist dabei allgemein von der Form: - 


n+2—m ee t =e. 2 
(-1 yeaah + 3 pNZayO Ok ee. = Lf Aman ( (2 a) Be Hees }3 
Die Lésungen dieses Gleichungssystems lassen sich durch folgenden Ausdruck 
darstellen: | 


n n+1—-x 
i ig 1 v+e ony eee T+x#x—1 ] 
(22)’ a, = (-1)"* (Po) amat + 3 (-1) (22)* k; amd at ), (7) 
x%=1 =] 
; . 
wo y=2bisw+ 2 und a, =f. 


Die + 2 Koeffizienten a, sind also durch k;* und die n Koehiienten Re 
eindeutig bestimmt. Da in den k,, nur die beiden Gréssen « und f stecken, inal 
auch die a, nur von diesen Konstanten abhangig, also a, = a »(«%; 6). Nunmehr ~ 
liefert uns (5) die Gleichung zur Bestimmung von % (6), namlich: 


: 
B n+2 a | 
ree eel a 


Setzen wir nun (7) in (8) ein, summieren und ordnen nach k,,, so erhalten wir 
schhisslich : 


B ts 2 TciRt 
Oo (n + 2) (n + 3) 


14 2 ae 
td lmaema (n+ 2) (n+ 3) _ ak 22)"k, =0. 


: 
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Es ware viel zu kompliziert, in (9) die k, durch die Gréssen « und 6 zu ersetzen. 
Wenn wir zum Beispiel bis n = 5 gehen, so lautet (9) 


B 


Xp 


- + 0,112.20 kj — 1,68299 ky + 4,6998 Ay — 13,2883 Ry 
433,397 k- — 58,28 koa 0. 


Zur Berechnung von %» bei festem Wert 6 nehmen wir vorerst eine Naherung | 
fiir %», berechnen daraus die Werte k,, nach (6) und setzen diese in (9) ein. Er- 
halten wir dabei nicht den Wert Null, so muss ein neuer Wert fiir «) gesucht 
werden. Meist geniigen zwei Versuche, um mit Hilfe der Regula falsi den 
Wert «, auf geniigend viele Stellen zu erhalten. Fiir n = 5 wurde «, fir0<$<1 
in Schritten des Wertes £ von 0,1 mit einer Handrechenmaschine berechnet. 
Durch Vergleich mit noch héheren Approximationen lasst sich erkennen, dass 
der fiir 7 = 5 erhaltene Wert «) héchstens um 1 Promille vom richtigen Grenz- 
wert abweicht. In der folgenden Tabelle sind die berechneten Werte aufgetra- 
gen, wobei die letzte Ziffer einen Fehler von ungefahr einer Einheit aufweisen 
kann. 


B | 0 | 0,1 | 0, 3 | 0,4 | 0,5 | 0,6 | 0,7 | 0,8 | 0,9 | 0.95 |0,98 | 1,0 


2 | 0, 
0,0787| 0,158 : 0,239 | 0,323 0,410 | 0,503 | 0,605 | 0,720 | 0,865 |0.964 | 1,050 | 1,193 


ay | 0 


Zur Kontrolle der Konvergenz des Verfahrens wurde fiir 6 > 0 und a > 0 
der Quotient £/x, und fiir 6 = 1 der entsprechende Wert « noch héher appro- 
ximiert. Die folgende Tabelle zeigt das Verhalten der Werte bei zunehmender 
Hohe der Naherung bis 7 = 8. Fiir £/«) kann iibrigens aus der Differential- 
gleichung exakt der Wert 4/2 gefunden werden. 


i : 
n ay |B 0 | pu 
| Xo —>O 
iG 1,05 1,38 
2 1,26 fe 
3 1,29 1,190 
4 1,276 1195 
5 1,272 1,194 
6 L2729 1,1936 
‘L 1 BEN 1,19314 
8 4273.26 alike tlts) 
exakt 1,273 24 


Die angegebenen Werte zeigen, dass das Verfahren konvergent ist. In 
Figur 3 ist der Verlauf «(@) aufgetragen. Das Gebiet zwischen der %-Kurve 
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und der «-Achse ist der Fangbereich, wo keine periodischen Losungen moglich 
sind, also stets Synchronismus eintritt. Zwischen der Kurve Xo und, f= 
liegt ein Teil des Haltebereiches, wo die Frage der Synchronisation stark von 
den Anfangsbedingungen abhangt. Fir « SP U19SL ist stets fiir B x i Syn- 
chronisation gewihrleistet. Die Abszisse von Figur 3 gilt aus physikalisch leicht 


Figur 3 
Verlauf von %p fiir die Differentialgleichung d?O0/dt? + «d@/dt + sinO = B. 


ersichtlichen Griinden auch fiir negative Werte von £, weshalb sie mit dem 
Absolutwert von £ angeschrieben ist. 

In Figur 3 sind noch die von HayEs [9] angegebenen beiden Kurven ein- 
gezeichnet, welche die (damals noch unbekannten) richtigen Werte eingrenzen 
sollen. Auch SEIFERT [10], Tricomi [11] und Amero [12] haben Grenzen an- 
gegeben, die aber sehr weit entfernt von der exakten Kurve verlaufen. PRESTON 
und TELLIER [4] fanden mit Hilfe der graphischen Isoklinenmethode fir 6 < 0,5 
die approximative Beziehung £/x) = 1,25. Auch in einer andern Arbeit [13] 
wird % nach der Isoklinenmethode bestimmt, und zudem werden Messresultate 
angegeben. Nur ist die erzielte Genauigkeit gering. 

Die vorliegende Arbeit ist ein Teil einer theoretischen und experimentellen 
Untersuchung iiber die Frequenzmodulation eines Reflexklystrons mittels 


‘ 
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Phasengegenkopplung. Sie entstand am Institut fiir Hochfrequenztechnik an 
der ETH unter der Leitung von Herrn Prof. Dr. F. TANK. 
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Summary 


The nonlinear differential equation 


2 
> ie sin@ = B 
is encountered in the treatment of several technical problems, e.g. the oscillations 
of a synchronous motor or the synchronization of a high frequency oscillator. 
The equation has two types of solutions @(t), one of them periodic with period 
22 in © and the other one with a constant final value 0). The periodic solution is 
usually not desired in practice because it means the system has fallen out of 
synchronism. Therefore it exists a need for knowing the boundary between these 
two kinds of solutions. For | | > 1 a periodic solution does always exist and 
for 0 <| | <1 there exists an «(f) such that the condition « > % implies 


the nonexistence of periodic solutions independent of initial conditions. In this 
- paper the values «(8) are calculated with the aid of an analytical approximation 


an 


method. The accuracy of the approximation can be increased to any desired 
degree. This is shown in two particular cases. 


(Eingegangen: 15. Februar 1955.) 


' ZAMP VII/9 


ZAMP 


Stress Distributions Within Solids Bounded 
by One or Two Cones 


By Antony J. A. MorGan?), Pasadena, California 


Abstract 


The Mellin transform is used to obtain a formal solution for the state of 
stress, under quite general (e.g., piece-wise continuous) loading conditions, 
within solids bounded by one or two cones. For certain ranges of the cone 
semi-vertex angles, these solids, all of linearly varying thickness and semi- 


infinite in extent, may be descriptively classified as: (I) circular plates, (II) coni- 
cal ‘shells’, and (III) solid cones. Under case III are subsumed: (a) the semi- — 


infinite solid, (b) the semi-infinite solid with conical depression, and (c) the 
infinite solid with conical exclusion. As an example the solution for case Illa, 


where a normal loading is uniformly applied over a circular area, is worked out — 


in detail and shown to check with previously known results. 


1. Introduction 


By making use of the Mellin transform formal solutions of the governing 
partial differential equation for the stress function are obtained in the case 
where the solid is bounded by one or two co-axial cones and is subjected to 


axially symmetric external loadings. If 6 > « are the semi-vertex angles of © 


the two bounding cones, then the solids (Figure 1) under consideration may 
be descriptively classified as follows: 

(I) 6 > x/2 > « circular plate of infinite radius and linearly varying thick- 

ness ; 

z/2 > 6 > « semi-infinite conical ‘shell’; 
z > B > a, one bounding cone: semi-infinite solid cone; 
p= 1/2 semi-infinite solid; 
Cem 2 semi-infinite solid with conical depression; 

(c) 6 = s¢ infinite solid with conical exclusion. 
It is emphasized that the word ‘shell’ in case II is not used in the usual sense 
since it does not carry its common connotation of ‘small’ thickness. 

The work given in this paper constitutes an extension of, and complements, 


that reported upon in [1]?), where polynomial load distributions were con- . 


1) AER, Inc. 
*) Numbers in brackets refer to References, page 145. 
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sidered’), since the present solutions cover loading cases (e.g., discontinuities) 
which were excluded from those given in [1]. However, the investigation in [1] 
is not entirely unconnected with (and, in fact, motivates) that given herein 
since in the former and latter instances, respectively, the stress function ® 
and its Mellin transform @ satisfy fourth order ordinary differential equations 
which are very closely related. Also, as mentioned, in [1], a further motivation 
for pursuing the present work was provided by TRANTER’S [2] success in using 


Toy Tig 
Toy 


r 


W (P,2)= Sal) 


a 
Zant | 


avant = mEN of. 
y (P:@) = f(p) 


Figure 1 


The semi-infinite conical shell of linearly varying thickness. 


the Mellin transform for the determination of the state of stress within semi- 
infinite two-dimensional wedges. ; 

The title problems are formulated in section 3 and their formal solution 
by use of the Mellin transform is given in section 4. The method of specializing 
the solutions for cases I and II so that they will be applicable to case III is 
discussed in section 5 and, in section 6, is applied to obtain solutions for case 
Illa. The advantage of the method becomes evident in section 6, where it 
appears to give a more direct result than could be obtained by an extension 
of the methods elucidated upon by TIMOsHENKO [3, p. 333]. As an example, 
in section 7, theory is applied for a particular loading-type to case IIIa, the 
"results are found to correspond to those given by TIMosHENKO [3, p. 335]. 


3) Precisely speaking, normal and shear loadings proportional to an arbitrary real number 


power (positive or negative) of the distance from the common apex of the two co-axial bounding 


cones. 
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2. Notation 


radius of circle over which uniform load is applied to face of semi- — 


infinite solid; 

constants of integration ; 

parameter locating path of line integral in Mellin inversion formula ; 
differential operator d/dC; 

Mellin transforms of 9? 7,(0) and 0? /,(@) respectively ; 

normal loading (force/unit area) distributions applied to conical 
surfaces at y = B and y = « respectively ; 


functions of p and ¢ defined by equations (40) to (43) respectively; — 


integrand of integral in equation (64) ; 
y—1; 


functions of p and € defined by equations (36) to Se aie 


variable ranging over positive integers: 0, 1, 2, 3, 


magnitude of normal loading over circle of acice a on the face © 


of semi-infinite solid, a constant; 
Legendre function of first kind of order w and argument ¢; 
complex parameter associated with Mellin transform ; 


Legendre function of the second kind of order w and argument C;_ 


cylindrical coordinates ; 

Mellin transforms of 03 s,(0) and 0° s,(@) respectively ; 

shear loading distribution applied to conical surface at y= f; 
shear loading distribution applied to conical surface at yp = «; 
semi-vertex angle of inner bounding conical surface; 
semi-vertex angle of outer bounding conical surface; 

Gamma function of complex variable p; 

PPA; 

cos f; 

cosy; 

COSA; 

pA; 

PoIsson’s ratio; 

spherical coordinates ; 

arbitrary order for Legendre function P,(¢); 


Jo, Sy, G9, Toy Stresses in spherical coordinates, axial symmetry imposed; 
O,, O;, 09, T,, Stresses in cylindrical coordinates, axial symmetry imposed ; 


Plo, y) 


stress function in spherical coordinates: 


Bar over variable(s): Mellin transform of indicated variable(s) ; 
Letter subscripts: partial differentiation with respect to indicated letters, un- 


less otherwise indicated by context. 
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3. Formulation of the Problem 


Both of the title problems are axially symmetric. Taking spherical coordi- 
nates (0, y, 6) the required stress distributions are determined, under the 
assumptions of classical elasticity, by a stress function ®(o, p) satisfying the 
partial differential equation [3], [4] 


2 0? 2 m) 1 n) 1 «92 \2 
V40 = V2V7?2O — et ese verss SS ee = 
(sar += car Baad ax Panes ) Cm Oe 

OSes, «<ySf), | 


and the boundary conditions 


o,=f,(o) (y=), (2) 
= flo) (p=), (3) 
and 
Toy» = Silo) (p= 8), (4) 
=s,(0) (p=2), 6) 


where /,(@), /.(@) are the normal stresses and s,(@), s,(o) the shear stresses 
applied to the faces of the conical ‘shell’ (Figure 1), all supposed given func- 
tions of o and such that their Mellin transform exists. Once ® has been deter- 
mined, the stress components in cylindrical coordinates (7, 2, 9) are given by: 


o) ‘ 0?@ 
0; = oy [» ea =r | , (6) 
t) ie yo’: 
c=. [pV - 5-1 (7) 
0?@ 
0, = — [ (2 — y) V2d — a , (8) 
and rs 
con Ba ro 30), 9 


where y is Porsson’s ratio. 
The stress components in spherical coordinates are then given by the trans- 


formation relations: 


a, = 0, sin?y + 7, Sin2 y + 0, Cos?y, re) 
Oy = Or cos*y + T,, sin2 y + 0; sin? y , (11) 
iO, (12) 
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and 7 
tos (o, — 6,) sin2 p+ T,, cos2y. (13) : 


oy 2 


4. Solution of the Problem 


Let 
C=cosy, (14) 


then equation (1) can be written as: 


dD 2 | . 8@ 


20 sap t OS ran 


Assume now that @ is such that 


n 
grin (n= 0, 1,2,3,4), @ sa lke Score 
220 : 4 
pis Dh | see ei es cg cae aA oe 
0 FOG,” and 0 00? 06? 


all vanish as 9 > oo. Then, on denoting the Mellin transform of ® by ®, that is: 


P(p, £) = | Blo, 0) e?-*do, (16) 
0 : 
integration by parts gives: 
Pant Md fa@ ab 
fe 002 Of” odo = p(p + 1) “dtm? / ocr e? 1 dg = dtm’ 
0 
fons, 02 Oe ieee = 
ee cage Ce = PP a aes (17) 
6 


and Be 


[ot eo or-1do=p (b+1) (6 +2) (+3) 8. 


On multiplying equation (15) by 9?*8 and integrating it with respect to 

from 0 to oo the use of relations (17) gives the ordinary differential equation: 
[(2? 1) D?+ 2¢D—(p—1)6] 

: on 18 

[(c*—1) D+ 2¢D—(p +1) 6+. 2)] GE; 4) =0, a 


where the operator notation D = d/df has been used. As shown in [1], the 
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solution of equation (18) is given by: 


DE; p) = Cy P,(0) + Cy Q,(0) + CeP,(E) + C1. 0,10), (19) 
where 
fe PoE tok 


P,(¢) = Legendre function of the first kind of order ~ with argument € , 


and 


Q,,(¢) = Legendre function of the second kind of order yw with argument ¢ . 


The stress components are expressible in terms of @ and its partial deriva- 
tives with respect to a and ¢. The substitution of equations (6) to (9) in equations 
(10) to (13) and the use of the transformation (14) indicates that: 


ns =(1-)0®,..+ —22—o® 


0 


ate ew ce, 


—* (1-9) (1-0) O..+ 52+ 7-324 ©, 


(20) 
s a i, 
+ = 2-9) £(L— 09) Be, — se 2[L + v— G—9) 2) B, 
= 222+) 61-04) G+ oe vb 04) Bere, 
Oy = 95 Byoq+ — (1 +29) £G,.+ oy (1-20) CO, 
+ + (2-2) (1-09) ® oe + ap (1-29) ©; 
: (21) 
2 il} 
= Se (L= EL 0) Dep Ge (1-29) L404) ®, 
— 4, B= 2950-09) Bee + Ze 1 =O ree, 
i 1 
Tp = VED og — (1-2) CPt oe (L— 20% 
ee De ol 20D, 


i < 
eA ED goer! ome en ee 


— To 
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Gea) pee Pe) Dee - (L— 2») ®,, 


ay C D, c 


+ os (1— 2) ©, +5 VE Deep — 0? 


+e? Dp ie Ae eee ar id (25)) 


4, 2[2tr—9 (+ 24]¢ 


il - es 
+, (1-1-2) Ge, 


03 66s 


+ 


, 


where the subscripts on ® denote partial differentiation with respect to the 


indicated variables. 
Denoting the Mellin transform by a superimposed ‘bar’, the followings 
relations are obtained: . | 


(0% o,) = [ 084, 0° "de, (e%0)y =| 9%, 0 * do, 
; : 
3 (24) | 


(0° os) =i 0° 0,0° ‘do, and (0? z,y) =| Q* Toy oe? * do; 
hy 


0 


furthermore, let 
= [0° fule) o? 2 de, and S,(6) = [0 sy(o) oft do (n= 1,2). (25) 
0 0 


Then, from equations (20) to (23), the following results are obtained: 


(0° o,) = » (1 — 6%)? O,.. —[(2— ) p+ 2 (2+ »)] 6 (1 — 04) @,, | 
+[6{p + (1—»)}—{(1—») p?+ (7—») p+. 2(5—»)}(1— £2] ®, | (26) 
—[(L—») (+1) (6+ 2)-2(2—») p] p29, | 


(9 o,) = (1—v) (1 — 09)? @,.. + [1 —») p- 8-2] C(1— £4) @,, 
+[(2—»)p(@+1) 1-0) -(1— 2») @ 414 2], 
+[p-—v (+1) (64+ 2] p69, 


ay a ee ee oe 
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(e%o,) =» (1-02)? G,.. + [1-» + 0]C0—6)6 | 
+[»p (6 +1) (L—£4) + (1-29) p(1— 262) + (1 — 29) (130m), (28) 
—[pp+(L—v)]p(6+2)00, | | 

and 


GH O-™ (o'r, )=— —)t 1-296, 
+[{4—» (p+ 2)}02-»(p + 2)]@ 
+[3p+2+r{p (+1) +2}]6@ 
— [+ 2)—2]p (6+ 2) 


where now the subscripts on the ® denote ordinary differentiation with re- 
spect to ¢. 

By use of equations (25) to (29), and the Legendre differential equation to 
express the second and higher derivatives of the functions P, ,(¢) and Q,, (2) 
in terms of these functions and their first derivatives, the boundary conditions 
(2) to (5) can be expressed as: 


F,(p) = [Ky (6; A) P(A) + K2(p; A) a 
+ [K,(p; A) 0),(4) + K2(p; A) Q,,(A)] Co a 
+ [K3(p; A) P(A) + Kip; A) P ne 
+ [K3(p; A) Q)(A) + Kg (6; 4) Q,(A)] Ca, 
F,(p) = [Ki (6; A) Pi (A) + Ke(p; A) P me 
+ [K,(p; A) Q),(A) + Ka(p; A) Q,(A)] Ce = 
+ [K,(p; A) P)(A) + Kalb; A) & an 
+ [K3(6; A) Q;(A) + Ka(p; A) Q,(A)] Ca, - 
(1 — A2)! S,(p) =[H,(p; A) Pi(A) + Ha(p; A) P,(A)] Cr 
peas se A) Q/,(4) + H2(p; A) Q,(A)] Ce Bs 
+ [H,(p; A) Pi(A) + Half; 4) P(A)] Cs 


+ [H;(p; 4) Q;(A) + Ha(b; 4) Q,(A)] Ca, 
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ns (1 — A)? S,(p) = [Hy(6; A) Pi(A) + Help; A) P(A) C, 
+ [Hs (6:4) (A) + Halb; A) QA Cs |. 
+ [H,(p; A) P’(A) + H,(p; A) BA) Cz 
+ [H4(p; A) Q,(A) + Hilp; A) Q,(A)] Ca, 
where 
A = cosa , (34) 
A = cosp (35) 
K,(6; 0) = (@+ 1)? -[p @ +1) +1], (36) 
Kp; )=-(*-p-I pe, (7) 
Kp; ) =—[p — 2 (1—»)] [6 — 2 — 2 »)] 6*@- 0-2») B — 29) 2 | Bee 
he = Geari20) PP nade) eke | | 
K,(p; 0) =—[2? @+ 1) + 1-2) (6 +2) (2p+ 1] 6, (39) | 
Hyp; 0)=—4 (+204 p@4+HC0-24, (40) 
Hp; t) =2 (6) (6-1) (2) +P +1 42-29 | a 
—p(p+220- 2), 
H,(p; C) = 10 08-2 (2p+1) C4 (627 4+6p4+2 C(1— 24), (42) 
Ha(p; 0) =2¥ (6-2) (@+1) (6 +2)—¥ (+2) (P*-7H—6) | a 
x (1-0) —p (6 +2270 - 0), 
and, for instance: P(A) H! [ae] ? : me 
The Mellin transform ® of the stress function ® is thus, in principle, 
completely known since the arbitrary constants Cy, ..., C, can be determined 


by use of equations (30) to (33). Once this is done the stress function @ is 
determined by use of MELLIN’s inversion formula RAE 
€+100 


Be, = 55; | Ob: De Pap. (45) 


€—1t00 


: 
: 
: 
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Since more interest, however, is attached to the stress components themselves, 
a further application of MELLIN’s inversion formula yields the following ex- 
pressions: 


€+tco Pa aoe 
1 f= 1 : 
So = 3 =P L of 
@ % i‘; 2:56 4 fi (@ 0.) g dp ‘ Q° Oy cae 2 71 | (o° 0) Q ‘ dp , 
c¢-—t0o c een 
(46) 
e¢+tco €+100 
he 2... f ee 1 aed 
a ae 3 -p Ce 2 
Q % J 2% 1 j (@ o) Q dp . and Q° Toy = 2 Ti / (e° Toy) Q a dp , 
c= 180 c—100 


where it is understood the barred quantities in equations (46) are to be replaced 
by their corresponding expressions from equations (26) to (29). This completes 
the formal solution of the problem. 


5. The Stress Distributions Within Semi-Infinite Cones 


On a cone only the two quantities o, = /,;(e) and t,, = s,(@) can be pre- 
scribed on its surface. The general solution (19) has four constants of integration, 
two of which must then be disposed of. To do this note [6] that the functions 
P,,,(¢) and their derivatives are regular while, on the other hand, the functions 
Q,,,,(¢) and their derivatives are singular at € = 1 (y = 0). Thus, the physically 
reasonable ‘over-all’ restriction that the stresses must remain finite on the 
cone’s axis (€ = 1) can only be realized by setting C,= C,=0 in equation (19)*). 
The applicable boundary conditions are then given by equations (30) and (32). 


6. Stress Distributions Within Semi-Infinite Solids 
Subjected to Mellin Transformable Axially Symmetric Loadings 


In this case the cone’s semi-vertex angle f is taken to be /2. As indicated 
in section 5 the Mellin transform of the stress function must be given by: 


D(C; p) = C, P,(C) + Cz P,P (6) . (47) 


It is simpler to work with the stress components o, and T,,,, when expressed 
in terms of the spherical coordinates (0, ¢ = cosy). The expressions for o,(@, ¢) 
and t,,(0, ¢) are obtained by applying the following transformation operators 


0 m) i 
OF Be 6 


0) m) il 
(1 — E} ae and yji— ti . mae FY: 0 ¢ Ya 


to equations (8) and (9). On taking the Mellin transform of 0% 0,(0, €) and 


4) However, some initial stress situations would be amenable to analysis by adopting the 
opposite point of view and setting C; = Cz = 9. 
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0? T,-(@, ¢) m the manner indicated by equations (24) and using the relations 
(17), it is possible to write 


0° a,(¢; p) | 
=(2—»[—C-1) p6+2)9 | 
+{(1— 0) (@-2) +1) +207 @+2)}@, 

| —£(1— 0%) @ +6) Bp + (L—L9)? Bee Keene 
+[e2@+H—-30 CV] le G4+20 | 
—3 (1 2) [08 +2) 1-09 + 1)| ® 
+ 30(1— 292 (p +3) D,— (1-29 G, 
and 
a: Tre(C; B) 
=(1—») [— (6-1) p+ 2) ©—-26 ? —p — 3) @, 

—{(1— 0%) (p+ Bienes = Bikitigs ae (49) 
Se ne ee (Pe oy Pape ae 
— —C[0— 22) {2p (6 +5) + 93-02 (H+ 2) (6 +3) ©, 
+(1-2)[G— 29 @4+2)— 294710, eco. | 


Now, equations (47) to (49) when evaluated at €=cosy=0 (w=a/2), 
the face of the semi-infinite solid, will yield the two required boundary con- 
ditions: 


0.0; 6) =[292h+ N+P 1G" ++I Sel 
a (50) — 
ae Pt) 8 GF lpg ~P) . 


and 
Q° tr 2(0; p) = [2 » (26 + 1) +f?) Cy BO) + 2 C, P0) = S,6), (61) 
where 44(p) and S,() are as defined by equation (25). The values of a Legendre 


function P,(¢) and its first derivative at ¢ = 0 can be expressed as [5, p. 63] 
or [6, p. 144]: 


fo) = FEE oat P= af’ 2) 
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and 
Pa) 


at 


RQ lata ahd (Lad, 
lt-0 «=©9F[(1 + o)/2) A{—o/2] ’ (53) 


where o can be either a real or complex quantity. On substituting equa- 
tions (52) and (53) in equations (50) and (51) (remembering that in one case 
o = «= p—1 and in the other o = y = p+1) and using the property of the 
Gamma function I’(1+ 0) = oT oe), the two boundary conditions can be 
written as: 


[2 (2p+ 1) + p*]€,— 9 (6 +1) C,— (2~4+- 1) Cy | 
ee: rh teepy are Ba 
Fea ee 
and 
[2%(2p+1)+ p7]C,—p(Ot+)Q 
(55) 


2a 


The constants C, and C;, parametrically dependent on #, can now be deter- 
mined from equations (54) and (55); they are given by the expressions: 


1 1 
2p+1 2a 


Eos 


[pbx ES) r(B) noon 4%) 71-8) 50], 


and, in terms of C3, 


C= “pay [29 2h +1) + PCs = 41-34) r( AOE 


b(p+ 
(57) 


Thus, the constants C, and C, have been determined for the case where (arbi- 
trary) Mellin transformable axially symmetric normal and shear loadings are 
applied to the face of the semi-infinite solid. : 

Focusing attention on the stress o, at the symmetry axis” (p ==/0 or 
¢=cosy=1) of the external loadings and within the semi-infinite solid, 
equations (47) and (48), with the aid of the relations 

By =. 
and 
dPo(6) 1 


gael o(o+-1), 
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yield the following expression for 0? 0,(1; p) 
@ o,(1; p) = + 2) [{2> 2p +1) + P7}Cs—PO+NG] | 


(58) 
+ (P42) O44) 2D+N Cy. 
By equations (55) and (56), equation (58) can finally be written as: 
1a) = cL [Ct O+4 4 (t=?) r(P) Rp 


— +1) +2) +3) PZ") (4) S)]. 


The stress o,(0, € = cosy) on the axis (y = 0) is then given by applying the 
Mellin inversion formula to equation (59), that is: 


c+i00 


ole, = zaz | ealtiP)e? ae, (60) 


€—1t0 


where 0? o,(1; £) is given by equation (59). 
Thus, equations (59) and (60) give a general relation connecting oa,(g, 1) 


with the arbitrary external loadings, which are characterized by A(p) and — 


S,(p). Of course, as previously mentioned, the external loadings are arbitrary 
to the extent that their Mellin transform exists. In the following section 
equations (59) and (60) are used to determine o,(0, 1) in a particular case. 


7. Example: Semi-Infinite Solid 
Subjected to a Uniform Normal Loading Applied Over a Circular Area 


In this case, equations (2) and (4) specialize to: 


POs Oaag 
fi(e) = (61) 
OF as ose 
and 


, . (62) 


where a is the radius of the circle over which the uniform load per unit area 


- is applied. The quantities F,(p) and S,(s) are determined by equation (25) 
that is: 


, 


P qPr3 


F,(p) Hee end anSa ie eae (63) 


ooo 
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Hence, on using the results (63) in equations (59) and (60), 


(2) ale, 1) = - 7 ii i “er ers? | (64) 
SA raya. | 


The line integral in equation (64), taken in the complex p-plane, can be 
evaluated by application of CAucHY’s theorem of residues. The only singu- 


larities of the Gamma function J"(z) are simple poles at the points z= —n 
(a = 0, 1, 2, ...). Since [7, p. 207] 
‘P\ _wi(P\ 15 yy _ (-1)" 
I5) = P(5) +2 nl (p+ 2n)’ e 


where Y’(p/2) is an integral function, it is seen that the singularities of I"(p/2) 
at p= 0, —2, and —4 in the integrand .7 of equation (64) are only apparent 
and, therefore, removable. Thus, \7has simple poles at p= —1, —3, —(2n+6), 
and (2”+1),n=0,1, 2,.... For —3<c¢< —1 and o/a <1, the line integral 
of 7 taken over a path composed of a circular arc of radius R to the left of 
(and including) the vertical line with end points at c—7R andc+iR will 
vanish on the circular arc as R + oo. Therefore, the sum of the residues of ? 
at the points f = —3 and — (2 +6), gives: 


é : =f oe eas Vx eh [1 & ed pier yan | ; 2 <1." *{(66) 


The same expression will be obtained for e/a > 1 by evaluating the line 
integral of .7 over a path composed of a circular arc to the right of (and 
including) the vertical line with end points c—1R and c+7RK and letting 
R >oo. On substituting equation (66) in equation (64) one finally obtains: 


o,(0, 1) = —P [-1 r (@ rar gual (67) 


which, except for a change in the sign of P due to the difference in the direction 
of the applied load, is exactly the same expression obtained by TIMOSHENKO 
(3, p. 335] on using an entirely different method. 


8. Discussion and Conjectures 


It has been shown that by use of the Mellin transform, it is possible to 
determine the stress distributions within conical ‘shells’ of linearly varying 
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thickness, cones, and semi-infinite solids under quite general conditions of. 


surface loading. Once the stress distributions within the aforementioned solids 
are known, it is of course possible to determine the corresponding displace- 
ments by use of the stress-strain relationships [4]. The generality of the method 
of solution is evidenced by the number of cases (see section 1) which can be 
subsumed thereunder. 

Attention should be directed to the fact that all of the solids considered are 
semi-infinite in extent so that boundary conditions need not be specified at 
the surfaces produced by cutting these solids by a plane perpendicular to their 
axes of symmetry or by a spherical surface centered at the origin. Such surfaces 
will, for convenience, be called ‘cross-wise’. The solids obtained as indicated 
above would be the ‘finite’ counterparts of those considered in this paper. To 
determine the stress distributions within such finite cones or ‘shells’ additional 
axially-symmetric boundary conditions, in static equilibrium with those speci- 
fied on the bounding cones, would have to be given at the ‘cross-wise’ surfaces. 
It may be possible to obtain solutions for such ‘finite’ problems by using 
superposition methods. That is, the stress function would be decomposed into 
two terms: ®, (corresponding to the solutions given in this paper) satisfying 
the boundary conditions imposed on the conical surfaces, and ®, a correction 
term which satisfies homogeneous boundary conditions at the conical surfaces 
and the boundary conditions at the cross-wise surface of the finite conical solid. 
The stress function ®, may in some cases be obtained by superposing an infinite 
number of ‘eigen’ functions determined by use of the solutions (16) of [1] or, 
possibly, by a variational method similar to that used by SILVERMAN [8] in 
working with the analogous two-dimensional finite wedge problems. However, 
an objective evaluation of the above-suggested methods must await further 
detailed analytical investigations. 

A problem of technical importance to which the theory developed in this 
paper might be applied is that of the determination of the state of stress 
within conical pivots. Related to this problem are the stress concentrations 
produced by applied loads over the faces of conical depressions in semi-infinite 
solids (case IIIb, section 1). In this instance, a limiting situation is that of 
the stresses concentrations produced in an infinite solid by loadings over the 
surface of a conical exclusion (case IIIc, section 1). Such investigations might 
find applications in determinations of the stress distribution within the material 


surrounding conical pivot-seats or certain types of lubrication passages within 
machine elements. 


9. Acknowledgments 


The author is indebted to Dr. Y. T. Wu, at the California Institute of 


Technology, for having generously given of his time in the discussion of certain 
aspects of the problems considered in this paper. ‘ 


OO  ———— 


Vol. VIT, 1956 Stress Distributions Within Solids Bounded by One or Two Cones 145 


REFERENCES 


[1] A. J. A. MorGan, Stress Distributions in Semi-Infinite Solids of Revolution, 
ZAMP 5, No. 4, 330-341 (1954). 

[2] C. J. TRanTER, The Use of the Mellin Tvansform in Finding the Stress Distvi- 
bution in an Infinite Wedge, Quart. J. Mech. appl. Math. 7, Part 2, 125-130 
(1948). 

[3] S. TIMOSHENKO, Theory of Elasticity (McGraw-Hill Book Co., Inc. 
York 1934). 

[4] A. E. H. Love, A Treatise on the Mathematical Theory of Elasticity, Fourth 
revised edition (Dover Publications, New York 1944), pp. 274-276. 

[5] W. Maanus and F. OBERHETTINGHER, Formulas and Theorems for the Func- 
tions of Mathematical Physics (Chelsea Publishing Co., New York 1954), 

[6] A. ErpEty1, W. Macnus, F. OBERHETTINGHER, and F. G. Tricom1, Higher 
Transcendental Functions, Bateman Manuscript Project, California Institute 
of Technology 7 (McGraw-Hill Book Co., Inc., New York 1953). 

[7] E. T. Copson, An Introduction to the Theory of Functions of a Complex 
Variable (Oxtord University Press, London 1946). 

[8] I. K. SitverMaAn, Approximate Stress Functions for Triangular Wedges, J. 
appl. Mech. 22, No. 1, 123-128 (1955). 


, New 


Résumé 


L’intégrale Mellin de transformation est utilisée pour obtenir une solution du 
probléme de distribution des tensions dans les solides délimités par un ou deux 
cones et dans des conditions tout 4 fait générales comme, par exemple, en cas 
de tension continue mais sporadique. Pour certaines limites des demi-sommets 
des angles des cénes, ces solides, tous d’épaisseur a variation linéaire et semi- 
infinis, peuvent étre classés descriptivement comme suit: 

I° plaques circulaires, II° enveloppes coniques et III° cénes solides. Dans le 
troisiéme cas, il y a lieu de subdiviser en a) solides semi-infinis, b) solides semi- 
infinis 4 dépression conique et c) solides semi-infinis, a exclusion conique. A titre 
d’exemple, dans le troisi¢me cas, oi une tension normale est appliquée uni- 
formément a une surface circulaire, la solution est détaillée ici et il est prouvé 
qu’elle est vérifiée par d’autres résultats déja connus. 


(Received: May 9, 1955.) 
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Sur Pintégration approchée d’équations du type parabolique 


Par CHarRLES BLaNc, Lausanne?) 


1. Introduction 


Le but de ce travail est l’étude au point de vue stochastique de l’erreur 
commise dans l’intégration approchée d’équations du type parabolique; la 
méthode mise en ceuvre est celle quia été utilisée pour des problémes semblables?).. 

Pour ne pas alourdir inutilement les calculs, on supposera que Véquation a la 


forme . 
OX OX 

= se Ny 1 
19. G Ot q Ox? (%, t) ( ) 
et que l’on se propose de déterminer d’une fagon approchée une intégrale satis- 
faisant a une condition initiale ; 
(7, O) eG (aye (1. 2)] 


3 Dt Sure 1 
les méthodes approchées seront des méthodes de différences finies. On se limitera 
du reste a l’étude de l’erreur sur wn pas a partir de valeurs supposées exactes. | 

: 
a 


2. Méthode du petit opérateur 


Désignons par X(x, t) une solution exacte du probleme constitué par (1. 1) 
et (1. 2)3); nous cherchons une valeur (approchée) X°® de cette solution, pour 
t= k (> 0), cette valeur étant calculée a partir de valeurs de G(x) pour ¢ = 0; 
pour cela, nous commencons par le procédé le plus simple, qui consiste 4 rem- 
placer dans l’opérateur L de (1. 1) les dérivées par des quotients de différences; 
nous choisissons : 
+ [X(u, #) — Xx, 8 — A] 


ILD BS 


— IX (+ hy t— hk) — 2X [x,t B) + X(x— ht B] 


1) Institut de mathématiques appliquées de l’Ecole polytechnique de l’ Université de Lausanne. 
Etude subventionnée en partie par le Fonds national de la recherche scientifique. M. W. Lrn1GER, 
math. dipl., a contribué dans une importante mesure a l’élaboration de ce travail et aux calculs 
numériques. | 

*) Voir Cu. Banc, Etude stochastique de erreur dans un calcul numérique approché, Comment. 
math. helvet 26, 225—241 (1952); Sur les formules d’intégration approchée d’ équations différentielles, 
Arch. Math. 5, 301-308 (1954). | 


Voir également: Cu. Branc und W. Linicer, Stochastische Fehlerauswertung bei numerischen 
Methoden, Z. angew. Meth. Mech. 35, 121-130 (1955). | 

3) Il n’importe pas, dans ce qui suit, que les données initiales soient relatives A toute valeur 
de x, ou seulement a un intervalle (a, b), avec conditions pour x = aet x = b, tf > 0. On suppose 
simplement que les données sont telles que le probléme a au moins une solution, et que les valeurs 
considérées de (x, t) sont dans le domaine de définition de cette solution. 
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1 
ox y [LX(4, 1) — pb X(x = h,t— hk) + (2p— 1) X(#,t—k) —pX(%+h,t—k)] 


avec ase 
ee 
a 2 F5 

la solution approchée satisfait A la relation 

) Lo XX F (2. 2) 


avec 
(4; 0) = G(x), 
as pei méthode, car nous négligerons ici celles qui résultent de 
| ee ee (253) 
avec 
PY=(L—L)X pour t=k, Y=0 pour ¢=0. (2.4) 


" Considérons maintenant X(x,¢) comme une fonction aléatoire, stationnaire 
d’ordre deux; supposons de plus que 


ee. 2) 00) 
et que le spectre s(x, 8B) de X(x, ¢) est de la forme 
S(a, B) = s,(x) s9(B) (2. 5) 


(cette hypothése a pour unique but de simplifier les calculs numériques qui 
Suivront; elle n’est pas nécessaire pour une étude générale); en d’autres termes, 
ona 


co 
EX (sy by) Xe t) = ff sy(a) 54(B) fl *O-9 aa ap 
—oo 


(2.6) 
~ °° 
=| s(x) ef *)2 dae | sa(B) et) ap; 
—~ —oo 


comme X est aléatoire, F, G, X® et Y le sont aussi; puisque X° est une combi- 
maison linéaire de valeurs de X, on a EX®= 0, d’ou, par (2.3), EY =0. 
Soient maintenant deux points P(%,, t,) et O(%2, t,); posons 
u(P,Q)= E Y(P) Y(Q); 


Vétude stochastique de Verreur consiste en premier lieu (et nous nous bornerons a 
cela) a déterminer la variance de l’erreur, donc cette fonction u pour P = Q; or ona 


4 L4, Li, u(P, Q) = (L°— L)p (L° — L)g E X(P) X(Q)*) (2.7) 


avec 
u(P, Q) = 0 (2.8) 


si P ou Q sont sur? = 0. 


4) Nous supposons pour cela que les opérateurs L et E sont permutables, ce qui implique cer- 
taines conditions de régularité pour s, et sy, conditions remplies par les fonctions qui seront choisies 
plus loin. 


M 
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Posons encore 


w(P, Q) = (L° — L)p (L° — Lig E X(P) X(Q), (2.9 
Vor 
dou 1s L), uP 0) = w(F, ©)5 (2.10 


la solution z = k w pour ¢ = k; il en résulte immédiatement que Yona 


UP AP) = BOP 2) (2.19 
sit= k au pomt P. 
a dp ; 
vere 


remarquons que si une fonction 2(P) est nulle pour ¢ = 0, l’équation L° z= w ; 


Calculons alors w(P, Q); par (2.9), ona 
0° 
w(P, Q) = ffs 3,(B) (L8 = Typ OOO (Te 
00 
posons 


1 MR p th + be 4.1 28) — Ph ahd BB oe Gite) alee 


alors 


(Lo — L)peilmathA) = ei2th) oly, By, . 
(L? — LD) et tanrael = - ei 2+hP) Gy, B) , 
GONG Sia 4 — i) a eo, 
2 | | 
i ” ay 
i elf sa(a) 59(B) ef" *2*"*) | (a, B) IP dae dB , : 
Jd, : 
Gy SO 
co 
1 
w(P, P) = se ff se) S2(B) | p(a, B) |? da dB, (2.13) 
as J 
donc, par (2.11) et (2.13), i 
co 
B v¥x, b) = ff (2) 54(8) | ole, B)|8 da dB (2.14 
= Co: 


Choisissons maintenant le spectre particulier) 


_1 1 
n= |For PT IS. iy [ape Pow Melson 
0 pour ja) > "mss | 0 pour |f| >a. 


5 ; eres 

) Ce choix Concust évidemment aux calculs les plus simples possibles; il a en outre l’avantage 
de correspondre a l hypothése que fait en général le calculateur, A savoir que des variations trés 
rapides de la solution sont trés improbables. On pourrait poursuivre les calculs avec un autre 


choix pour le spectre, mais l'information que l’on retire dans ce cas est déja suffisante pour une 
comparaison des diverses méthodes. 
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| On a alors simplement 


EY*(x, ky = — oak -/ \o(e B) |? dB (2.15) 


Wy, Ws 


} (en tenant compte du fait que la fonction p est paire dans les deux variables). 
Nous reviendrons sur la mise en ceuvre de cette formule. 


’ 3. Emploi d’un grand opérateur 


On peut évidemment substituer aux dérivées des opérateurs aux différences 
plus fins que ceux qui ont été considérés au § 2; cela conduit a tenir compte d’un 
“plus grand nombre de valeurs de la fonction (dou Vexpression de «grand» 
opérateur). Il est utile de rechercher dans quelle mesure la solution ainsi obtenue 
est meilleure que celle que donne le «petit» opérateur. 
Choisissons par exemple l’opérateur défini par 


Rei — + [3 Xie. ) = 2 1 ee ee — 8 | 


} 


BoE 


? 


(3.1) 


=} Ale + hi, t— kit 7 X(t 24)] | 


qui suppose donc connue la solution exacte sur les deux lignes t — k et t— 2h. 
Les considérations précédentes s’appliquent presque sans modification. On a ici, 


a la place de (2.11), 


Pe Pls eo WLP, PY, (3. 2) 


puis, en posant 


é s , ee he 
lo, B) = > — oF (p eit + pete +2 2p) 4+ 56 20hB pb Atar®—ik B, 


(3.3) 
il vient 00 
w(P, P) = gx [f sala) sy(8) | ole, B) |? de a8 (3.4) 


puis). Y 2 (x,:k)-. 


4. Méthode plurilocale 


Considérons maintenant un procédé plurilocal®); remplagons l’équation (1.1) 
Bear 


avec 


SOX° = TF (4.1) 


ASPX = 2 X(x,1) — (3p + (Xe — hy t— B+ Xv st AY] | 
£2 (3 p— 1+ o) X(x,t— 2) — (F + 2e) X(xF- 28) | (4.2) 


4 e(X(*¥ —h,t— 2k) + X(w# +h, t— 28) 


8) Mehrstellenverfahren; voir L. COLLATzZ, Numerische Behandlung von Differentialglerchungen 
(Springer-Verlag, Berlin 1951), p. 235 et suiv. 


4 
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f 

et 
¥ 


TF = F(z, t) — (e— p) (F(x — h, t— b) + F(x +h, t— By) (4.3) 
—2(p—1-— 2) Fi, t= F)s 


« étant un paramétre qui peut étre choisi arbitrairement; la solution approchée 
¥X° introduit une erreur Y = X° — X;; ici, pour la ligne calculée, | 
SV=(TL—S)X, (4.4) 


Y étant nulle sur les lignes ot X est connu exactement (lignes de départ). ; 
Ici, nous poserons 


w(P, Q) = (T?L — Sp (T9L — S%)p E X(P) X(Q), (4.5). 
t 
; u(P, Q) = EY(P) Y(Q); 
on a alors 
59, 58, u(P, Q) = w(P, Q); | (4.6) 


on a simplement ns | 


u(P, P) = == w(P, P); (4.7) 
on a ici encore © : 
w(P, P) = gx ff se) selB) lola, 8) 1? da ap (4.8), 
mais avec (si ¢ = 0) ie 
pa, f) => — phta—ikB | 
+ 2[3p(1— cosh a) —1+4+ (ph? a2+ik B) (b— pcosha— 2) e~'*F (4.9) 
= > g-2kB | 


5. Résultats numériques 


La mise en ceuvre des formules développées ici permet de répondre a quelques 
questions d’intérét pratique. Considérons par exemple la question suivante: 
une équation étant donnée, avec des conditions initiales fixant la solution, quelle 
est l’influence du choix de la méthode, ainsi que deh et k, sur la précision du résultat ? 
Nous pouvons pour cela fixer w, et w,, ce qui ne restreint pas la généralité 
(puisque cela ne fait qu’impliquer un choix d’unités pour » et t). Faisons également 
(ce qui est restrictif) g = 1; on verrait que l’allure générale des résultats reste la 
méme pour d’autres valeurs de ce paramétre. Cela étant, on calcule u(P, P) pour 
chacune des trois méthodes, pour h = 0,0 (0,2) 1,0, ki= 0,2 (0)2)°1,0. Remar 
quons en passant que la valeur h = 0 correspond en fait A un cas limite, obtenu 


en faisant tendre vers zéro le pas selon O x. On obtient ainsi les valeurs reportées 
sur les tableaux ci-dessous. 
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Valeurs de Vu- 104 


h | Rk = 0,2 | 0,4 0,6 | 0,8 1.0 
a eee cies OO ne — Satan Sa a 
I. Petit opérateur 

el eee Se ae ee aS 
0,0 135 534 1184 2071 3176 
0,2 134 a32 1181 2067 3168 
0,4 132 527 Ahilyas: 2055 3155 
0,6 130 520 1161 2036 3130 
0,8 128 512 1145 2011 3095 
1,0 128 503 1126 1942 3044 
Il. Grand opérateur 
0,0 64 241 ou. 861 1291 
0,2 64 241 510 860 1286 
0,4 64 240 508 855 1282 
0,6 65 239 505 850 1273 
0,8 67 240 503 843 1261 
1,0 72 243 502 838 1248 
III. Méthode plurilocale 
0,0 56 231 538 986 1584 
0,2 53 224 526 969 1561 
0,4 44 208 501 935 1516 
0,6 = 178 458 861 1433 
0,8 KZ 138 393 783 1318 
ihe, 34 88 amis) 639 ENT: 


Ces tableaux montrent tout d’abord que la précision des résultats est plus 
grande avec les méthodes II et III qu’avec I; ce fait pouvait étre prévu qualitati- 
vement; par contre, on remarque aussi, et c’est plus inattendu, que l’erreur 
(au sens ow nous la considérons ici) ne diminue pas nécessairement lorsque h 
diminue; au contraire, si h > 0, et cela surtout dans la méthode plurilocale, cette 
erreur augmente sensiblement; on voit ainsi que si l’on veut accroitre la précision 

des calculs (pour une méthode donnée), il faut, ayant choisi k, chosir judicieuse- 
ment la valeur de h. D’autre part, la précision augmente évidemment si & diminue, 
4 vrai dire pas dans une mesure aussi grande que semblent 1’indiquer les résultats 
ci-dessus; en effet, en diminuant # de moitié (par exemple), on double le nombre 
de pas nécessaires pour atteindre une valeur donnée de #; l’étude de V’erreur 
pour une valeur donnée de ¢ (> 0), en fonction de h et k, exige l’étude de l’erreur 
commise sur plusieurs pas de calcul; cette question sera abordée dans un autre 


travail. 
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Zusammenfassung i 


Es wird in dieser Arbeit die stochastische Methode der Fehlerauswertung auff 
Naherungslésungen durch Differenzenmethoden fiir die Warmeleitungsgleichung 
behandelt. Die Methode erlaubt die Genauigkeit der verschiedenen Differenzen- 
methoden zu vergleichen sowie den Einfluss der Schrittlangen in den beiden” 


Variablen abzuschatzen. é 
4 


(Regu: le 27 septembre 1955.) 


Uber den Druck in Ultrazentrifugen 


Von HEINRICH GREINACHER, Bern?) 


ee 


Die Druckverh4ltnisse lassen sich am einfachsten an einem mit Fliissigkeit ge- 
fiillten Hohlzylinder iibersehen. Wird ein solcher in Rotation. versetzt, so ent- 
stehen durch Wandreibung und innere Reibung der Fliissigkeit die in Figur 1- 
angedeuteten Strémungen, die iibrigens infolge der Coriolis-Krafte nicht, wie 
gezeichnet, in radialer Ebene verlaufen, sondern aus ihr heraustreten. Sie haben — 
bekanntlich ihre Anwendung in der hydraulischen Kupplung gefunden. In dem 
Masse, wie die Fliissigkeit allmahlich in Mitrotation versetzt wird, nehmen diese . 
Str6mungen ab und miissten im Idealfall verschwinden. Das heisst die Fliissigkeit 


wiirde in allen Teilen mit der Winkelgeschwindigkeit wm des Hohlzylinders rotieren. 


Stro6mung und Kompression in rotierender Berechnung des Druckes in rotierenden 
Flissigkeit. Flissigkeiten und Gasen. 


In diesem Falle lasst sich der durch die Zentrifugalkrafte bewirkte Fliissigkeits- 
druck angeben. j 


Bezeichnen wir den Druck im Abstand x von der Achse mit fy (Figur 2), so 
haben wir fiir seine Zunahme auf der Strecke dx % 


dp, = dm, Cra 


1) Physikalisches Institut der Universitat. 


Vol. VII, 1956 Kurze Mitteilungen — Brief Reports — Communications bréves 153 


wo dm, das Massenelement der Flacheneinhéit bedeutet. Da dm, = odx 
(e = Dichte), so ist 


dps o100* Kae 1 
und daher pik “) 
ow* a 
Px =e 2 = oe (2) 
Fiir den Wanddruck erhalt man entsprechend 
@ 
p= Se (3) 


Dieser Wanddruck entspricht iibrigens durchaus nicht etwa der gesamten, der 
Fliissigkeit innewohnenden Zentrifugalkraft, bezogen auf die Flacheneinheit der 
Wandung. Hierfiir ergibt sich der kleinere Wert 9 w? 72/3. 

Die Kenntnis der Druckverteilung kann nun offenbar dazu beniitzt werden, 
um die Kompression der Flissigkeit zu berechnen. Insoweit eine lineare Abhan- 
gigkeit der Dichte vom Druck angenommen werden darf, das heisst unter An- 


setzung von 
@ = Q (1 + & p,) (4) 


(@) = Dichte bei Atmospharendruck, k = Kompressibilitatskoeffizient), haben 
wir zunachst 


2 
tgs tea: 


Die im Zylinder enthaltene Masse m ergibt sich aus 


dm=22xhdx 0 


(dm = Masse in einem Hohlzylinder der Dicke dy und der Lange A) zu 
iG 
m=2a%h i oxdx, 


0 


das heisst 
4 


m=2ahoo| xdx(1+h g 
0 


Indem man im Klammerausdruck unter Vernachlassigung hdherer Glieder 
0 = Q, setzt, ergibt sich 
k 
m=rhov?+ ah a 4 Oy a 
Da zr2h gy die Fliissigkeitsmasse m, vor der Rotation bedeutet, so ist 
k 2 2 
mM = My + My 4 Qo WT". 
Das heisst die Massenzunahme infolge Kompression ist 


k 
M — Mg = ~My Qo iT? . (6) 
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—_— ee 


Da sich m — m, experimentell bestimmen lidsst, bietet sich auf Grund von (6) — 
eine einfache Methode zur Bestimmung des Kompvressibilitatshoeffizienten, Indem — 
man auf die Achse des Rotors eine Kapillare aufsetzt (Figur 1), lasst sich die ~ 
Volumverminderung AV aus der Senkung des Meniskus von 4 auf B entnehmen. — 


Da m — m) = 0 AV, so folgt aus (6) 


AV = Em ote? (7) 


Beispielsweise ergibt sich bei Wasserfiillung mit & = 50 - 10~1* dyn~* cm? und 


den Werten A = 2cm, r= 5cm und w= 22-100s-!: 4dV = 0,019 cm’. Dies J 


entspricht bei einem Kapillardurchmesser von 1 mm einer Senkung 4 > B von 


24mm. Der Wanddruck betragt in diesem Falle 5 kp/cm~?, immer unter der — 


Voraussetzung, dass sich in der Fliissigkeit vdllig statische Verhaltnisse ausge- 
bildet haben und dass nicht durch restliche Stromungen ein teilweiser Druckaus- 


gleich erfolgt. Ob dies zutrifft, lasst sich im Prinzip mit Hilfe von (6) priifen. — 
Man verwendet eine Fliissigkeit mit bekanntem # und priift experimentell, ob — 


die gemessene Volumdifferenz AV den Wert der aus (6) berechneten erreicht. In 


Wirklichkeit wird man eine etwaige Ausdehnung des Rotors mitberiicksichtigen — 


miissen, was durch Ausfiihrung einer Differenzmessung mit zwei Fliissigkeiten 
geschehen k6énnte. 


Wesentlich einfacher gestaltet sich die Priifung der Druckverhaltnisse bei den ~ 
Gas-Ultrazentrifugen. Hier werden auch tatsachlich hohlzylindrische Rotoren — 


verwendet, wahrend in der Praxis fiir Fliissigkeiten Ultrazentrifugen mit zylin- — 


drischen Aussparungen in einem Vollzylinder bentitzt werden, bei denen zum 
vorneherein statische Druckverhaltnisse angenommen werden diirfen. Die Druck- 


verteilung in Gas-Ultrazentrifugen erhalt man analog wie im friiheren Fall aus” 


Beziehung (1), indem man jetzt fiir @ nicht einen konstanten Wert, sondern den 


aus der Gasgleichung entnommenen einsetzt. Es ist e = M p,/RT (M = Mole- | 


kulargewicht) und somit 
Chey (06 


M 
RT Deg Het 
Hieraus folgt 

Px a Po ew 42/2 RT (8) 


(Py = Druck in der Achse). Diese Beziehung lasst sich bei Gasgemischen auch 
auf die einzelnen Partialdrucke anwenden. Demgemass ist fiir die Trennung von 
isotopen Gasen mit den Molekulargewichten M und M’ der Quotient 


eM wo? 7/2 RT 


PES Ea eh IES e(M - M’) w? f° /2 RT 
eM’ wr /2RT 


massgebend. Bei solchen Versuchen hat sich nun das Auftreten von unkontrol- 
lierbaren Gasstromungen herausgestellt. Es gelang zwar, solche durch Beimi- 
schung von H, herabzusetzen2). Doch scheint es angezeigt, einmal zu priifen, 
ob die fiir statische Verhiltnisse giiltige Beziehung (8) iiberhaupt erfiillt ist. 

Es sei daher im folgenden eine einfache Methode angegeben, welche die Prii- 
fung der Zulassigkeit von (8) erlaubt. In den Rotor wird ein einheitliches Gas 
von beliebigem Druck pa bzw. der Dichte oy eingefiillt. Man misst nun die Druck- 
erniedrigung, die in der Achse bei der Rotation entsteht. Dies geschieht sehr 


*) G. Hertz und E. Nann, Z. Elektrochem. 58, 612 (1954). 
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einfach durch Anbringen einer axialen kleinen Funkenstrecke F (Figur 3), 
Insoweit das Funkenpotential V proportional dem Druck verlauft, hat man 
einfach p,/p4 = V,/V4. Besser noch wird man V fiir das verwendete Gas in 
Funktion des Gasdruckes bestimmen und p,/p4 dann aus einer Eichkurve ent- 
nehmen. f,/p4 ergibt sich folgendermassen: Es ist die Gasmenge im ruhenden 


Figur 3 
Prifung und Anwendung der Druckformel fiir rotierende Gase mit Hilfe des Funkenpotentials. 


und im rotierenden Zustand dieselbe. Im ersten Falle ist 


4 M 
m=a2v? hog =a? h 7 Pa » (9) 
im zweiten 
r a 
M 
m= | 2axhdxo= [2nvhdx Oe by. 
0 0 


Setzt man fiir p, seinen Wert aus (8) ein, so ergibt sich 


= Po Fiber ecaae Dy 1) , (10) 


Durch Gleichsetzen von (9) und (10) erhalt man unter Verwendung der Abkiir- 
gung @= VM w*r2/2 RT 


Pas BOT (11) 
Po if 

Fiir kleine Werte von a kann man in erster Naherung schreiben 
Digweed “(a 


Wahlt man beispielsweise als Fiillung N, (M = 28, T = 300°, y= 5cm und 
w = 2- 500s-} [also fiinfmal grésser als oben]), so erhalt man fiir a = 0,139, 
~ und nach (11) wird p4/f) = 1,073. Das heisst, es entsteht eine Druckerniedrigung 
von 6,8%. Es ist also bereits bei Tourenzahlen, die noch wesentlich unter den bei 
Ultrazentrifugen angewendeten liegen, ein gut messbarer Druckeffekt zu er- 
- warten, Gleichung (11) lasst sich somit leicht experimentell priifen, zumal hier 
nicht, wie bei Fliissigkeitsfiillung, eine Deformation des Rotors berticksichtigt 
zu werden braucht. Sofern sich die Anwendbarkeit der Formel (11) ergibt, kann 
diese auch zur Bestimmung des Molekulargewichtes von Gasen verwendet werden. 
Ausser der Tourenzahl waren nur noch Spannungen zu messen. 
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Auf alle Falle lasst sich die Anordnung von Figur 3 dazu verwenden, um die 


Gaszentrifugierung bzw. die Druckabnahme in der Achse des Rotors zu demon- 
strieren. Hierzu kann irgendein transparentes Gefass von drehrunder Form mit 
einer axialen Funkenstrecke dienen. An diese wird man zum Beispiel eine Induk- 
torspannung anlegen, die zur Ziindung des Funkens gerade noch nicht ausreicht. 
Tritt nun bei zunehmender Rotationsgeschwindigkeit eine Druckabnahme ein, 
so wird man dies unmittelbar am Einsetzen der Entladung erkennen. Zur Erzie- 
lung eines Effektes bei kleinerer Tourenzahl kann man statt Luft ein schweres 
Gas, zum Beispiel Xenon, verwenden. 


Résumé 


La répartition de la pression et la compression sont calculées pour un liquide 
enfermé dans un vase cylindrique en rotation. Pour la diminution du volume on 
trouve AV = k/4 m, w2 72, supposé que toute la masse posséde une vitesse angu- 


laire uniforme. Cette condition se vérifie en comparant le AV calculé avec le © 


AV mesuré. D’autre part le coéfficient de compressibilité k se détermine en ~ 
mesurant AV. De suite la répartition de la pression dans un centrifugeur 4 gaz — 


est envisagée. La formule déduite sous la supposition d’un mouvement purement ~ 


rotatoire, vérifiée par l’expérience, peut servir a déterminer le poid moléculaire 
des gaz. Pour ce but on introduit a l’axe du centrifugeur deux électrodes qui 
permettent 4 mesurer le potentiel explosif. Il suffit d’exécuter deux mesures, 
une pour le gaz en repos et une pour le gaz en rotation. 


(Eingegangen: 10. September 1955.) 


Zur Theorie der dreidimensionalen Instabilitat 
laminarer Grenzschichten 


Von GUNTHER HAMMERLIN, Karlsruhe’) 


EKinleitung 


In dieser Arbeit soll iiber den Einfluss einer Vernachlassigung berichtet wer-_ 


den, die bei der Aufstellung der Stérungsdifferentialgleichungen der dreidimen- 
sionalen Instabilitat laminarer, inkompressibler Grenzschichten gemacht wurde. 


Wir sprechen dabei von den Stérungen langs gekriimmten Wanden, die GOrT-_ 


LER [1]?) im Jahre 1940 ansetzte. Es handelt sich dabei um Wirbel, deren Achsen 
in Richtung der Grundstrémung verlaufen; im Gegensatz zu den in Stro6mungs- 
richtung fortschreitenden Tollmien-Schlichtingschen zweidimensionalen Schwin- 
gungen sind dies stehende Wirbel, ahnlich den bekannten Taylorschen Wirbeln, 
die zwischen zwei rotierenden koaxialen Zylindern auftreten. 

Der mathematischen Formulierung des Problems liegen die Navier-Stokes- 
schen Differentialgleichungen fiir eine Stroémung langs gekriimmter Wand zu- 
grunde, in die der Ansatz dreidimensionaler St6rungen eingefiihrt wird. Diese 
Gleichungen vereinfachen sich zundchst durch die Voraussetzung, dass sich die 


) Diese Mitteilung ist im wesentlichen wahrend der Arbeit des Verfassers am Mathematischen 
Institut der Universitat Freiburg i. Br. entstanden. 


2) Die Ziftern in eckigen Klammern verweisen auf das Literaturverzeichnis, Seite 163. 
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ungestorte laminare Grundstr6mung nur wenig langs eines betrachteten Wand- 
stiicks andern soll und dass die Instabilitat ebenfalls als von der Wandbogenlange 
unabhangig angenommen werden darf. Einige Vernachlassigungen, die man 
neben dem Linearisieren des Problems vornehmen kann, beruhen darauf, dass 
nur solche Wande betrachtet werden, deren Kriimmungsradius gross gegeniiber 
der Grenzschichtdicke ist, und aufder Annahme, dass die angesetzten Stérungen 
sich nicht wesentlich tiber die Grenzschicht hinaus erstrecken sollen. 

Das so entstehende System von Differentialgleichungen stellt unter den dem 
Problem gemassen Randbedingungen eine Eigenwertaufgabe dar. Darin haben 
die Eigenwerte die Bedeutung einer Zahl Re(6/R)!/2, wobei 6 die Grenzschicht-. 
dicke, R den Kriimmungsradius der angestrémten Wand (R wird positiv fiir 
konkave Wande gewahlt) und Re die mit der Grenzschichtdicke 6 gebildete 
Reynoldssche Zahl bezeichnen. Es interessieren insbesondere die kleinsten Eigen- 
werte in Abhangigkeit von der Wellenlange 4 der angesetzten Instabilitat und 
unter diesen vor allem diejenigen neutraler StGrungen — das heisst solcher St6- 
rungen, die weder gedampft noch angefacht werden —, welche die «kritische» 
Kurve liefern. Jede St6rung der angesetzten Art wird gedampft, wenn der ihr 
zugehorige kleinste Eigenwert Re(6/R)! unterhalb des kritischen Werts liegt; 
liegt er dariiber, so wird sie angefacht und kann schliesslich zum laminar-turbu- 
lenten Umschlag fiihren. 

Nachdem man mit Hilfe von Satzen aus der Theorie der Integralgleichungen 
nachgewiesen hat, dass ein kleinster, positiver Eigenwert des Problems existiert — 
das heisst, dass die dreidimensionale Instabilitat an konkaven Wanden auftritt -, 
kann man auf verschiedenen Wegen eine strenge Lésung finden [2]. Diese Lésung 
liefert die kleinsten Eigenwerte neutraler und angefachter Stérungen wie auch 
die zugehérigen Eigenfunktionen. 

Es zeigt sich, dass die kritische Kurve Re(d/R)!, aufgetragen tiber o = a 6 
(a = 22/4), einen endlichen Grenzwert fiir o > 0 und einen fiir grosser werdende 
o monoton wachsenden Verlauf hat (vgl. in Figur 1 die Kurve mit kurz ge- 
stricheltem Ast). 


of 02 03 
Figur 1 
Die verbesserte kritische Kurve. 
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Bei der Berechnung der zum kleinsten Eigenwert neutraler Storungen geho- 


renden Eigenfunktionen stellt sich nun heraus, dass die als Stérung in Richtung ~ 


der Geschwindigkeitskomponente senkrecht zur Wand angesetzte Instabilitat mit 
kleiner werdendem o immer weiter tiber die Grenzschicht hinausgreift; damit wird 
im Bereich kleiner o der Rahmen gesprengt, der durch die Annahmen gesteckt 
worden war, die der Theorie zugrunde liegen (vgl. [2], Figuren 3 und 4). Die da- 
malige Fragestellung zielte ausschliesslich auf eine konsequente und vollstandige 
mathematische Behandlung der aufgestellten St6rungsdifferentialgleichungen 
hin, und in der Folge taucht nun die Frage nach einer Erganzung jener Ergebnisse 


durch Einschrankung gewisser Annahmen im Hinblick auf die Beriicksichtigung ~ 
der oben genannten Tatsache auf. Es soll das Ziel dieser Arbeit sein, eine ver- . 


feinerte Betrachtung fiir den Bereich kleiner Werte des Parameters o vorzunehmen. 


Die Storungsdifferentialgleichungen 


Vorgelegt sind die Navier-Stokesschen Gleichungen fiir eine Str6mung langs © 


einer gekriimmten Wand, deren Kriimmungsradius Rk (R> 0 fiir zur Stromung 
konkave Wande) konstant sei. Auf Grund der schon genannten Annahme, dass 


sowohl Grundstr6mung als auch Stérglieder unabhangig von der Wandbogen- — 
lange x sein sollen, werden die Ableitungen nach x von vornherein gestrichen. — 


y sei die Koordinate senkrecht zur Wand, z die zu ¥ und y senkrechte Koordinate. — 


Da sich die uns interessierenden Stro6mungsvorgange nur so weit iiber die Grenz- — 


schicht hinaus erstrecken sollen, dass y/R gegen 1 vernachlassigt werden darf’), 


ergeben sich die entsprechenden Vereinfachungen bei geometrischer Entwicklung 
der Ausdriicke 1/(R — y) und R/(R — y). Wir kommen damit zu diesen Diffe-— 


rentialgleichungen: 
Uu 
4 AN —_ 5 y 
Up, + U (u, — ze) + WU, = v (yy aa eg =) ; | 
u> Ba v 

ieee ; ~ y 

COO eae = hah a eee 
i Riri), (1) 

PES Is =, z y / 

Wi UW, + WW, = —— +» (yy + W22— 3Y), 2 

@ 
v 
Vy 5 a — Ww, == (y) a 


Mit der Grundstromung u = u(y), v= 0, w= 0, p= Po(y) setzen wir an — 


U= Udy) + U*(y, 2,0), U= v*(y, 2,1), w= w*(y, 2,2), p= Poly) + p*(y,z, 2); (2) 


da auch die ungestérte Grundstromung den Gleichungen (1) geniigt, entsteht bei 
Einfithren des Ansatzes (2) 


du U * ; * 
4*  y® OP. Be ‘0 u uU | 
peel (7; STR Fo Rag nieene v (uy + ut, — =P). 
2 Uy u* + 4*2 * * 
v* + y* y* 1 0 sai * U 
a ine R WW mn Mi ee 
(3) 
wt + u* w* + wt wt Ba ge z A wy 
y Z =——-t4 py Wi yt Wes op) 
* 
U 
Gok Ae oe ee * 
ly R ~ Wz = 0. 


——— 


5 ree, aoe 
) Dabei wird -ausdriicklich gesagt, dass y ein Vielfaches von 6 sein darf. 
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Es ist uj ~ u*/y, u*,~ u*/y?, und damit verhalt sich ux/R zu ux, wie v/R 
zu 1. Nach unserer vorhin gemachten Annahme darf also ut/R gegeniiber u*, 
vernachlassigt werden. Ebenso verfahren wir auf der linken Seite der ersten er 
Gleichungen (3) mit w*/R gegeniiber w*. Nachdem man in allen Gleichungen (3) in 
dieser Weise vorgegangen ist und gleichzeitig linearisiert hat, bleibt das System 


du u 
* yk oO pal 1S gen 
Uy + U ( Poe _ R ) = y (u%, + Ux.) 
u Dm 
* 0 4)* at. ] ; 
De aa = 2h to iy + Uh)» 
(4) 
* 
w% Pag ok + »(w*, + w*,) 
¢ 0 yy 22) » 
Qe 
a * ae 
Vi + we =0. 


Die anzunehmenden Grenzschichtgrundprofile u(y) der ungestérten Laminar- 
str6mung sollen so beschaffen sein, dass der Anstieg vom Wert Null an der Wand 
bis zum vollen Wert der Aussenstr6mung innerhalb der Grenzschicht erfolgt. 
Damit ist du,/dy ~ u,/d in yS 6 und du,/dy = 0 fiir y> 6. Greift die angesetzte 
St6rung v* nicht wesentlich iiber die Grenzschicht hinaus, spielen sich also die 
Stdérungsvorgange in y< 6 und damit in dem Gebiet ab, wo du,/dy gegeniiber 
u,/R weitaus tiberwiegt, kann man u,/R vernachlassigen. Das so entstehende 
System wurde bisher behandelt und ist die Grundlage der Theorie in [2]. Wie 
sich aber dort zeigt, erfiillt v* fiir kleine Parameterwerte o diese Forderung nicht 
mehr. Wir werden im folgenden diesen Parameterbereich unter Beibehaltung des 
Glieds u,/R untersuchen. 

Der St6rungsansatz i 
u = U,(y) Cosa z ePt 
v = 0,(y) cosa zeP', 
w = w,(y) sina ze! , 

p = Py(y) cosa ze", 
a = 2a/A (A die Wellenlange der Storung) , 
B = Dampfungskonstante 


. . ae . . 
fiihrt nach Eliminieren von /, und w, und nach Einfiihren der dimensionslosen 


_ Groéssen 
2 2 2\1/2 2 
n= U= o=ao, ra (SEO pa 2(—2*) 
i, O\o2 a be sven O3 
Be ca ley eg 
auf 
out tun (S70), vy — 2(r2?+ 62%) 07+ tot v=—porUu (5) 
1] 


mit den bekannten Randbedingungen 
u(0) = v(0) = v’(0) = 0, u(oo) = v(00) = U.(co) = 0. 
Fiir neutrale Stérungen wird daraus (f = 0, t = 0): 
dU 


wt —otu=(G-—7u)y, vY — 2o%v" + o4u=—po? Un. (6) 


- Diese Gleichungen unterscheiden sich von den entsprechenden friiher behandelten 


5. is te 


dadurch, dass das Glied — y U v hinzugekommen ist. 
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Bestimmung des kleinsten Eigenwerts fur kleine 0 


Die folgenden Betrachtungen zur Lésung des Systems (6) haben von vorn- 
herein das eng begrenzte Ziel, den kleinsten Eigenwert 1, in Abhangigkeit von 
kleinen Werten des Parameters o und insbesondere sein Verhalten fiir o > 0 zu 
ermitteln. Die Idee liegt nahe, hier einen Versuch mit Potenzreihenansatzen fiir 
u(n) und v(y) zu machen, wobei in den Koeffizienten dieser Potenzreihen der 
Parameter o erscheinen wird. Im Bereich 0S 7S 1, wo sich eine exakte Loésung 
von (6) nicht angeben lasst, soll dieser Weg beschritten werden. Anders liegen 
die Verhaltnisse fiir 1< , wo eine exakte Lésung von (6) bestimmt werden 
kann. Wir werden von den Potenzreihen der Innenlésungen zu verlangen haben, 
dass sie die Randbedingungen bei 7 = 0 erfiillen und bei 7 = 1 mit der exakten 
Aussenlésung iibereinstimmen sollen. Diese Zusammenschlussbedingung wird uns 
eine Bestimmungsgleichung ftir “(o) liefern. 

Fiir 7 2 1 wird aus dem System (6) [U(y) = 1): 


Lful=—yo, L{v]=—potu, L-3,-# | (7) 
fe(col = 05 Vut(co) =O, 2c) Ol | 


Also ist 
Liv] = yworv 


mit den Fundamentallésungen e°”” 


Qre= + (0? + Vy no?) ?, ogg & (02 + 07718 Vy py 02)", 
sg = & (0? + 88 Yor ot)? 


: Da die Randbedingung v(co) = 0 zu erfiillen ist, setzt sich die gesuchte 
Lésung aus denjenigen Fundamentallésungen zusammen, deren Realteil kleiner 
als Null ist. Das ist wegen «> 04) zunachst fiir g, der Fall. . 

Die Wurzeln g, bis g, sind zu je zweien konjugiert komplex. Nach einfach ew 


Umformung ergeben sich als Fundamentallésungen mit abklingendem Exponen- 
tialfaktor “- 


Be a : 
Tcosax n, v3(y) =e 420 sin Ae . 


r= + (o?+ e212 (2 = Vy o? p) : 


1 } 
2 V(o? — 6%)? + 62 e@ + - (2 ie e*)]" (8) 
il 


04(7) = Le V9(q) = e- 


i Selbstverstandlich erfiillt die Lésung v,(7), die sich aus den Fundamental- 
i Sa (8) Zusammensetzt, auch die Randbedingung v’(co) = 0. Man kénnte 
nee * ee Pee bat: der Gleichungen (7), u4(7) aus v,(7) berechnen, 
ne weiteres, dass diese Lésung u,(n) di i 7 

> mont =. er ) n) die Randbedingung 
u(oo) 0 befriedigt. Die drei im Unendlichen zu erfiillenden Ravid bedi ee 
ae ee, 


4) In [2] wurde u > 0 gezeigt; be 


alsoraibh tsier ss = 0. = 0 kann nicht eintreten, das bringt nur die triviale Lésung ; 


| 


€ 


nS eB eed Ce 
_ 
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lassen also von den sechs urspriinglich verfiigbaven Konstanten noch drei un- 


bestimmt. 
Die Lésung v,(7) geniigt den Differentialgleichungen 


UN + Gg" + av’ + agv=De w+ ayy” + av" +a,v'=0, 
oY F-ac ye au" 4 a)0" = 0 
mit 
— 2 
Ay = A, (AS + AF?) = 3 2+ a 
@, = AB+ AFF + 24,4, 
{ 7 o2\2. oh , o2 o2\2 2791/2 os om 
oe Aa Mee | et > Ge | o it 
. V/( =] e2 /2 (1 + >) VIG <=) =P r| e2 a 1: 
ine) 24, ae | r+ te yzy/{(i BE igh odd OP data wal Tk 
A 1 2 | se EM call | =| res a ea 


Unter Zuhilfenahme des Systems (7) lassen sich daraus folgende Gleichungen 


herleiten, in denen nur w(n), u3(n), Va(y) , ..., vg (n) auftreten: 


(a, — a8 + 2 o?) uv” + (a) — a, ay) v’ — (a, a,+ ot) —poru 
(a, — a3 + 2 o?) vu" + (ay — ay a.) v” — (Ay a, + Go) v0’ — wo? wu’ 
oy + a + av + av 


I 


Von den noch unbekannten Innenlésungen u,;(7) und v;(7) fordern wir nun als 
_ Zusammenschlussbedingung, dass sie bei 7 = 1 diesen Differentialgleichungen 
geniigen sollen. Im Hinblick darauf, dass wir uns nur fiir die Lésungen bei kleinen 
Werten des Parameters o interessieren, werden die Koeffizienten von (9) in dieser 


Form geschrieben: 


ao = a] 1+ jo (z)} (fo(0) = 0], 
a, -3{ 24+4(G)} (0 =a. 
i ~2{ 24h (Z)} 0) = 9, 


Tl 
oo 

| 
is) 
+ 


i 
ay — 4, 4 = 6 {- 3+ 1 (3) [1(0) = 0], 
ay 4, + o* = e4 = 2+ Yo (3)} [ (9) = 9] 
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A* = 
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Mit diesen Bezeichnungen wird aus (9) 
; a 
(2— ps) v" Fe S) pi) Rie OES 96) Oe ee | 
te 10 
(2 pa) OE eS gah OC ete Pa) 7 erates eae | i 
ot e(2t flout ett hl’ t elt hjv—o0. 


Von den Betrachtungen in [2] her ist pekannt, dass die spezielle Wahl der 


Grenzschichtgrundprofile keinen wesentlichen Einfluss auf den Verlauf der kri- — 


tischen Kurve hat, wenn wir alle Profile auf gleiche Impulsverlustdicken be- 
ziehen ([2], S. 314). Wir kénnen das auch hier annehmen und uns darauf be- 
schranken, mit dem Streckenprofil 


y fiir 


0<n<1 
U(n) = : 


ibe ale<y)) 


zu arbeiten. Dann lauten die Differentialgleichungen (6) in 0< 7 <1, aus denen 
die Innenlosungen hervorgehen: 


u”“— Pu=(1l—yyn)v, vV — 262v"+ ctv =—po?nu | (11) 
(a(S 0 w(0) = 04! voy Store hast 
Potenzreihensatz: 
u(n) == a2 MEE CRE AAA n\) = aie Bs 2 Z 1 
te Oy OT es Ole ek ene eee (12) 


Die Randbedingungen werden erfillt durcha,=0, Bp) = 0, By =0. 

Lasst man «,, 6, und f, noch unbestimmt, so sind die iibrigen Koeffizienten 
der Potenzreihen in Abhangigkeit von diesen auszudriicken. Die Forderung, 
dass die Lésungen (12) bei 7 = 1 das Gleichungssystem (10) erfiillen sollen, fiihrt 
auf ein homogenes Gleichungssystem fiir die Koeffizienten a,, B, und f,; die 


Determinante dieses Systems muss also Null werden. Diese Determinante er- 
rechnet sich zu 


A = — A*, 
ie 
jo era i cise at ae ee 2p! ae 
3! ue alae 812 ayig ee 348 ea ie eae 81s 
4 1 4 
é 24 E 1 4 
ay oak eta Pane ate Sages 2436+ 8+ Set oe 


Zeer 2er+ gos 1+ 26+ €%+ Bo3 


66 
gis = Pix (0% + , 7), Polynome. P,,(0, 0, 0) =0. 


A = Oist zunachst erfiillt fur e4/y = 0, das heisst fiir 4 = 0. Dieser Wert scheidet 
aus, da er nur die trivialen Lésungen u = 0 und v = 0 mit sich bringt. 
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nee Berechnung méglicher Eigenwerte sind also die Nullstellen von A* aufzu- 
suchen. 


4 : 
teksts Airs a Call Ee ae 101 
ah ie na oer ee kaa © 


1 Ege \2 1 8 
a: at a1 at ow o(24 e+ 5 et) + @ (a2 ; y) = 0. 


Man findet die kleinste Nullstelle unter der Voraussetzung o < 1 angendhert zu 
e=2y'; ® bleibt unberiicksichtigt, da es nur die hGheren Potenzen von o und 
e* enthalt, wahrend y von der Gréssenordnung 10-4 ist). 
Also ist 
lime = 2 yt!4 | 
o—>0 


das bedeutet 


uu 6 \1/2 32 yl 
2 = Re (>) we 2M er oft)}; 


bezogen auf die Impulsverlustdicke #@ (5/4 = 6 beim Streckenprofil) : 


U,8 1/9 dp. 
a8 Yije © 00064 2 (1 + o(1)}. 


Von der kritischen Kurve kénnen wir uns also dieses Bild machen: Die hier 
gewonnene asymptotische Darstellung ist sicher brauchbar fiir o < 0,05, wah- 
rend der Einfluss des in (6) neu hinzugenommenen Glieds fiir etwa o > 0,4 keine 
Rolle mehr spielt. Verbinden wir diese Kurvendste wie in Figur 1, so besitzt die 
damit gewonnene kritische Kurve ein Minimum bei o = 0,33; das entspricht 
einer St6rungswellenlange 2 = 19 6. Die Daten in Figur 1 beziehen sich auf die 
Impulsverlustdicke #, das Minimum der kritischen Kurve liegt bei « # = 0,055 
(y = 10-4). Fiir kleinere Werte von y liegen die Minima entsprechend weiter links. 

Diese Ergebnisse erganzen die Theorie der dreidimensionalen Instabilitat. 
Die kritischen Werte werden bestatigt, ebenso die Tatsache, dass die kritischen 
Wellenlangen ausserordentlich gross sind; das ergab sich in [2] aus der Betrach- 
tung der Kurven geringer Anfachung. Der Einfluss des hier mitberiicksichtigten 
Glieds macht sich nur fiir sehr kleine o bemerkbar und andert nichts Wesentliches 
an dem Gesamtbild, das sich in [2] ergab. 


LITERATURVERZEICHNIS 


[1] H. Gortier, Uber eine dreidimensionale I nstabilitat laminarer Grenzschichten 
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5) Vel. die in [2], S. 317, zitierten Versuchsdaten von H.W. LI=PMANN. 
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Summary 


In the case of large wave lengths a refined treatment of the three-dimen- 
sional instability of laminar boundary layers is made possible by taking into 
account a further term of the differential equations of the instability. The 
improved curve of the Gértler-Parameter Re (#/R)1/? for neutral disturbances has 
a minimum with an instability wave length approximately equal to 19 6, 6 the 
boundary layer thickness. The results of [2] are not changed essentially, especially 
the critical values of the Gértler parameter are confirmed and also the fact, 
that the critical wave length is large. 


(Eingegangen: 7. Oktober 1955.) 


Eine Forme! von Wronski 
und ihre Bedeutung fiir den Quotienten-Differenzen-Algorithmus 


Von Heinz RuTISHAUSER, Ziirich?) 


1. Einleitung 


In einer friiheren Arbeit ([3], Abschnitt 6) hat der Verfasser gezeigt, wie man 
die Nullstellen eines Polynoms N(z) (vom Grade m) mit Hilfe des QD-Algorithmus 
bestimmen kann. Man baut dabei das QD-Schema einer rationalen Funktion 
f(z) = N,(2)/N(z) [mit einem beliebigen Polynom (m — 1)-ten Grades als Zahler] 
von oben nach unten auf, wobei man die fiir den Beginn des Prozesses notwendige 
oberste Schragreihe, das heisst die Werte qe), el), gq’, s Panes el) re qs) aus der 
S-Kettenbruchentwicklung von f(z) entnehmen kann. Es ist natiirlich ein Nach- 
teil des Verfahrens, dass zuerst eine Funktion f(z) in einen Kettenbruch ent- 
wickelt werden muss. 

Es diivfte daher von Intevesse sein, dass man die Kettenbruchentwicklung ver- 
meiden kann, indem sich das QD-Schema der Funktion z"-1/N(z) auf Grund einer 
Formel von WRonskKI in besonders einfacher Weise aufstellen léisst. 

Zu dieser Vereinfachung wurde der Verfasser durch eine Arbeit von A. C. 
AITKEN [1] angeregt, in der eine entsprechende Vereinfachung des Aitkenschen 
Algorithmus behandelt wird. Den Hinweis auf die Formel von WRONSKI ver- 
danke ich Herrn Dr. F. L. Bauer, Technische Hochschule Miinchen. 


2. Die Wronskische Formel 


co 
Seley (7) k= 2 a, x eine Funktion, die in einer gewissen Umgebung des Null- 


punkts regular sei. Dann ist (falls ay + 0) die Potenzreihe der reziproken Funk- 


tion 1/f(¥) durch folgende Reihe gegeben (siehe WRONSKI [6]; WHITTAKER [4], 
Abschnitt 60): 


OR ania 
i: if y 1 a, ie | ao ay x2 fy 01 3 
i ' => By Gq Ge ba ee 
f(*) ay ag os a,| as sien ar its (1) 
| 41 4p as, 


1) Institut fiir angewandte Mathematik der ETH. 
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mit dem allgemeinen Glied 


OF 0) ar OL ayaa" | 
Oe 0T or Be ty as 

ite ; 

( 2 on 

(~1) ; gerd (2) 

0 a 0 
py nw oe ay dig | 
Cag snr a4 a, 


Dies erlaubt die Berechnung der kleinsten Nullstelle von f(x) als Grenzwert des 
Quotienten aufeinanderfolgender Koeffizienten der Wronskischen Reihe (vgl. 
hierzu WHITTAKER [4] oder [5], Abschnitt 60). 

Wir wollen nun aber die Formel (1) zum Aufbau des QD-Schemas der ratio- 
nalen Funktion 


1 So Sy Sg Sn—1 
We ae ee eh (3) 
n 
mit N(z) = » c, 2* verwenden. Dabei wird natiirlich s)='s,;=s,=---=s, ,.= 0, 
0 


so dass der Anfang des QD-Schemas von 1/N(z) zunachst gar nicht definiert ist. 
Im allgemeinen wird es aber von der mit »y =  — 1 numerierten Schragzeile an 
abwarts existieren?) und kann auch noch ein Stiick weit nach oben fortgesetzt 
werden. 

Substituiert man in (3) z = 1/%, so resultiert 


es) = So%¥ +51 24 we + 5) 4 i peg BY Sg BETT eee (4a) 
= (0) 
oder 
SMTA Tt te aes (4b) 
oder 
im N (ifs) =o, + Cy) ¥ tees + og = : (5) 


+= rece 
Sn—-1 ap Sn eaG Sn+1 we 


Da nun die rechte Seite von (5) nach (1) in der Form 


2 } x3 
1 s Sn—1 Sn i | : A 
: =" ‘9 SH TT Sn-i Sn Sata) Ga + 
Sn-1 n—1 Sn Sn41i Sn—1 Se ae ee n-1 


2) Das QD-Schema von 1/N(z) ist von dieser Stelle an nichts anderes als das QD-Schema 
von z”—1/N(z). 
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geschrieben werden kann, folgt durch Koeffizientenvergleich . 3 
4 
“a+ ) al > 
. ( 2 (wm+1—“) = = = 
Cnn (=) Ty gett (x Ol ae ahserens n) ; 


wobei fiir alle » 


ae te mee 


Sy Spy. Sppx—t 


joe= 


| 
Sp44 Spie Sytse | ‘ ae ae _ a. 
Ny | (x> 0, mit sy = Sy = ++ = Spe = 0). 


Sp4n—-1 Spex Sv42%—2 


Damit haben wir ein Mittel gefunden, um die mit dem QD-Algorithmus zu-_ 
sammenhangenden Hankel-Determinanten H aa der s, explizite anzugeben. — 
Ausserdem haben wir ja 


i 
Sie eS eecein eed OR OMS OS ae. 
Sn—xt1 Sn—n+2 +++ Sn 0 O «++ Spt Sn : 
H"-” At wey 5 = 1 (2) x 7 , 
oy = ae e Fy (= ) Sn-1 ( ) 
0 oe : 
tie WS iia OSA 5 aS hah. eae eee 
und entsprechend 
| Sn—n—1 Sn—x +++ Sn—2 
1 Sn—x 
edi : =0 (¥=1,2,...,0-1) (8) 
pees tee Sin-tn—3 


3. Der Aufbau des QD-Schemas von 1/N(z) 
Fiir das Weitere beniitzen wir die Tatsache, dass samtliche im QD-Schema 
(komprimierte Form, siehe [3], Abschnitt 1) auftretenden Grdssen q”) und e) 


durch die mit den Koeffizienten s, gebildeten Hankel-Determinanten H ) aus- 


gedriickt werden kénnen. Es gilt naémlich mit den gleichen Bezeichnungen wie 
in (6) (vgl. hierzu auch [2], Abschnitt 5): 


+1) Fr) 1 
g® = sa eres pe) ey eee (9) 
Hee Se é H®) He*) 


In unserem Fall gelten diese Formeln zunachst nur fiir y= — 1, wo das QD- 
Schema von 1/N in normaler Weise definiert ist. Nun ist aber N (2) ein Polynom 
vom Grade m und daher e”) = 0 fur alle y (vgl. [2], Satz 2), so dass sich das QD- 
Schema von 1/N(z) auf Grund der Formeln (4) in [3] auch riickwAarts fortsetzen 
lasst, soweit die Gréssen q, 2”) aus den Formeln (9) sinnvoll berechnet werden 
konnen (keine Nullen im Nenner), Dies ist offenbar der Fall, solange in den im 
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Nenner auftretenden H die Summe von oberem und unterem Index = n ist; 
das heisst fiir »> n — x. 

_Somit lasst sich der im untenstehenden QD-Schema von 1 /N(z) eingerahmte 
Teil durch Fortsetzung von der mit » = » — 1 numerierten Schragzeile aus nach 
oben sinnvoll erganzen: 


= 
ge a 


Es zeigt sich nun, dass die am Kopf des so erganzten QD-Schemas stehenden 
Werte g~*) und e”~*) auf Grund der Formeln (6) bis (9) leicht berechnet werden 
k6nnen; es ist namlich: 


i 4 2 
qn) Hy” = a) Snel pod =-—C Ss — Ct 
ie rams =f \0 Bel ed ia 
Ht" (—1) SiS Cn 
Ht+1-#) A» 
a i er eae wegen Bo == 0m sei Zne3), aon dh)ie 
* H”-” H®+1-*# % 
x x-1 (10) 
5 (4 = i) = 
= $i— 2 “+2 = 
elt —*) | oe Be s (4) seer ee el eae (is Sn Cn=n-1 
x ¥ (n—x) zy(n+1—x) x atl Cc 
a ae peti tela d ett * alee 


4. Zusammenfassung 
n 
Zur Bestimmung der Nullstellen des Polynoms N(z) = »'¢, 2% stellen wir das 
0 


QD-Schema der Funktion 1/N(z) auf, indem wir mit den in einer horizontalen 
Zickzacklinie stehenden Werten 


(w—1) _ __ _©n—-1 
q as 


92) 25 Oo rePee ln By 5. HM) 5 


x” 


elt-*) — nA (ae Ie Doe dt mille 
3° % Cc 
n—x 


168 Kurze Mitteilungen — Brief Reports — Communications bréves ZAMP 


oder schematisch: 


(41) (2) (42) (29) (93) (In—-1) (€,-1) (Qn) 
on—1 0) 0, Au) 0 
Cn a ee 3 
Cn—2 Cn-38 \ i oo 
ns B 4 = 
Cn-1 Cyn—2 . a Cy 


OE  —————s 


beginnen und dann unter Beachtung von ee) = e) = 0 in der iiblichen Weise 
nach unten fortsetzen; die Formeln hierzu lauten (siehe [3], Abschnitt 5): 


os 


ge th = gq’ ae e”) ae Cee : A f (12) 
: 
) 0) ; . 
oD = Pn Tusa (13) 
“ (vy +1) 
qs, 


a 


Dabei ist es mnemotechnisch vorteilhaft, die Formeln (12) und (13) als sogenannte 
Rhombenregeln zu formulieren’) : | 

Die Formeln (12) und (13) verkniipfen immer je vier in einem Rhombus_ 
angeordnete Werte des QD-Schemas (in komprimierter Form), und zwar unter- 
scheidet man g- und e-Rhomben. : 


spate 
~~ 


Fiir die g-Rhomben: OMe ee buat | 
(v+1) os 
q. S 


Die Summe der tiber der Trennungslinie stehenden Werte ist gleich der Summe der 
beiden unteren Werte (entsprechend dey Formel (12)]. 


Sua el) 
Fiir die e-Rhomben: Got) tos Socata 
ss tx41 
ae ee SS 
en 


Das Produkt dev tibey der Trennungslinie stehenden Werte ist gleich dem Produkt 
der beiden unteren Werte [entspricht Formel (13)]. Die Nullstellen ergeben sich 
dann in der in [3], Abschnitte 6 und 7, angegebenen Weise. 


3) Nach einem Vorschlag von Prof. Dr, E. STIEFEL, ETH. 
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Numerisches Beispiel: x*— 9x2 — 8% +2=0 (wie in [3]). 


val ey Gs es G3 
++ 0,8889 — 0,25 

9,8889 <Srngs6 +025 
= 2024 -- 0,0549 

9,7865 — 0,9816 0,1951 
-+ 0,0103 — 0,0109 

9,7968 —+1 0028 0,2060 
— 0,0010 + 0,0022 

9,7958 — 0,9996 0,2038 
+ 0,0001 — 0,0004 

at 299 — 1,0001 0,2042 
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Summary 


In [3], paragraph 6, the author showed how to find the zeros of a polynomial 
with the aid of the QD-Algorithm. The method required a certain rational 
function to be developed into a continued fraction of Stieltjes type to give the 
starting values for the QD-Algorithm. In the present paper the author shows 
how starting values can be obtained without computing a continued fraction. 


(Eingegangen: 11. Juni 1955.) 
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Meeting of the Association for Computing Machinery, Los Angeles, 
August 27-29, 1956 


Annual meeting of the Association for Computing Machinery will be held’ 
on the University of California, Westwood Campus, Los Angeles, August 27-29, 
1956. For information write G. W. King, Box 3251, Olympic Station, Beverly 
Hills, California. Submit papers (abstract and four page manuscript in triplicate) 
by May 15 to J. P. Nash, University of Illinois, Urbana, Illinois. (See January 
issue Journal of Association for Computing Machinery for further details.) 3 

E. L. BRowN 


Congrés international de cybernétique 


La Province de Namur organise, sous le Haut Patronage du Ministére de 
V'instruction publique et de l’U.N.E.S.C.O., et sous la Présidence d’Honneur de 
Monsieur le Gouverneur de la Province de Namur, un congrés international de 
cybernétique qui se tiendra 4 Namur (Belgique), du 26 au 29 juin 1956. 

Les personnes qui auraient l’intention de participer au congrés au titre 
d’auditeur, d’auteur de communication ou d’exposant de matériel, sont priées de 
se faire connaitre en écrivant au Secrétariat du congrés international de cyber- 
nétique, 13, rue Basse-Marcelle, Namur (Belgique). Des informations détaillées 
leur seront envoyées. E. STIEFEL 


Erratum — | 


Corrections to the paper On Conditions for Stability of Solutions of 
Pendulum-Type Equations by James C. Litto and GEroRGE SEIFERT, 
ZAMP 6, fasc. 3, 239-243 (1955). : a 

On page 240 in the third sentence of the section entitled The Lower Bound, 
the inequalities restricting b should obviously be reversed. 

On page 241 the inequality in (5) should be reversed. | 

On page 242, line 16, the second inequality should read: y,(6) < h(0); i.e., 
a subscript ¢ should be added to the function (6) 


ey 


Buchbesprechungen — Book Reviews — Notices bibliographiques 


Millimicrosecond Pulse Techniques. Von I. A.D.Lewis und F.H.WeELts 
(Pergamon Press Ltd., London 1954). 310 S., 150 Fig.; 40s. 

Die vorliegende Monographie befasst sich mit der Theorie und dem prakti- 
schen Entwurf von Schaltungen der modernen Impulstechnik. Der iiberwiegende 
Teil der Anwendungen ist der Kernphysik entnommen, doch diirften manche der 
untersuchten Probleme auch in anderen Gebieten der Messtechnik von Interesse 
sein. Insbesondere werden Impulstransformatoren, Breitbandverstarker, «distri- 
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‘buted amplifier», Impulsgeneratoren sowie Kathodenstrahloszillographen und 
Koinzidenzschaltungen mit extrem hohem zeitlichem Auflésungsvermégen ein- 
gehend behandelt. Zahlreiche Literaturhinweise erleichtern das Studium der 
Originalarbeiten. 

In diesem Zusammenhang sei erwihnt, dass Sekundaremissionspentoden zur 
_Erzeugung von kurzzeitigen Rechteckimpulsen grosse Vorteile bieten. Es ist des- 
halb zu begriissen, dass solche Schaltungen in dem vorliegenden Werke Aufnahme 

gefunden haben. Die beiden Verfasser sind in einem bekannten Entwicklungs- 
pe ratorians (A. E.R. E. HARWELL) tatig, und ihre grosse Erfahrung zeigt sich in 
der tibersichtlichen und klaren Darstellung sowie in der sorgfaltigen Beriicksich- 
tigung der jiingsten Entwicklungen in diesem Spezialgebiet der elektronischen 
-Messtechnik. Das Buch entspricht einem Bediirfnis und darf Physikern und 
Ingenieuren sehr empfohlen werden. E. Baldinger 


Advances in Geophysics, Bd. 1. Herausgegeben von H. E. LANDSBERG 
(Academic Press Inc., New York 1952). 362 S., 65 Fig.; $7.80. 

Das Interesse an Geophysik hat in den letzten Jahren um ein Vielfaches zu- 

genommen. Neue Zeitschriften sind erschienen und bekannte haben ihren Umfang 

-beachtlich erh6ht. Entsprechend wurde in relativ kurzer Zeit eine ungeahnte 
Fille von Tatsachenmaterial zusammengetragen. Auffallend wenig wurde aber 
bis heute getan, um dieses Material kritisch zu sichten und zu ordnen und damit 
eine konsequente Weiterentwicklung zu erleichtern. Der Gedanke von LANDSBERG, 
nun auch in der Geophysik von Zeit zu Zeit eine Reihe von zusammenfassenden 
Berichten, die den jeweiligen Stand der Erkenntnisse auf Teilgebieten darstellen, 
in Buchform herauszugeben, ist daher sehr zu begriissen. 

Der vorliegende Band enthalt folgende Beitrage: 

Automatic Processing of Geophysical Data. Von JOHN C. BELLAmy. In der Geo- 
physik sind vorwiegend einmalige Vorgange zu beobachten, die oft von zahl- 
reichen Parametern abhangen, oder solche, bei denen zum Beispiel die 4usseren 
Bedingungen dauernd geandert werden. Einwandfreie Schliisse k6nnen daher nur 
‘aus einem umfangreichen Beobachtungsmaterial gezogen werden. In diesem Be- 
richt wird nun dargelegt, wie die neuen Hilfsmittel, wie Tonbander, Mikrofilme, 
Lochkartenmaschinen, Analogierechengerate und elektronische Rechenmaschi- 
nen, mit Erfolg zur Verarbeitung eines umfangreichen Beobachtungsmaterials 

eingesetzt werden kénnen. Die Darstellung ist insofern etwas einseitig, als nur 
Beispiele aus der Meteorologie und insbesondere die Messungen mit Radiosonden 
behandelt werden. 

Some New Statistical Techniques in Geophysics. Von ARNOLD Court. Dieser 
Beitrag behandelt vor allem einige neue Gesichtspunkte der Statistik, wie zum 
Beispiel die Kreisverteilungen unter besonderer Beriicksichtigung der Belange 
der Geophysik. Die grundlegenden Theorien sowie Beweise oder illustrative Bei- 
spiele, die auch die Grenzen der Anwendungsméglichkeiten leichter erkennen 
lassen, sind weggelassen. Die Darstellung ist sehr klar, und vom Leser werden 
‘nur einige Vorkenntnisse verlangt, wie zum Beispiel die Methode der kleinsten 
Quadrate, einfache Korrelationen usw. 

Studies of the General Circulation of the Atmosphere. Von Bert Borin. Hier 
wird das Verhalten der Atmosphare im Grossen betrachtet. Es werden insbe- 
sondere folgende Fragen untersucht: Die mittlere Bewegung und ihre Stérung 
infolge der unregelmdssigen Erdoberflache im Zusammenhang mit den physika- 
lischen Ursachen, die Energiebilanz und die beiden bekannten Modelle der Atmo- 


-sphare. 
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a 
Explovation of the Upper Atmosphere by Meteoritic Techniques. Von Frep L, 
WHIPPLE. In diesem Beitrag werden vorerst die verschiedenen Beobachtungs- 
methoden der Astrophysik erlautert und dann vier verschiedene Methoden zur 
Bestimmung der Dichte der héheren Schichten der Atmosphare aus photogra- 
phischen Beobachtungen von Meteoren dargelegt. Bemerkenswert ist, dass alle 
Methoden annahernd gleiche Resultate liefern. ; 
Unsolved Problems in Physics of the High Atmosphere. Von N.C. GERSON, 
Dieser Bericht gibt vorweg eine Zusammenfassung der heute bekannten physi- 
kalischen Eigenschaften der Atmosphare. Die Messtechnik sowie Probleme im 
Zusammenhang mit der Nachrichteniibermittlung werden tibergangen. Im ganzen 
gesehen, steht hier die Forschung noch im Anfangsstadium. ; : 
Estuarine Hydrography. Von D. W. PritcHARD. Der Verfasser berichtet im 
wesentlichen iiber die neuen Untersuchungen, die in der Chesapeake Bay und im 
James River durchgefiihrt wurden. Ferner werden die Theorien von STOMMEL 
und CAMERON dargelegt. ; 
The Earth’s Gravitational Field and Its Exploitation. Von GEORGE Prior WooL- 
LARD. Nach einem historischen Uberblick wird die Messung und Berechnung der 
Erdbeschleunigung kurz behandelt. Hierauf wird das Schwerefeld der Erde be- 
schrieben und die Ursachen von Schwerestérungen und die Genauigkeit von 
Schweremessungen diskutiert. Weiter wird auf die praktische Bedeutung der 
gravimetrischen Messungen fiir die Geodasie eingegangen, wahrend die eigentliche 
angewandte Gravimetrik kaum gestreift wird. Es ist daher zu hoffen, dass dieses” 
wichtige Teilgebiet, in dem in den letzten Jahren bedeutende Fortschritte erzielt 
wurden, in einem weiteren Beitrag ausfiihrlich bearbeitet wird. . 
Aeromagnetic Surveying. Von JAMES R. BAtsLey. In diesem Aufsatz werden 
vorerst die heute gebrauchlichen Messgerate, deren Empfindlichkeit und Genauig- 
keit besonders hervorzuheben ist, sowie die verschiedenen Hilfsgerate, wie Spe- 
zialkameras fiir Luftaufahmen, elektronische Navigationsgerate, Altimeter usw., 
beschrieben. Dann folgt eine Darstellung der eigentlichen Messtechnik und der 
erzielten Ergebnisse, die ohne Zweifel, wie allein das beigegebene Kartenmaterial 
zeigt, als hervorragend bezeichnet werden miissen. 
Allgemein ist noch zu sagen, dass allen Beitragen ein ausfiihrliches Literatur- 
verzeichnis beigegeben ist. M. Weber 
Hohere Mathematik, Teil V: Formelsammlung. Von G. DoBBRACK (Ver- 
lagsgesellschaft B.G.Teubner, Stuttgart 1954); 124 S., 74 Abb.; DM 4.80. 
Nachdem die meisten Bande des Lehrwerks von R. Rotue in neuer Auflage 
herausgekommen sind, erscheint nun auch die zu den drei ersten Banden (im. 
wesentlichen: Differential- und Integralrechnung, Elemente der Differential- 
geometrie, gewohnliche Differentialgleichungen reeller Veranderlicher) gehdrende 
und diese zusammenfassende Formelsammlung zum dritten Male. Es sind eine 
Reihe von kleinern Anderungen vorgenommen worden, einerseits um die Zusatze 
der iibrigen Teile zu beriicksichtigen, andrerseits um die Brauchbarkeit zu er- 
hohen. Es ist nicht zu bezweifeln, dass die Formelsammlung dem Studierenden 
wie auch dem Praktiker niitzliche Dienste leisten wird. E. Roth-Desmeules 


Berechnung magnetischer Felder. Von Franz OLLENDORFF (Springer- 
Verlag, Wien 1952). 432 S., 287 Abb.; Fr. 67.50. 

Das Buch setzt die Kenntnis der Maxwellschen Gleichungen sowie eine sehr 
gute mathematische Schulung voraus. Zweck des Buches ist die mathematische 
und zahlenmassige Erfassung von magnetischen Feldern, soweit diese fiir die 
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Elektrotechnik von Interesse sind. So werden beispielsweise behandelt: Vollpol- 
und Schenkelpolmaschinen, auch mit genutetem Auker, die Ankerriickwirkung, 
die Stirnverbindungen, das Nutenproblem vermittels der Schwarz-Christoffel- 
schen Differentialgleichung, der Wendepol, die Jochstreuung von Transforma- 
toren, die Induktivitat technischer Mehrfachleitersysteme, das Nutenstreufeld, 
die Streufelder von Scheiben- und Zylinderwicklungen bei Transformatoren, die 
Induktivitat von Kafigringen, der Zugmagnet, der Drehmagnet, magnetische 
Abscheider, Messwerke und elektrodynamische Krafte bei Sammelschienen, 
Trennern und Kontakten. 

Zur Bewiltigung dieser Aufgabe werden alle einschlagigen mathematischen 
Methoden herangezogen, wie skalare Potentiale und Vektorpotentiale, geradlinige 
und krummlinige Koordinaten, die durch analytische Funktionen vermittelten 
konformen Abbildungen fiir ebene Probleme, ferner die Kugel- und Zylinder- 
funktionen, Hankelsche und Besselsche Funktionen, Jakobische elliptische Punk- 
tionen, Fourier-Reihen usw. Es wird also eingehend Gebrauch gemacht von dem 
mathematischen Riistzeug des Physikers. Die Berechnungen werden so weit 
durchgefiihrt, dass die zahlenmissige Auswertung méglich wird. Auch die gra- 
phische Darstellung wird eingehend gepflegt. 

Das Buch ist ausgezeichnet geschrieben und behandelt die Probleme von 
hoher Warte aus. Es enthalt die erstmalige Lésung verschiedener Probleme; auch 
der Druck ist hervorragend. Das Buch wendet sich an Mathematiker, Physiker 
und theoretische Elektroingenieure. Ed. Gerecke 


Studies in Mathematics and Mechanics. RICHARD von MIsEs gewidmet, 

von Freunden, Kollegen und Schiilern (Academic Press, Inc., Publishers, New 
York 1954). 353 S., 16 Fig.; $9.00. 
_ Am 19, April 1953 hat RricHARD voN MiIsEs seinen 70. Geburtstag gefeiert; 
am 14. Juli 1953 hat ihn der Tod ereilt. So ist die als Ehrung des Siebzigjahrigen 
gedachte, von G. BirKHOFF, G. KuERTI und G. SzEcG6 veranlasste und zusam- 
mengestellte Sammlung von Arbeiten, welche die Vielseitigkeit des grossen 
Forschers und Lehrers aufs eindriicklichste dokumentiert, zum Gedenkband 
geworden. 

Der Band enthdlt eine kurze Wiirdigung von RICHARD von MisEs’ Pers6nlich- 
keit und Werk aus der Feder von Pu. FRANK sowie ein vollstandiges Verzeichnis 
‘seiner Publikationen mit ihren bekannten Héhepunkten in: Die Differential- und 
Integralgleichungen der Mechanik und Physik; Wahrscheinlichkeit, Statistik wnd 

Wahrheit; Wahrscheinlichkeitsrechnung und thre Anwendungen in dey Statistik und 
theoretischen Physik; Theory of Flight und Positivism, a Study in Human Under- 
standing. Als besonderes Charakteristikum dieses Verzeichnisses sei erwahnt, dass 
es auch eine Anzahl von Ver6ffentlichungen iiber RAINER Marta RILKE enthalt. 
Die R1cHARD von Mises gewidmeten Arbeiten verteilen sich auf die Gebiete 
Algebra und Zahlentheorie, Geometrie, Analysis, Theoretische und Angewandte 
Mechanik, Wahrscheinlichkeit und Statistik. 

Die Beitrage aus dem Gebiete der Mechanik sind: J. M. BurcERs, Some 
Considerations on Turbulent Flow with Shear; P. R. GARABEDIAN, An Example 
of Axially Symmetric Flow with a Free Surface; G. KUERTI, On a Class of Spheri- 
cally Symmetric Flows; C. C. Lin, On Periodically Oscillating Wakes in the Oseen 
Approximation; CH. LOEWNER, Some Bounds for the Critical Free Stream M ach 
Number of a Compressible Flow Around an Obstacle; G.S.S.LUDFORD, Two Topics 
in One-Dimensional Gas Dynamics; N. Minorsky, On the Stroboscopic Method ; 
G.Pétya, Estimates for Eigenvalues; W.PRAGER, On Slow Visco-Plastic Flow; 
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J.L.Syvneex, Relativistically Rigid Surfaces; P.W.BripeMan, Certain Aspects of 
Plastic Flow Under High Stresses; H.W.Emmons, Natural Convection Heat Trans- 
fer Correlation; K, FEDERHOFER, Dev senkvecht zu seiner Ebene belastete, elastisch 
gebettete Kreistriger; A.M. FREUDENTHAL, On Inelastic Thermal Stresses; H. and 
E. REISSNER, Torsion of a Circular Cylindrical Body by Means of Tvactions’ 
Exerted Upon the Cylindrical Boundary; A.S1cNortNt, Una espressiva applica- 
zione della proprieta di media dello stress communi a tutti i sistemt continut. 

Die Beitrage aus dem Gebiete der Wahrscheinlichkeit und Statistik sind: 
A.H.CopELanp, Sr., A Finite Frequency Theory of Probability; M.PRECHET, 
Interdépendance du centre et du vayon empiviques de variation de n observations 
indépendantes; F.N.FRENKIEL and J.W.FOoLLIN, jr., On Multivariate Normal 
Probability Distributions; Y.Garti et T.Consoxt, Sur la densité de probabilité du 
produit de variables aléatoives de Pearson du type III; H.GEIRINGER, On the 
Statistical Investigation of Transcendental Numbers; C. Gin, Estensioni e portata’ 
della teovia della dispersione; J.NeymaN and E.L.Scott, On the Problem of 
Expansion of Clusters of Galaxies; N.T.UzG6rEN, The Asymptotic Development 
of the Distribution of the Extreme Values of a Sample. H. Ziegler 


An Introduction to Power System Analysis. Von FREDERIK S. Roraw 
(John Wiley and Sons, New York 1953). 187 S.; $5.-. 

Das Buch behandelt in sehr moderner Form das Verhalten der Stréme, Span- 
nungen und Leistungen bei Einphasen- und Drehstromsystemen von Industrie- 
frequenz im stationdren und nichtstationaren Zustand. Insbesondere werden der 
Transformator, der Synchrongenerator und die Freileitung besprochen. Die kom- 
plexe Rechnungsweise findet durchgehend Anwendung; Zeigerdiagramme sind 
jedoch selten. Jedem Kapitel sind am Schlusse zur Einiibung des Stoffes Zahlen- 
beispiele beigegeben. 

Nach der Theorie des Transformators werden die symmetrischen Komponen- 
ten bei Drehstrom besprochen und die unsymmetrische Belastung bei Stern- und 
Dreieckschaltung berechnet. Es folgt die Anwendung auf den unsymmetrisch 
belasteten Drehstromtransformator mit Beriicksichtigung der Nullimpedanz. Die 
anschliessende Behandlung der Synchronmaschine im nichtstationéren Zustand 
ist ausserordentlich elegant und klar. Es werden zundchst die in Europa noch 
wenig bekannten Parkschen Gleichungen fiir die mit Dampferwicklungen ausge- 
riistete Schenkelpolmaschine abgeleitet. Die Parkschen Gleichungen werden nun 
auf typische Belastungsfalle angewandt, wie zum Beispiel das plétzliche Zu- 
schalten von symmetrischen und unsymmetrischen Belastungsimpedanzen und 
bei Kurzschliissen. Es folgt die Beschreibung einer «A-C Network Analyzer instal- 
lation ». Fur stark unsymmetrisch belastete Drehstromsysteme zeigen sich ferner. 
die nun erlduterten « Clarke-Komponenten » als nititzlich. Es folgt die Berechnung 
der Konstanten einer Freileitung, wobei auch die Hankelsche Funktion beniitzt 
wird. Die letzten zwei Kapitel sind der elektrischen Stabilitat von Freileitungen 
und den zugehGrigen Generatoren gewidmet, unter Verwendung des «power 
circle diagram ». 

_Der Text ist sehr klar und modern, die Ausstattung des Buches ist ausge- 
zeichnet. Es kann allen Ingenieuren, die sich mit der Energietibertragung be- 
schaftigen, aufs warmste empfohlen werden. Ed. Gevecke 


Analytische Geometrie. Von W. BLASCHKE (Birkha 

= 4 duser Verlag, Basel 1954), 

190 S. mit 67 Abb., geb. sFr. 19.60, brosch. sFr. 16.65. 4 

s Der Verfasser hat sich zum Ziele gesetzt, eine Einfiihrung in die rechnerische 
ehandlung der Geometrie zu schreiben. Wahrend die klassischen Biicher iiber 
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analytische Geometrie in der Regel nicht iiber Fragen der affinen und projektiven 
Geometrie hinausgehen, tragt BLAscHKE im vorliegenden Werk eine Fiille von 
Problemstellungen aus allen Gebieten der Geometrie zusammen, um daran die 
Methoden der analytischen Geometrie auseinanderzusetzen. Die Grundlage bildet 
der Koordinatenbegriff; an Beispielen werden speziell besprochen: die kartesi- 
schen, die homogenen, die pentazyklischen und die elliptischen Koordinaten. 
Besondere Aufmerksamkeit gilt durch das ganze Werk hindurch dem Gruppen- 
begriff und den Invarianten beziiglich einer Gruppe. Ein lingerer Abschnitt ist 
der Kreis- und Kugelgeometrie gewidmet, ein Gebiet, das dem Verfasser durch 
eigene Arbeiten besonders naheliegt. Die Darstellung der Liniengeometrie und 
der Theorie der quadratischen Formen verdient speziell hervorgehoben zu werden, 
da diese beiden Gegenstande von Problemstellungen der Mechanik (Statik und 
Kinematik des starren Kérpers) aus entwickelt werden. Das Problem der kiir- 
zesten Verbindungen auf Flachen 2. Ordnung bietet Gelegenheit zu einem Bei- 
spiel fiir die Anwendung der Infinitesimalrechnung in der Geometrie. 

Das Werk gliedert sich in folgende Abschnitte: I. Grundbegriffe, Vektoren 
und Matrizen, II. Kugeln, III. Stabe (Liniengeometrie), IV. Tragheitsmomente 
(quadratische Formen), V. Quadriken, VI. Konfokale Quadriken, VII. Formel- 
sammlung. — Die Formelsammlung beriihrt in Form eines Anhanges ohne Beweise 
auch einige weiterreichende Fragen. 

Das vorliegende Werk BLascHkKEs breitet auf relativ engem Raum eine Fiille 
von bekannteren und weniger bekannteren geometrischen Erkenntnissen aus und 
wird dadurch zu einer wahren Fundgrube. Die zahlreich eingeflochtenen histo- 
rischen Bemerkungen lassen die Namen jener wieder hervortreten (CAYLEY, 
M6sivus, CHASLES u.a.), die mit der Entwicklung der analytischen Geometrie 
aufs engste verbunden sind. Das Buch ist vortrefflich abgefasst und kann zur 
-Einfiihrung und zur Vertiefung in die analytische Geometrie bestens empfohlen 
werden. M. Jeger 


Mathieusche Funktionen und Spharoidfunktionen (Die Grundlehren der 
mathematischen Wissenschaften, Band 71). Von J. MEIXNER und F. W. SCHAFKE 
(Springer-Verlag, Berlin, Géttingen und Heidelberg 1954). 414 S., 29 Abb.; 
DM 49.-. 

Besonders in den letzten Jahren haben die Mathieuschen Funktionen und 

-Spharoidfunktionen eine ganze Reihe von bemerkenswerten Anwendungen ge- 
funden. Es ist daher zu begriissen, dass die beiden Verfasser hier ein Werk mit 
-weitgehend lehrbuchm4ssigem Charakter vorlegen, das eine abgerundete und zu- 
‘sammenfassende Ubersicht bietet iiber die wichtigsten bekannten Eigenschaften 
dieser und verwandter Funktionen. Zudem sind noch eine Anzahl bis dahin un- 
-veréffentlichter Ergebnisse aufgenommen worden. Die Darstellung ist in mathe- 
-matischer Strenge durchgefiihrt worden, und zwar vom funktionentheoretischen 
Standpunkt aus, indem sowohl fiir die Variablen als auch fiir die Parameter belie- 
‘bige komplexe Werte zugelassen werden. 

Im ersten Kapitel sind die mathematischen Grundlagen kurz dargestellt : 
‘Abriss der Theorie ganzer Funktionen endlicher Ordnung, Eigenwerttheorie fiir 
Probleme mit einem und mit zwei Parametern, dreigliedrige Rekursionen. Das 
zweite Kapitel enthalt die Theorie der Mathieuschen Funktionen, namlich Eigen- 
schaften der Lésungen, charakteristische Exponenten, Stabilitatskarte, Addi- 
tionstheoreme, Reihenentwicklungen, Integralrelationen, Asymptotik. Die ent- 

sprechende Theorie der Spharoidfunktionen wird im folgenden Kapitel darge- 
stellt. Endlich bietet das letzte Kapitel einen trefflichen Ausschnitt aus den man- 
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nigfachen Anwendungen dieser Funktionen auf mechanische, akustische und) 
elektromagnetische Probleme, wie Schwingungen, Eigenschwingungen, Abstrah- 
lung, Beugung und ferner auf wellenmechanische Probleme. — Methoden Zur 
numerischen Berechnung der Lésungen von Eigenwerten oder charakteristischen 
Exponenten sind héchstens angedeutet; solche sind ja in andern Werken bereits 
ausfiihrlich dargestellt. Den Abschluss des Buches bildet ein ziemlich vollstan-_ 
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diges Literaturverzeichnis. — Das ausgezeichnete Werk, das in verschiedener 
Hinsicht eine Liicke in der Literatur ausfiillt, wird allen, die mit diesen Funktio- 
nen zu tun haben, sehr niitzliche Dienste leisten. E. Roth-Desmeules 


Introduction to Mathematical Statistics. Von P. G. Hort, 2. Aufl. (J. 

Wiley & Sons, New York 1954), 331 S, 57-Fig.; 35,00. j 
Dieses Lehrbuch gibt trotz seines elementaren Charakters eine ausgezeichnete ~ 
Einfiihrung in die wichtigsten Methoden der modernen mathematischen Sta- 
“tistik. Vor allem sind die klaren Definitionen der benutzten statistischen Be-— 
griffe und die entsprechenden «tests» fiir die Schatzung statistischer Gréssen und — 
statistischer Hypothesen hervorzuheben. Ein umfangreiches Ubungsmaterial und © 
wichtige statistische Tabellen sorgen dafiir, dass dieses Lehrbuch zweifellos von © 
Praktikern und Studenten sehr viel benutzt wird. | W. Saxer” 


Vorlesungen iiber Differential- und Integralrechnung. Dritter Band: 
‘Integralrechnung auf dem Gebtete mehrerer Variablen. Von A. OSTROWSKI Cee 
hauser Verlag, Basel 1954). 475 S., 36 Fig.; Fr. 78.-. ees: 

Mit dem nun erschienenen dritten Band, Integralvechnung auf dem Gebiete 
’ mehrevey Variablen, ist das bereits 1945 von A. OSTROWSKI begonnene Werk Vor- 
lesungen tuber Differential- und Integralrechnung abgeschlossen worden. 4 

- Dieser Band bringt zunachst noch einige Erganzungen zur Technik des Inte- 
grierens, so die Partialbruchzerlegung und gewisse integrable Typen irrationaler 
und transzendenter Funktionen, samt einem Beweis der Transzendenz von é. } 
Anschliessend folgt schrittweise eine sorgfaltige Einfiihrung des Begriffs des | 
mehrfachen Integrals; seine Berechnung wird im nadchsten Kapitel durchgefiihrt, 
wo auch die totalen Differentiale und die Variablensubstitution zur Sprache 
kommen. Das vierte Kapitel ist den Anwendungen mehrfacher Integrale gewidmet, 
speziell den Oberflachenintegralen, den Satzen von OSTROGRADSKI und STOKES — 
diese auch in vektorieller Formulierung — ferner Volumenberechnungen und 
Schwerpunktsbestimmungen. Im folgenden Kapitel werden die verschiedenen 
Typen von einfachen uneigentlichen Integralen, deren Eigenschaften und Berech= 
nung dargelegt, wahrend das nachste Kapitel den mehrfachen uneigentlichen 
Integralen gewidmet ist; hier wird unter anderem auch die Paulische Unbe- 
stimmtheitsgleichung abgeleitet sowie die Gammafunktion und die Betafunktion 
eingefiihrt. Den Abschluss bildet ein Abschnitt iiber Fourier-Reihen und Fourier- 
Integrale. - Den einzelnen Kapiteln sind zahlreiche — im ganzen sind es tiber | 
tausend — ausgewdhlte Aufgaben beigefiigt, eine Fundgrube interessanter Pro- 
bleme. Es ist nur zu begriissen, dass der Verfasser beabsichtigt, die Lésungen der 
Aufgaben noch separat herauszugeben (das Buch soll ja auch zum Selbststudium 
_geeignet sein). ; i 

Das Buch ‘zeichnet sich durch eine grosse Sorgfalt, eine lebendige Darstellung 
sowie cine ausgesprochene Originalitat aus. Die Beweise sind wohldurchdacht, 
und das Ganze zeugt von der grossen Erfahrung des Verfassers. Damit liegt nun 
ein ausgezeichnetes Werk der Differential- und Integralrechnung vor, das wohl 
auf langere Zeit hinaus nicht so leicht libertroffen werden wird, cl 
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NEUERSCHEINUNG 


Speiser 
Die Theorie der Gruppen von endlicher Ordnung 


Mit Anwendungen auf algebraische Zahlen und Gleichungen sowie auf die Kristallographie 


Von Prof.Dr. Andreas Speiser, em. Professor an der Universitat Basel. Vierte, um einen Anhang 
und eine Farbtafel erweiterte Auflage, 284 Seiten mit 43 Figuren und einer Farbtafel. Preis 
gebunden Fr. 26.— (DM 26.-), broschiert Fr. 23.- (DM 23.-). Sammlung «Lehrbiicher und 
Monogtaphien aus dem Gebiete der exakten Wissenschaften», Mathematische Reihe, Band 22. 


Inhalt: Die Grundlagen — Normalteiler und Faktorgruppen — Abelsche Gruppen — Konjugiette © 
Untergruppen — Sylow-Gruppen und p-Gruppen — Symmettien der Ornamente — Die Kristall- 
klassen — Permutationsgruppen — Automorphismen — Monomiale Gruppen — Darstellung der 
Gruppen dutch lineare homogene Substitutionen — Gruppencharaktere — Anwendungen der 
Theorie der Gruppencharaktete — Arithmetische Untersuchungen iiber Substitutionsgruppen — 
Gruppen von gegebenem Grade — Die allgemeinen linearen homogenen Substitutionen und ihre 
Inyarianten und Kovyatrianten — Gleichungstheorie — Anhang — Namenverzeichnis — Sachver- 
zeichnis. 


In ihten elementaren Teilen besteht die Gruppentheorie aus vielleicht nicht immer organisch 
zusammenhangenden Methoden und Begtiffen, und die Gliederung des Stoffes ist schon in 
hohem Grade festgelegt. Erst mit der Theorie der Substitutionsgruppen setzt eine weittragende 
und systematische Theorie ein, die im Grunde auf zahlentheoretische Betrachtungen hinauslauft. 
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